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A Study of Network Performance and Customer 
Behavior During Direct-Distance-Dialing Call 
Attempts in the U.S.A. 


By F. P. DUFFY and R. A. MERCER 
(Manuscript received April 22, 1977) 


A survey was conducted throughout the Bell System in October 1974 
to gather detailed information about Direct-Distance-Dialing call at- 
tempts. The dispositions, setup times, and customer abandonment 
times associated with DDD attempts are discussed in detail in this 
article to provide network performance and customer behavior char- 
acteristics to network planners and administrators and to designers 
of equipment and systems which use, and interact with, the telephone 
network. It is shown that both network performance and customer 
behavior affect the call dispositions and the total call setup time; 
however, customer-dependent failures to complete account for 85 
percent of all failures, and customer-determined components of the 
call setup time make up 71 percent of the total setup time. It is found 
that traffic composition in terms of the relative mix of business and 
residential originations exerts a strong influence on call dispositions. 
Network performance affects the probability of equipment blockages 
and failures and the interval from end of dialing to receipt of a network 
response. These are both found to depend on calling distance, while 
the latter is also affected by the types of originating and terminating 
local switching. 


1. INTRODUCTION 


A complex sequence of interactions and reactions is initiated each time 
a person or machine attempts to call another person or machine via the 
switched public telephone network. In the case of local calls, the setup 
process involves station equipment, subscriber loops, at least one local 
switching office (end office) with its multitude of equipment, and per- 
haps interoffice trunks and local tandem offices with possible local al- 
ternate routing capabilities. Several local switching arrangements are 
illustrated in Fig. la and b. In the case of long distance (toll) calls, the 
switching arrangements are more complex because of a five-level 
switching hierarchy, and they are more flexible because of the extensive 
use of alternate routing.! A standard toll switching arrangement is il- 
lustrated in Fig. 1c. At one extreme of this switching arrangement, a toll 
call may encounter ten switching offices interconnected by seven final 
intertoll trunks and two toll connecting trunks; at the other extreme, a 
toll call may encounter two end offices connected by a single direct in- 
tertoll trunk. In between these extremes, a toll call may be established 
through several switching offices interconnected by a combination of 
toll connecting, final intertoll, and high-usage intertoll trunks. 

In general, a toll call may be established in several different ways be- 
tween the same two originating and terminating stations. In Fig. 1c, there 
are four separate routes which may be traversed. The toll switching al- 
gorithm establishes calls on a “link-by-link” basis. At a given office in 
the switching hierarchy, precedence is given to the trunk group which 
provides the most direct route to the terminating end office. If the trunk 
group associated with that route is busy, an alternate route is sought by 
assigning precedence to the next most direct trunk group. This process 
is repeated at each consecutive office in the switching path until a route 
is established or until the setup process cannot proceed due to a system 
blockage. 

A telephone customer making a call attempt has no perception of this 
process as it occurs, but is aware of the results: namely, the “system setup 
time” required for the network to provide an identifiable response after 
completion of dialing, and the probability that the calls are routed to the 
number dialed (‘‘system completion’). Beyond these network effects, 
calling and called customer characteristics such as the time required to 
dial a number, the abandonment behavior of calling customers during 
the call setup process, and the time required for the called customer to 
answer impact strongly upon the overall completion probability and the 
call setup time experienced by telephone customers. 

The combined influence of network performance and customer be- 
havior tends to mold customer attitudes about the telephone network 
and plays an important role in determining the amount of equipment 
needed to serve customer requests and the revenues which are realized 
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Fig. 1—Local and long-distance switching arrangements. 


from successfully completed calls. As a result, the Bell System operating 
telephone companies continually collect detailed data on network per- 
formance and, to a lesser extent, customer behavior to identify and 
correct network problems and to be able to engineer the network prop- 
erly to meet specified performance criteria. 

Bell Laboratories has gathered a sample of these data using 
a sampling plan that provides a statistically valid systemwide charac- 
terization of the call setup process as seen by the calling customer. This 
article presents an analysis of the sampled data in the expectation that 
the information on overall network call setup performance and on cus- 
tomer behavioral characteristics which affect the call setup process will 
be useful to the planners and administrators of telephone networks and 
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to designers of equipment which utilizes and systems which interact with 
the telephone network. For instance, the probability of a fast retrial 
following the occurrence of a no circuit/reorder signal or announcement, 
as presented in Table III, can be an important parameter in trunk group 
size engineering. 

Several previous studies by Bell Laboratories have dealt with the call 
setup process.2-° All but the first of these are of such different emphasis 
that the points of comparison with the current study are limited to a few 
statements on customer behavior. The Larsen and McGill study? is the 
most similar of the studies, and like the current one, is based on service 
observing data. However, the data were obtained from a number of sites 
chosen arbitrarily and not according to scientific sampling procedures. 
Hence, the representativeness of its results is unknown; and in particular, 
the confidence interval calculations do not properly account for the 
complexity of the service observing process (see Section II). That study 
emphasizes time of day and day of week trends in the completion per- 
centage and call setup time results; and while a direct comparison with 
the current study is difficult, the qualitative statements on these trends 
are borne out by the present work. The overall completion percentage 
in the Larsen and McGill study is slightly over 1 percent lower and the 
call setup time about 1 second higher than currently, but the confidence 
intervals in the current study overlap these results. 

The discussion in this article pertains to Direct-Distance-Dialing 
(DDD) toll call attempts. Those attempts which are directly dialed by 
the calling party represent approximately 88 percent of all the long- 
distance calls established within the Bell System. The sample design 
used to gather the data is described in Section II. The network and 
customer call attempt characteristics derived from this survey are pre- 
sented in Sections III and IV. 


ll. SAMPLE DESIGN 


The term “sample design” denotes the entire process of planning a 
sample survey. That process includes (7) the definition of the population 
to be studied, (ii) the adoption of a sampling plan to collect the data, and 
(iii) the derivation or selection of appropriate statistical estimation 
formulas.§ The Bell System call attempt survey described in this article 
was designed to provide a systemwide characterization of the call setup 
process utilizing data collected through an existing Bell System proce- 
dure. Therefore, considerable care was required in the design of the 
sample survey to properly account for the complex sampling scheme 
inherent in the established procedure. This design is discussed in detail 
in terms of the three features state above in the following subsections 
in order to demonstrate the statistical validity of the results. 
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2.1 Population 


A distinction must be made between the target population and the 
sampled population. The target population is the collection of elements 
about which information is desired. In some situations it may be unat- 
tractive or even impossible to sample from the target population. In those 
cases, the population definition is modified for sampling purposes, and 
the modified population is referred to as the sampled population. The 
modifications should be carefully applied so that the sampled population 
does not deviate too greatly from the target population. 

The target population for the Bell System Call Attempt Survey was 
defined to be all DDD call attempts originated by Bell System subscribers 
including attempts to toll directory assistance operators. This population 
definition was modified to take advantage of an existing data source, Dial 
Line Service Observing (DLSO). The sampled population is restricted 
to call attempts which originate from local switching offices which serve 
more than 3000 subscriber lines and consists basically of calls which are 
directed to a terminating switching office more than 25 miles from the 
originating switching office or to a different Numbering Plan Area (NPA). 
Exceptions to this general description are as follows: (7) calls to a dif- 
ferent NPA which are dailed as a seven-digit number (permissible by 
protecting codes) are excluded unless the calling distance exceeds 25 
miles, and (ii) calls to Inward Wide Area Telephone Service (INWATS) 
subscribers, or calls to toll directory assistance, which are generally dialed 
as NPA-555-1212, are included without regard to any milage restrictions. 
Operator-dialed toll call attempts and customer-dialed, operator-ser- 
viced toll call attempts are not included in this study. In general, partially 
dialed attempts are also excluded from this study because of the classi- 
fication problems posed by insufficient or incorrect digits. 


2.2 Sampling plan 


The sampling plan defines the method used to gather the survey data. 
It was intended that the call attempt data include end-to-end informa- 
tion about the setup process and the final result of that process so as to 
enable the characterization of calling and called customer behavior and 
telephone network responses. Although test calling procedures and 
human factors tests can provide valid information about some aspects 
of network responses and customer behavior respectively, only “‘live”’ 
telephone traffic can provide real-life information of this type. 

Dial Line Service Observing (DLSO) is a function performed daily at 
425 bureaus geographically located throughout the Bell System to obtain 
nationwide information about the quality of the call setup process for 
telephone traffic. An ongoing sample of actual call attempts is observed 
in over 3500 local switching offices, and the data are summarized peri- 
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Fig. 2—Geographical stratification of the contiguous 48 states. 


odically. Each observation begins at the instant the call attempt is 
originated and ends when the called party answers or when the calling 
party abandons. In the course of an observation, the call disposition and 
the timing sequence of call setup and abandonment events relative to - 
the start of the attempt are recorded. These detailed records satisfy the 
survey criterion of end-to-end information, and the general application 
of DLSO satisfies the requirement for systemwide information. In ad- 
dition, the DLSO definition of DDD call attempts agrees reasonably well 
with the target population definition for this survey. Therefore the 
sampled population definition in the previous subsection was inten- 
tionally designed to agree closely with the DLSO-DDD population in order 
to adopt DLSO as the data source for the Bell System call attempt sur- 
vey. 

The sampling plan can be briefly described as a four-stage selection 
process with stratifications imposed prior to selecting the first- and 
second-stage sample units. The first-stage units are DLSO bureaus. Prior 
to selecting the bureaus, the area of the contiguous 48 states was parti- 
tioned into 30 mutually exclusive geographical areas which are called 
primary strata. This stratification is illustrated in Fig. 2. Each primary 
stratum was designed to consist of a continuous geographical area and 
to be approximately the same size as the other primary strata with re- 
spect to the number of toll-call attempts originated from the area an- 
nually. The primary strata are much larger geographically in sparsely 
populated areas than in densely populated areas due to the second design 
criterion. This stratification guaranteed a dispersion of the sampled 
DLSO bureaus throughout the Bell System. 

Five primary strata contained only one DLSO bureau. Those bureaus 
were selected as first-stage or primary sample units with a probability 
of one. Two DLSO bureaus were selected with replacement within each 
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Table | — Bell System call attempt survey primary or first-stage 


sample 
Primary sample unit Primary Primary sample unit Primary 
DLSO bureau stratum DLSO bureau stratum 

Salem, Ma. 1 Cleveland, Oh. 17 
Providence, R.I. 1 Kast Liverpool, Oh. 17 
Boston, Ma. 2 Plymouth, Mi. 18 
Bridgeport, Ct. 3 Pontiac, Mi. 18 
New Haven, Ct. 3 Grand Rapids, Mi. 19 
Glens Falls, N.Y. 4 Springfield, Oh. 19 
White Plains, N.Y. 4 Blue Island, IL. 20 
Manhattan, N.Y. 5 East Chicago, I]. 20 
Brooklyn, N.Y. 6 Louisville, Ky. 21 
Morristown, N.J. 6 Oshkosh, Wi. 21 
Paterson, N.J. 7 Monroe, La. 22 
Newark, N.J. 8 Saint Louis, Mo. 22 
New Brunswick, N.J. 9 Fort Worth, Tx. 23 
Asbury Park, N.J. 10 Houston, Tx. 23 
Camden, N.J. 10 Salina, Ks. 24 
Lansdowne, Pa. 11 Minneapolis, Mn. 24 
Norristown, Pa. 11 Phoenix, Az. 25 
Pittsburgh, Pa. 12 Las Cruces, N.M. 25 
Fairmont, W.V. 12 Spokane, Wa. 26 
Baltimore, Md. 13 Spokane, Wa. 26 
Baltimore, Md. 13 Compton, Ca. 27 
Athens, Ga. 14 Los Angeles, Ca. 27 
Charoltte, N.C. 14 Oakland, Ca. 28 
Delray Beach, FI. 15 Oakland, Ca. 28 
Melbourne, FI. 15 Santa Ana, Ca. 29 
Columbus, Ms. 16 Santa Ana, Ca. 29 
Memphis, Tn. 16 Sacramento, Ca. 30 

Santa Cruz, Ca. 30 


of the remaining 25 primary strata. Those selections were made with 
probabilities proportional to the sizes of the bureaus as measured by the 
number of annual toll call attempts which originated in the entities 
observed by the bureaus. Because sampling with replacement was em- 
ployed to simplify the estimation formulas, four bureaus were selected 
twice. Two separate bodies of data were collected from each of these four 
bureaus. As a result of the first stage of sampling, the first-stage sample 
consists of 55 primary units which correspond to 51 distinct DLSO bu- 
reaus. The first-stage sample is listed in Table I. The analysis presented 
in this paper is based on 11,146 DDD call attempts observed by these 
bureaus during the study period. 

Each DLSO bureau had the responsibility of observing between one 
and 124 local switching entities at the time the data were collected in 
October of 1974. A total of 3521 entities was being observed by the 425 
DLSO bureaus in existence at that time for an average of approximately 
eight entities per bureau. Since the larger bureaus had a greater chance 
to be selected because the probabilities of selection were proportional 
to the sizes of the bureaus, there is an average of approximately 23 en- 
tities per DLSO bureau selected for the first-stage sample. 
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Somewhere between 20 and 100 service-observing loops were randomly 
assigned to groups of subscriber line appearances within each of these 
local switching entities. This assignment represents the second stage of 
sampling. The second-stage sample unit is a line group which is simply 
a collection of physically adjacent subscriber line appearances within 
a local switching entity. A line group varies in size from 50 to 600 sub- 
scriber line appearances depending upon the type of switching machine. 
Since the selection of line groups was made independently in each entity 
and since the entities associated with a DLSO bureau represented a 
partitioning of the bureau into mutually exclusive strata, the entities 
were treated as substrata in the sampling plan. The term “‘substratifi- 
cation” denotes a partitioning imposed within the first-stage sample 
units prior to selecting the second-stage sample. The second-stage 
sample units were randomly selected without replacement. In each line 
group belonging to the second-stage sample, the assigned service-ob- 
serving loop was connected to a randomly selected subscriber line ap- 
pearance. This selection of third-stage sample units generally was re- 
peated within each local switching entity on a weekly basis during the 
month-long study. Observations of actual DDD call attempts on the se- 
lected subscriber lines formed the fourth-stage sample. These obser- 
vations are referred to as sample elements since they comprise the 
last-stage sample and as such are members of the sampled popula- 
tion. 


2.3 Estimation procedures 


Because the Bell System call attempt survey data are not self- 
weighting, individual data items contribute differently to a system es- 
timate. The complex sampling plan structure discussed in the previous 
subsection leads to varying probabilities of selection among the indi- 
vidual data items. The individual probabilities of selection are related 
to the amount of long distance telephone traffic a data item represents. 
Appropriate traffic weights, which are basically the inverse probabilities 
of selection for the data items, are used as sampling weights in all cal- 
culations to account for the individual contribution of a data item toa 
system estimate. 

Ratio estimators of the form r = x/y, where r is an estimate of a pop- 
ulation characteristic R and x and y are estimates of the population totals 
for two random variables X and Y, were used to calculate all system 
estimates. The distribution of such an estimate is approximately normal 
for sufficiently large samples. This assumption was used to calculate 
confidence intervals for the estimates presented in this study in the form 
r +to,, where r is the ratio estimate, o, is the standard error, and t is 
the appropriate normal deviate. Since 90 percent confidence intervals 
are used consistently throughout this study, t has the value 1.65. Con- 
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fidence intervals are not given when the sample size is too small for the 
normality assumption. 

The confidence intervals indicate the precision associated with the 
population estimates. Precision is a measurement of the probable sam- 
pling error associated with a statistical estimate of a population pa- 
rameter. It gives a measure of the reliability of an estimate with respect 
to sampling errors; it does not take into consideration nonsampling er- 
rors. The 90 percent confidence intervals in this article mean that if the 
study were repeated 100 times, then 90 out of the 100 confidence inter- 
vals which could be calculated for the 100 separate estimates of a given 
population parameter would be expected to contain the true value of the 
population parameter. 


lil. NETWORK RESULTS 


Call attempt performance is characterized through the presentation 
of disposition probabilities, setup and abandonment time distributions, 
and some aspects of customer retrial behavior. The call disposition re- 
sults summarize the outcomes of DDD attempts in terms of the relative 
frequencies of completion and various reasons for failure to complete. 
The setup time results show the speed of network response during the 
provision of system signals and intercepts. They also illustrate customer 
behavior while dialing and answering telephone calls. The abandonment 
time results illustrate customer behavior in terms of releasing telephone 
facilities prior to call completion. In some cases, the facilities are released 
after an explicit system indication of failure to complete, while in other 
cases, the facilities are released following a subjective judgment of failure 
to complete or a desire to terminate the attempt prematurely. The cus- 
tomer retrial results characterize retrials which are initiated within 60 
seconds of a failure to complete. Attention is focused on each of these 
attempt characteristics individually in Sections 3.1 through 3.3, re- 
spectively. 


3.1 Call disposition 


A summary of DDD call attempt dispositions is given in Fig. 3. An es- 
timated 70.7 percent of the attempts are completed satisfactorily. Ap- 
proximately one out of every four attempts are incomplete due to a called 
customer did not answer (DA) condition, called customer line busy (BY) 
condition, calling customer dialing irregularity (CDI), or called station 
number change or disconnect (C/D) situation. These customer-dependent 
failures to complete account for 24.8 percent of all DDD attempts and 
84.6 percent of all DDD failures to complete. The second largest category 
of failures is manifested by abandonments without system responses 
(NR—no response) after dialing valid telephone numbers; they account 
for 2.6 percent of all DDD attempts and 8.9 percent of all DDD failures 
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to complete. These abandonments occur because the calling customer 
aborts prematurely or the network fails to respond properly to the re- 
quest. The third and final category of incomplete attempts consists of 
equipment blockages and failures (EB & F). Since equipment blockages 
occur because of insufficient equipment to service the request and 
equipment failures are symptomatic of faulty equipment, the respon- 
sibility for these EB & Fs clearly belongs to the telephone network. These 
network failures account for 1.9 percent of all DDD attempts and 6.5 
percent of all DDD failures to complete. 

Complete attempts are calls properly processed by the switching and 
signaling network. They represent the ideal in call disposition since they 
fulfill the calling party service requests and generate revenues for the 
telephone companies. Customer-dependent failures (DA, BY, CDI, C/D) 
also are properly processed attempts from a switching and signaling 
viewpoint. While these attempts cannot be completed for reasons beyond 
network control, the calling party is informed about the nature of the 
failure. Therefore, the results in Fig. 3 show that at least 95.5 percent 
of all DDD attempts are properly processed by the telephone network. 
This percentage is a lower bound because a substantial number of ab- 
andonments without a system response are in fact premature disconnects 
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Fig. 4— Distribution of the time from end of dialing to abandonment without a system 
response. 


on the part of the calling customers. A comparison of the distribution 
of time from end of dialing to abandonment without a system response 
in Fig. 4 with the distribution of time from end of dialing to the first 
system response in Fig. 5 illustrates that attempts abandoned within five 
seconds after dialing probably were intentionally prematurely discon- 
nected by the calling parties, since the customers could have little ex- 
pectation that a response should have been received during that interval. 
Only a few system responses are received within that time interval. This 
part of the distribution represents 30 percent of the NR conditions. At 
the other extreme of the distribution, it appears that about 14 percent 
of the NR conditions, or about 0.4 percent of all attempts, may be rea- 
sonably classified as network failures since very few system responses 
occur after 30 seconds. The remaining 56 percent of the NR distribution 
overlaps the system response distribution, and those abandonments 
cannot be easily classified as premature disconnects or network “high 
and dry” conditions. 

While at least 95.5 percent of all DDD attempts are properly processed 
from a switching and signaling viewpoint, the disposition results also 
show that 1.9 percent fail to complete due to system blockages or 
equipment problems. This EB & F estimate is a lower bound for network 
failures since some of the abandonments without a system response, 
which were discussed above, represent network “high and dry” condi- 
tions. Individual estimates for blockages and equipment failures cannot 
be calculated because a common network signal, which is a tone inter- 
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_Fig. 5—Distribution of the time from end of dialing to the first system response after 
dialing a valid telephone number. 


rupted 120 times per minute, is often issued for both failures. The tele- 
phone companies are responsible for providing and maintaining suffi- 
cient equipment to handle service requests within reasonable blockage 
and failure bounds. Since traffic loads are statistical in nature and 
equipment failures are random in nature, the system is continuously 
monitored and tuned to assure reasonable service within proper eco- 
nomic constraints. The EB & F estimate is an indication of the overall 
impact of system procedures for the provisioning and maintenance of 
equipment upon the DDD network response to customer service re- 
quests. 

The reciprocal of the completion ratio gives the average number of 
attempts per successfully completed call. Using the DDD completion ratio 
of 0.707, this average is computed to be 1.4 attempts per successful 
completion. Substantial improvement in this average can only come 
about by a reduction in failures which are dependent upon calling and 
called customer behavior, since such failures account for about 85 percent 
of all DDD failures to complete. The occurrence of such failures can be 
reduced through the mutual cooperation of telephone companies and 
telephone customers. Reasonable service agreements in terms of the 
types and amounts of equipment required for individual customer 
communication requirements in conjunction with the application of new 
services made possible through the advent of Electronic Switching 
Systems (ESS) such as call forwarding and call waiting can produce a 
substantial impact. 
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Fig. 6—Call attempt time for incomplete DDD attempts. 


Complete eradication of all network failures would only decrease the 
average number of attempts per completed call by about 5 percent. 
Reduction of network failures is an area of continuing interest, and the 
historical record of EB & F reduction shows continuous improvement. 
Whenever practical, new technological advances continue to reduce 
network failures. For example, modern ESS offices will make a second 
attempt automatically if they are able to detect an irregularity in es- 
tablishing the next link in a call. 


3.2 Call attempt time 


Call attempt time is defined as the nonconversational period which 
commences with calling station off-hook and terminates with called 
station answer for completed calls and with calling station return to 
on-hook for failures to complete. This usage of network facilities during 
the call setup process is absorbed as network overhead since revenues 
are not directly obtained. The average call attempt time with accom- 
panying 90 percent confidence interval is 45.1 + 1.8 seconds for incom- 
plete attempts and 32.8 + 0.7 seconds for completed calls. Call attempt 
times for incomplete attempts are not only longer on the average than 
those for completed calls, they are also more variable as is seen by the 
standard deviations of 24.7 and 11.0 seconds respectively. The distri- 
butions of call attempt time for incomplete attempts and completed calls 
are presented in Figs. 6 and 7 respectively. The overhead represented 
by these call attempt times amounts to an average of 50.8 seconds per 
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Fig. 7—Call attempt time (nonconversational) for successfully completed DDD at- 
tempts. 


completed call, since 1.4 attempts are initiated for each successful at- 
tempt (32.8 + 0.4 X 45.1). If the time from off-hook to start of dialing 
plus start of dialing to end of dialing (13.8 sec) is excluded on the basis 
that it impacts mainly upon local facilities, the toll network overhead 
is 31.5 seconds per completed call (50.8 — 1.4 X 138.8). 

Further insight into call attempt time is gained by considering the 
customer-controlled and network-controlled components. The com- 
ponents of the call attempt time interval are schematically illustrated 
in Fig. 8, which also contains the average holding time associated with 
each component. A more detailed statistical characterization, which 
includes the mean, standard deviation, and 10, 50, and 90 percent points 
of the cumulative distribution function, is given for each component in 
Table II. The calling customer components of call attempt time include 
the intervals from off-hook to start of dialing and from start of dialing 
to end of dialing, which are common to all attempts, and the intervals 
from the beginning of a system response to disconnect, from the begin- 
ning of an answer at the wrong station to disconnect, and from the end 
of dialing to abandonment without a system response for incomplete 
attempts. Called customer behavior is directly illustrated in two call 
attempt time components; namely, the intervals from the start of ringing 
to answer and from the end of dialing to answer with no audible ringing. 
These customer controlled components of call attempt time are discussed 
in Section 3.2.1 below. The network controlled components include the 
intervals from end of dialing to any of the various network signals, an- 
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Fig. 8—Average DDD call attempt setup and abandonment times (seconds) with 90 
percent confidence intervals. 


nouncements, or manual intercepts, and they are discussed in Section 
3.2.2 below. 


3.2.1 Customer-controlled components of call attempt time 


The intervals from off-hook to start of dialing and from start of dial- 
ing to end of dialing constitute the service request period of a call at- 
tempt. On the average, a request consumes 13.8 seconds of the overhead 
time per attempt. The type of dialing equipment at the originating 
location substantially impacts upon this service request time, since the 
average dialing time per DDD attempt is 13.7 seconds at rotary dial sta- 
tions and 7.0 seconds at TOUCH-TONE*® dialing stations. The interval 
from off-hook to start of dialing is also slightly longer for rotary dial 
customers than for TOUCH-TONE customers. The end result of these 
differences due to dialing equipment are rotary-dial service request times 
on the average of 15.5 seconds which are almost twice as long as the 
TOUCH-TONE service request times, which average 8.6 seconds. 

In addition to influencing the service request phase of the call attempt 
process, calling customers also influence the speed of abandonment as- 
sociated with incomplete or incorrectly completed attempts. Customers 
generally abandon faster when there is a positive network indication of 
failure to complete than when a subjective judgment of failure is re- 
quired. This finding is borne out by a comparison of the abandonment 
times following network indications of line busy (BY), no circuit or re- 
corder (NC/RO), and no-such-number (NSN) conditions with the aban- 
donment times associated with disconnects after receipt of audible 
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Table II—DDD call setup and abandonment time statistics 








Cumulative 
distribution 
percentage points 
(sec) 
Call setup or abandonment interval Mean (sec) Std. dev. (sec) 10% 50% 90% 
Off-hook to start of dialing 1.7+0.1 2.0 1 2 4 
Start of dialing to end of dialing 12.1+40.3 5.1 6 12 18 
End of dialing to ring before answer or disconnect 10.9 + 0.5 5.0 5 11 17 
End of dialing to answer without a ring signal 11.1+0.9 5.3 5 10 17 
Start of ringing to answer 8.5 + 0.3 7.2 3 6 15 
Start of ringing to disconnect without an answer 38.1 + 2.5 20.7 18 33 59 
Answer to disconnect for wrong numbers 18.7 8.9 — — = 
End of dialing to busy (60 IPM) signal 10.5 + 0.9 5.2 5 10 17 
Start of busy to disconnect 46+ 0.3 3.7 2 3 8 
End of dialing to no circuit/reorder (120 IPM) signal 7.4 + 2.0 7A 0 6 15 
Start of no circuit/reorder to disconnect 3.6 + 0.4 2.3 2 3 6 
End of dialing to no-such-number (NSN) tone 1.2 1.4 — _ — 
Start of NSN tone to disconnect 2.6 1.7 —_ — — 
End of dialing to ring prior to NSN Announcement 5.3 4,1 — _ _ 
Start of ringing to NSN announcement 7.1 4.9 
End of dialing to NSN announcement without a ring signal 7.3416 5.4 1 6 11 
Start of NSN announcement to disconnect 76418 10.4 2 5 13 
End of dialing to ring prior to intercept for number change or disconnect (INT) 10.9 + 1.1 4.6 6 10 19 
Start of ringing to INT 8.8 + 1.8 6.3 3 9 15 
End of dialing to INT without a ring signal 13.9 4.4 — — — 
Start of INT to disconnect 24.7 + 14.2 24.6 4 15 81 
End of dialing to customer abandonment without a system response 15.5 + 2.8 17.8 1 12 35 


ringing for DA conditions and with disconnects without a system response 
for network “high and dry” conditions. The reasons for the relatively 
slow abandonment times associated with intercepts for number changes 
or disconnects (C/D) will be explained later. The BY, NC/RO, and NSN 
indications immediately signal failure and encourage the calling party 
to abandon; while in other cases, the calling party abandons only after 
judging that the called party had sufficient time to answer or that the 
network has had sufficient time to respond. 

Calling customers also appear to differentiate between the various 
definitive network indications for failure. The NSN tone, which contin- 
uously varies in frequency, is quite different from the BY and NC/RO 
signals. It consistently encourages rather fast abandonments as can be 
seen by the estimates for the mean and standard deviation in Table II. 
Hither the characteristics of the NSN tone itself or the unfamiliarity of 
the tone may cause this behavior. More surprisingly, however, customers 
also seem to distinguish between BY and NC/RO signals since both the 
mean and the standard deviations are smaller for abandonment times 
after NC/RO signals compared with those after BY signals. The main 
distinction between these signals is that the BY tone is interrupted 60 
times per minute and the NC/RO tone is interrupted 120 times per min- 
ute. 

Abandonments which follow definitive network signals for failure to 
complete occur faster than those which follow recorded, automatically 
composed (computer-generated) or manual intercepts. NSN an- 
nouncements and C/D intercepts convey verbal messages, and this 
transfer of information requires additional time. While all NSN an- 
nouncements were automated, one third of the C/D intercepts were 
manual. This partially explains the slower C/D abandonment time, since 
the average abandonment time for the manual C/D intercepts is ap- 
proximately 10 seconds slower than for the automated C/D intercepts. 
In addition, the average abandonment time following automated C/D 
intercepts is almost twice as slow as the average for automated NSN 
announcements. In some instances an automated C/D intercept is fol- 
lowed by an “operator intercept cut-through” to provide the calling party 
with an opportunity to discuss the problem with a person. In such cases 
the time to abandon could well be even greater than for an immediate 
manual intercept. The occurrence of such cut-throughs are not recorded 
in the data base; however, if it occurred during a significant number of 
the automated C/D intercepts, it might well account for slower aban- 
donment times. That some such cut-throughs occurred is suggested by 
the presence of several outlying values in the automated C/D abandon- 
ment time data. These outlying values along with several long holding 
times for the manual C/D intercepts inflate the confidence interval 
substantially and result in an average overall C/D abandonment time of 
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24.7 seconds which is quite larger than the median of 15 seconds. Ex- 
cluding the automated C/D intercepts which appear to be followed by 
operator intercept cut-throughs, the average abandonment time for 
automated C/D intercepts is still longer than for NSN announcements. 
Since the standard messages which are conveyed are similar in length, 
this difference must be attributed to customer reaction. 

While some type of system response is encountered for 97.4 percent 
of all DDD attempts, the remaining 2.6 percent are abandoned after 
dialing without a system response. They are abandoned because the 
calling party no longer desires to complete the call, is interrupted by some 
event, or perceives that the network has failed to respond properly to 
the dialed digits. Abandonments without a system response were dis- 
cussed in Section 3.1 above. The conclusions were (1) about 30 percent 
clearly appear to abort prematurely, (ii) about 14 percent clearly appear 
to be network failures, and (iii) about 56 percent of the distribution for 
these abandonment times overlaps the distribution for the times from 
end of dialing to first system response. 

The final calling customer abandonment characteristic to be discussed 
in this article pertains to the 19 wrong numbers which were observed 
during the survey. The average time from answer at a wrong station to 
abandonment is 18.7 seconds; however, this estimate is based upon a very 
small number of observations. The confidence interval for this estimate 
is omitted in Fig. 8 and Table II because the small number of observa- 
tions negates the meaningfulness of the confidence interval calcula- 
tion. 

The large differences in average abandonment time among the various 
reasons for abandoning cause the high variability seen earlier for the call 
attempt times associated with incomplete attempts. The longest average 
call attempt time is associated with the single most frequent reason for 
failure to complete, called customer did not answer. The average call 
attempt time for DA conditions is 62.8 seconds, of which 38.1 seconds 
are consumed by the calling party listening to ringing. An average of 
between six and seven rings are received before the calling party aban- 
dons. Since over 99 percent of all DDD answers occur prior to the 38.1 
second average abandonment time, it is meaningful to investigate pos- 
sible benefits associated with an effort to change customer behavior 
through educational programs. Not only may it be beneficial to dis- 
courage very long holding times, it may also be beneficial to discourage 
very short holding times to give called customers more opportunity to 
respond. 

At the other end of the abandonment time spectrum, abandonments 
following a busy signal occur rather fast. However, 30 percent of the 
customers who encounter a BY signal remain on the line over 5 seconds 
following the initial receipt of the tone. Since 10.1 percent of all DDD 
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attempts result in busies, this excessive holding time represents a drain 
on network facilities without any possibility of completion. 

Focusing attention at the terminating end of a call, one aspect of called 
customer behavior is provided by the following answer time results: (7) 
the average time from end of dialing to called station answer with audible 
ringing observed at the calling station is 19.4 seconds, and (ii) the average 
time from end of dialing to called station answer without the receipt of 
audible ringing at the calling station is 11.1 seconds. The substantial 
difference between estimates is mainly attributed to toll directory as- 
sistance calls. The first category, which represents 96.1 percent of all DDD 
calls (messages), contains relatively few calls to toll directory assistance 
(5 percent), while slightly over two-thirds of the calls in the second cat- 
egory are to toll directory assistance. The average time to answer for toll 
directory assistance operators is 7.4 seconds when ringing is heard at the 
calling station and 1.1 seconds when it is not. This suggests that in the 
former case the operators are busy when the call appears, and that in the 
latter case an idle operator answers almost immediately; therefore, 
ringing is not observed at the calling station. Customer answers take an 
average of 8.5 seconds when ringing is observed at the calling station and 
4.2 seconds when it is not. These results along with the composition of 
the two categories described above account for about two-thirds of the 
difference between the 19.4 and 11.1 second estimates. The remaining 
third is explained by differences in connect times for toll directory as- 
sistance calls and other DDD calls. The average time to establish a con- 
nection from the originating local switching office to the called station 
is 3.8 seconds shorter for toll directory assistance calls than for other 
DDD calls. . 

About 71 percent of the average call attempt time is under customer 
control during the service request or answer and abandonment phases 
of the call attempt process. The service request phase represents 33 
percent of the call attempt time, and the answer and abandonment phase 
represents 38 percent. The remaining 29 percent is under network control 
for the purpose of establishing a physical connection between the calling 
and called stations. 


3.2.2 Network controlled components of call attempt time 


Upon completion of dialing, the network attempts to route a call to 
the called station by establishing transmission and signaling paths be- 
tween a number of local and toll switching offices. A successful connec- 
tion may be established, whereupon audible ringing is normally returned 
to the calling station, or the network may be unable to reach the called 
station, whereupon a specific network indication of failure to complete 
is normally returned to the calling station. The first system response is 
defined as the first network signal, network intercept, or called station 
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answer without audible ringing at the calling station which follows the 
end of dialing. The distribution of these response times is given in Fig. 
5 and has an average time to respond of 10.9 + 0.5 seconds with a stan- 
dard deviation of 5.3 seconds. In Fig. 8 and Table II, summary statistics 
are listed for the individual system responses which include answer 
without audible ringing at the calling station, audible ringing, busy signal, 
no circuit or reorder signal, no-such-number tone or announcement, and 
intercept for number change or disconnect. 

The average times from end of dialing to ring prior to answer or dis- 
connect, busy, NC/RO, and NSN tone are 10.9, 10.5, 7.4, and 1.2 seconds 
respectively. The averages for ring and busy are relatively long because 
a connection must be established to the far-end local switching office 
and the called customer line must be tested before a ring or busy signal 
is returned to'the calling station. A NC/RO signal is returned to the calling 
station when an equipment blockage or failure occurs. The average time 
from end of dialing to receipt of the NC/RO signal is shorter than the ring 
and busy averages because the signal is applied at the office in which the 
EB & F is encountered. That office may be the originating local switching 
office, an intermediate toll switching office, or the terminating local 
switching office. The average time from end of dialing to a NSN tone is 
the shortest of the four because the tone generally can be applied at the 
originating local switching office or the first toll switching office en- 
countered in the switching path. 

The two remaining responses which qualify as a first system response 
after dialing are recorded, composed, or manual intercepts for no- 
such-number conditions and number change or disconnect conditions. 
The NSN intercept is encountered near the originating end of a con- 
nection in a manner similar to the NSN tone. The C/D intercept occurs 
at the far-end of a connection. As a result of these distinctions, it is clear 
why the average time from end of dialing to C/D intercepts is almost twice 
as long as the average for NSN intercepts. Both types of intercept often 
are preceded by audible ringing at the calling station. Note that the av- 
erage time from end of dialing to ringing prior to a NSN intercept is only 
5.3 seconds. This clearly reflects the fact that NSN intercepts occur 
towards the near end of a connection. The average time from end of 
dialing to ringing prior to a C/D intercept is 10.9 seconds which is and 
should be in agreement with the average for ringing prior to answer or 
disconnect. 


3.3 Fast retrials 


In the event that an attempt is not successfully completed, the calling 
party must decide to permanently or temporarily abandon or to try 
again. The retrial probability and time to retry are completely deter- 
mined by calling customer behavior under the influence of the various 
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Fig. 9—Time to retry for retrials which occur within 60 seconds of a DDD failure to 
complete. 


network indications for failure to complete. Retrial disposition is a 
function of both calling and called customer behavior, and to a lesser 
degree network behavior such as the persistence of NC/RO conditions. 

An estimated 19.1 + 2.5 percent of all DDD attempts which fail to 
complete are retried within 60 seconds of the failure, and 34.3 + 5.4 
percent of those retrials are successful. Over 50 percent of these fast 
retrials are initiated within two seconds of the abandonment time for 
the preceding failure. This time to retry is characterized by the distri- 
bution given in Fig. 9 for which the average is 6.7 + 1.3 seconds and the 
standard deviation is 10.7 seconds. The substantial difference between 
the average and median times is caused by the high degree of positive 
skewness associated with the distribution. 

Fast retrial statistics are summarized in Table III for the most frequent 
types of failures to complete. The fast retrial rates and the contribution 
to the overall fast retrial phenomenon for individual types of failures 
to complete are listed in the top part of the table. The major conditional 
disposition probabilities associated with retrials after failure to complete 
are given in the lower part of the table. Since these probabilities were 
derived from rather small samples, they are only intended to provide 
an indication of what happens to fast retrails, and are quoted without 
confidence limits. The sample size for the DA, BY, abandoned without 
a system response, NC/RO, and NSN results are 95, 183, 166, 66, and 85 
attempts respectively. 

Almost one-third of all DDD fast retrials follow a failure due to called 
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Table Ill — Retrial characteristics for retrials occurring with 60 
seconds of an unsuccessful attempt 


Unsuccessful Percent which retry* Percent of retrials* 
call attempt dispositions (within 60 seconds) (within 60 seconds) 
Called party did not answer (DA) 5.24£1.6 11.7 + 3.9 
Called party busy (By) 17.3 + 3.8 30.9 + 5.1 
Calling party abandoned without 53.44 6.5 24.4 + 4.2 
a system response (ABAN) 
No circuit/reorder condition 41.7 + 13.8 11.2 + 3.5 
(NC/RO) 
No-Such-Number condition 55.5 + 8.8 14.2 + 4.7 
(NSN) 
Unsuccessful Conditional disposition probabilities 
call attempt disposition (60 second retrials) 
Called party did not answer (DA) Complete: 0.23 DA:0.18 No order: 0.57 
Called party busy (By) Complete: 0.20 By: 0.61 


Calling party abandoned without Complete: 0.46 No order: 0.19 
a system response (ABAN) 


No circuit/reorder condition Complete: 0.49 NC/Ro: 0.40 
(NC/RO) 

No-such-number condition Complete: 0.31 NSN:0.23 No order: 0.24 
(NSN) 


Percent which retry—The percent of call attempts which fail to complete and are tried 
again within 60 seconds for the specific failures listed. 

Percent of retrials—The percent of all retrials occurring within 60 seconds of a failure 
to complete which are attributable to the specific failures listed. 


customer busy. These retrials are very likely to result in a second BY 
condition as is evidenced by the conditional probability for a BY following 
a BY of 0.61. The conditional probability for successful completion after 
a BY condition is only 0.20. While the fast retrial rate for called customer 
did not answer conditions is relatively low, such retrials account for 11.7 
percent of all fast retrials because DA conditions occur relatively fre- 
quently on first attempts. Over half of these retrials which follow DA 
conditions are no orders, i.e., the calling party goes off-hook and back 
on-hook without dialing any digits or dials a partial number and returns 
to the on-hook state within 10 seconds after the end of dialing. The 
conditional probability for successful completion following a DA con- 
ditions is 0.23. 

Fast retrial rates are relatively high for attempts which are abandoned 
without a system response or encounter NC/RO indications. The condi- 
tional probabilities for successful completion are also relatively high at 
0.46 and 0.49 respectively. While less than 10 percent of the retrials after 
an abandonment without a system response also are abandoned without 
a system response, the reoccurrence of NC/RO having already received 
a NC/RO indication is very likely as can be seen by the conditional 
probability of 0.40 in Table III. 

The fast retrial rate is also relatively high for attempts encountering 
NSN intercepts. The conditional disposition probabilities for fast retrial 
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following NSN intercepts present a mixed picture. These probabilities 
are 0.31, 0.23, and 0.24 for successful completion, NSN intercepts, and 
no orders respectively. The latter two probabilities indicate substantial 
confusion on the part of the calling party, since the first suggests the 
customer again dials the incorrect number while the second indicates 
the customer is too uncertain to continue dialing. 

Fast retrials comprise only a small segment of the complete retrial 
process, since hours and even days may elapse between a failure to 
complete and a legitimate retrial. The data collection scheme for the 
DLSO call attempt survey was not designed to capture information on 
other than immediate retrials; however, a systemwide survey to char- 
acterize the complete retrial process is currently in progress at Bell 
Laboratories. 


IV. CORRELATION OF NETWORK RESULTS WITH SEVERAL ATTEMPT 
CHARACTERISTICS 

Call disposition is strongly influenced by certain attempt character- 
istics such as class of subscriber service and calling distance. Call setup 
time also is influenced by these characteristics and by the type of 
switching which is encountered during the setup process. While the class 
of subscriber service effects are totally dependent upon customer in- 
fluences, the calling distance and switching machine effects are also 
dependent upon network influence. Call disposition characteristics vary 
by day of week and time of day; however, the major causative influence 
for these variations is the relative traffic composition by class of sub- 
scriber service. The roles which class of subscriber service, calling dis- 
tance, and type of switching play in the call attempt process, are explored 
in the following subsections. 


4.1 Class of subscriber service 


Service observing is restricted to originating attempts. While the 
originating class of service is known for most observations, the termi- 
nating class of service is unknown for all observations except for calls 
to directory assistance and INWATS calls. The survey results discussed 
in this subsection are predicated on originating class of service with the 
main service distinction being residential versus business. Substantial 
differences in customer behavior are clearly shown by the DDD com- 
pletion rates of 76.6 +2.2 and 66.9 £2.5 percent for business and resi- 
dential customers respectively. The primary cause of this 10 percentage 
point difference is attributed to called customer behavior in terms of did 
not answers and busies, since the business DA rate of 8.9 £1.6 percent 
is approximately half the residential rate of 16.4 +2.3 percent, and the 
business BY rate of 8.6 +1.4 percent is approximately three-fourths the 
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residential rate of 11.38 +2.1 percent. Of the remaining dispositions, only 
customer abandonments without a system response border upon dif- 
fering significantly. These estimates are 2.7 +0.7 and 1.7 +0.5 percent 
for business and residential attempts. Although the business abandon- 
ment rate is higher, business customers are as patient as residential 
customers in the sense that the median holding time prior to abandoning 
without a system response is 12 seconds in both cases. In addition, the 
distributions for the time from end of dialing to first system response 
are basically the same for business and residential traffic. These sys- 
temwide abandonment time and response time similarities for business 
and residential traffic seem to contradict the difference between the 
business and residential abandonment rates. The cause of this apparent 
contradiction is not known. It may be due to some characteristic such 
as routing for which the information is not available, or it may be due 
to sampling error. 

Traffic composition in terms of the relative percentages of business 
and residential originated attempts exerts a strong influence upon overall 
DDD dispostion characteristics. In the following paragraphs it is shown 
that variations ‘in the mixture of business and residential originations 
are responsible for many day of week, time of day, and calling distance 
trends associated with DDD call dispositions. In general, completion 
increases as the percentage of business originations increase because the 
occurrence of called station did not answer and busy conditions de- 
cline. 

Weekend traffic is primarily residential in nature. Weekday traffic 
is almost evenly split between business and residential originations. As 
a result of this difference in traffic composition, the weekend completion 
rate of 69.1 + 4.8 percent is lower than the weekday completion rate of 
71.8 + 1.8 percent. A high rate of DA conditions, which is typical of res- 
idential traffic, is the primary cause of the lower weekend completion 
rate. 

In addition to this weekday versus weekend variation, traffic com- 
position varies by hour of day. The percentage of business-originated 
DDD attempts by hour of day for weekday traffic is illustrated in Fig. 10 
for the hours of 8:00 A.M. to 10:00 P.M. The business-residential com- 
position is fairly constant during the morning hours of 8:00 A.M. to 12:00 
noon. During the noon hour the business originations drop by about 10 
percentage points. In the afternoon business originations initially jump 
to about 58 percent of the total traffic and then gradually decrease to 
about 49 percent of the traffic by 4:00 P.M. At this time a transitional 
period begins during which the business-residential mixture changes 
dramatically. By 7:00 P.M. the percentage of business originations is 
down to about 16 percent of the total traffic. The evening traffic from 
7:00 to 10:00 P.M. consists of almost completely residential originated 
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Fig. 10—Percent business-originated attempts with 90 percent confidence interval by 
hour of day. 


attempts. 

The variations in business-residential traffic composition are clearly 
evident in the disposition results for the five periods of day in italics 
above. The estimates for percent complete, which are listed in Table IV, 
increase and decrease along with the percentage estimates for business- 
originated attempts which are also listed in Table IV. These changes in 
completion are accompanied by corresponding changes in the opposite 
direction for the composite estimates of DA + BY conditions. While 
completion increases when the percent of business-originated attempts 
increases due to fewer DA and BY conditions, the increase in completion 
is not as great as the decrease in DA + BY conditions because the calling 
customer abandonment rate also increases with business-originated 
attempts. Considering the previous description of business and resi- 
dential call disposition characteristics, it is clear that the time of day 
variations in call disposition are essentially caused by variations in the 
business-residential traffic composition throughout the day. 

Business-residential traffic composition also varies with calling dis- 
tance, which is defined as the airline distance between the originating 
and terminating local switching offices. While the majority (56 percent) 
of DDD attempts with calling distances between 26 and 400 miles orig- 
inates from residential stations, the majority (also 56 percent) of the 
attempts with longer calling distances originates from business stations. 
This change in traffic composition contributes to several call disposition 
trends related to calling distances. Positive correlation of completion 
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Table IV — DDD call disposition estimates with 90 percent confidence intervals for five periods of the day for 
weekday traffic 


Disposition 8-12 A.M. 

Complete to desired station 74.14 2.3 

Called station did not answer 8.7+1.5 

Called station busy 96+ 1.4 

Calling customer abandoned without a 3.4 + 0.8 
system response 

Percent business traffic 53.6 + 6.2 


Percentage estimates by period of day 


12-1 P.M. 1-4 P.M. 4-7 P.M. 

67.5 + 8.3 720 4:2.9 66.5 + 5.0 

15.2 + 6.3 11.6 + 2.1 13.6 + 2.6 

11.3 + 6.2 9.4+ 1.6 15.1 + 5.5 
1.7 Zio Oe 1.4 

43.5 + 8.3 55.1 + 6.2 33.2 + 7.1 


with calling distance comes about because the business completion rate 
is higher than the residential rate and the percentage of business origi- 
nations increases as calling distance gets longer. In a similar manner, 
negative correlation of the DA + BY rate with calling distance comes 
about because the business DA + BY rate is lower than the residential 
rate and the percentage of business originations increases as calling 
distance gets longer. The abandonment rate without a system response 
also is positively correlated with calling distance. In this case, the busi- 
ness abandonment rate of 2.7 40.7 percent appears to be invariant with 
calling distance, while the residential abandonment rates of 1.5 +0.6 
percent for attempts with calling distances between 26 and 400 miles 
and 3.0 +2.1 percent for those with calling distances greater than 400 
miles appear to be positively correlated with calling distance. These 
trends, along with additional call disposition trends which are related 
to calling distances but are not influenced by class of subscriber service, 
are discussed more thoroughly in Section 4.2. 

Class of subscriber service plays only a minor role in determining call 
setup and abandonment times, since statistically, call setup and aban- 
donment time results for business- and residential-originated traffic are 
similar with the exception of customer dialing times. The average resi- 
dential rotary dialing time is 1.5 seconds longer than the average time 
for business customers. The average residential TOUCH-TONE dialing 
time is 1.3 seconds longer than the average time for business customers. 
Residential and business customers originate almost the same percentage 
of Home company, Home Numbering Plan Area (H-HNPA) attempts. 
Therefore, the differences in dialing times cannot be attributed to calling 
patterns in which one group is required to dial NPA codes more fre- 
quently than the other. Most likely, business customers have shorter 
dialing times because they use the telephone more frequently or are more 
familiar with the numbers they call. The only clear correlations of call 
setup and abandonment times with day of week or time of day once again 
is for customer dialing time for which the average dialing time on 
weekends is longer than the average for weekdays, and the average 
dialing time in the evening is longer than the average during the re- 
mainder of the day. These findings are in total agreement with the dif- 
ferences between business and residential dialing times and the traffic 
composition characteristics which were discussed previously. Namely, 
dialing times tend to be longer when the traffic is dominated by resi- 
dential originations and shorter when the traffic is dominated by busi- 
ness originations. 


4.2 Calling distance 


Beyond the customer influences discussed above, there are additional 
network effects upon the call attempt process that appear in the calling 
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Table V—DpD call disposition percentage estimates with 90 percent 
confidence intervals as a function of calling distance 


Disposition 26-400 miles >400 miles 
Complete to desired station 69.8 + 2.0 72.0 + 3.9 
Called station did not answer 13.9 + 1.6 10.9 + 2.8 
Called station busy 9841.3 9.0 + 2.3 
Calling customer abandoned 2.3 + 0.5 40+ 2.1 

without a system response 
No circuit/reorder condition 1.4 2.5 
No-such-number condition 1.740.7 0.5 
Other failures to complete 1240.5 1.1 


distance discussion in this section and in the type of switch discussion 
in Section 4.3 below. 

As stated previously, calling distance is defined as the airline distance 
between the originating and terminating local switching offices. Gen- 
erally, as this distance increases, the call setup process becomes more 
complex due to additional toll switching, intraoffice processing, and 
interoffice trunking. This added setup complexity has considerable 
impact upon call disposition and setup times. It does not have much 
impact upon customer abandonment times. 

Direct-Distance-Dialing attempts with calling distances of 25 miles 
or less are excluded from the following analyses because many such at- 
tempts do not belong to the sampled population for the DLSO call at- 
tempt survey (see Section 2.1). INWATS attempts are also excluded from 
the following analyses because the calling distances are not known for 
the intrastate, INWATS attempts. All remaining attempts to both cus- 
tomer stations and toll directory assistance operators are included in 
the calling distance analyses which follow. 

Disposition results for attempts with calling distances between 26 and 
400 miles and for attempts with calling distances greater than 400 miles 
are listed in Table V. As stated in the previous subsection, the completion 
and abandonment without a system response rates both increase and 
the DA + BY rate decreases as the calling distance becomes longer. The 
completion and DA + BY trends are primarily attributed to changes in 
business-residential traffic composition. The abandonment without a 
system response trend is attributed to a dependence of residential ab- 
andonments upon calling distance. 

The occurrence of no circuit/reorder conditions increases as the calling 
distance gets longer. This trend is caused by the added complexity to 
the call setup process. The increase in toll switching creates more op- 
portunities for failure to obtain the required facilities and for equipment 
irregularities. 

The occurrence of no-such-number conditions decreases as calling 
distance increases due to the substantial difference between the 3.3 +1.7 
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percent intra-NPA NSN rate and the 0.7 +0.3 inter-NPA rate. While over 
three-fourths of the intra-NPA NSN conditions are due to the omission 
of the toll access code 1, which is required in many areas of the country, 
only about one-third of the inter-NPA NSN conditions are due to that 
reason. Apparently the presence of an NPA code in the called number 
reminds the calling party to dial the access code. 

The impact of calling distance upon setup time is shown in Table VI. 
Average dialing times are greater for attempts with long calling distances 
because such attempts more frequently require an NPA code to establish 
a call. The interval from end of dialing to a ring or busy signal also in- 
creases as calling distances get longer because additional switching tends 
to slow down the setup time. The estimates for end of dialing to a NC/RO 
signal do not indicate any dependence upon calling distance. Because 
NC/RO conditions can occur at any intermediate office within the 
switching path and ring or busy are always returned from the terminating 
local switching office, the average time from end of dialing to NC/RO is 
shorter than to ring or busy, while the standard deviation is wider. 

The last four intervals listed in Table VI represent customer aban- 
donment times. The second and third of these intervals are for aban- 
donments after positive indications of failure to complete; namely, busy 
and NC/RO signals. The results show that there is no correlation between 
calling distance and the time to disconnect after hearing a NC/RO signal. 
The results also indicate that there may be a positive correlation between 
calling distance and the time to disconnect after hearing a busy signal; 
however, the evidence for such a correlation is rather weak. The two 
remaining abandonment times for disconnect after receipt of audible 
ringing with no answer and disconnect without a system response are 
dependent upon calling customer judgments that the attempts would 
not complete satisfactorily. The time to disconnect after receipt of au- 
dible ringing without an answer is not correlated with calling distance. 
The time to disconnect without a system response may be positively 
correlated with calling distance; however, once again, the evidence for 
such a correlation is rather weak. 


4.3 Local switching and call setup time 


The average setup time from end of dialing to ring or busy varies 
substantially for attempts classified by the types of local switching. The 
average setup times with accompanying 90 percent confidence intervals 
and standard deviations are listed in Table VII for originating and ter- 
minating local switching classifications. Estimates for attempts which 
originate from rotary dial and TOUCH-TONE dialing stations are given 
separately for the originating local switching classifications. 

Call attempts which originate from TOUCH-TONE dialing stations 
served by step-by-step (SXS) switching machines have a significantly 


NETWORK PERFORMANCE AND CUSTOMER BEHAVIOR 29 


8Z6L AYVANVE “IWNYNOP IWOINHOSL WSLSAS Tad SHL 0€& 


Table ViI—DDD call attempt setup and abandonment time statistics as a function of calling distance 


26-400 miles >400 miles 
Mean Std. dev. Mean Std. dev. 

Setup or abandonment interval (sec) (sec) (sec) (sec) 
Rotary dialing time 13.1+0.4 4.6 14.4 + 0.3 4,3 
TOUCH-TONE® dialing time 6.9 + 0.4 3.3 8.5 + 1.0 4.4 
End dialing to ring or busy 11.1+0.6 5.2 12.5+0.4 4.2 
End dialing to no circuit/reorder 8.6 + 3.4 8.7 6.9 + 3.1 5.6 
Start of audible ring to disconnect 39.9 + 3.3 20.7 39.3 + 4.2 21.1 

without an answer 

Start of busy to disconnect 4440.4 3.5 5.8 + 1.6 4.6 
Start of no circuit/reorder to disconnect 3.8 + 0.4 2.5 3.7 + 0.9 1.7 
End dialing to abandonment without a 14.84 3.1 13.8 20.2 + 7.3 27.9 


system response 


Table VII— DDD call attempt statistics for the time from end of 
dialing to ring or busy as a function of originating and terminating 
local switching 


Originating Rotary dial customers TOUCH-TONE® customers 

switch Mean (sec) Std. dev. (sec) Mean (sec) Std. dev. (sec) 
PAN 13.6 + 3.0 7.4 a = 

SXS 12.6+41.6 5.8 16.5 + 1.7 4.7 

5XB 0.14 0.4 4.4 10.8 + 0.6 4.6 

1XB 10.8 + 0.8 5.8 12.3+1.9 4.6 

ESS 9.2 + 0.6 4.6. 9.5 + 0.9 3.6 

Terminating All customers 
switch Mean (sec) Std. dev. (sec) Legend: 
PAN Panel 

SXS 12.7 + 0.6 4.6 SXS_ Step-by-step 

CDO 12.8 + 1.0 4.7 5XB_ No.5 crossbar 

5XB 9.7+0.7 5.0 1XB_ No. 1 crossbar 

1XB 9.5 + 0.6 3.8 ESS Electronic Switching System 
ESS 7.9+40.5 3.7 CDO Community Dial Office 


longer delay between the end of dialing and the beginning of audible ring 
or busy than attempts which originate from rotary dial stations served 
by SXS switching machines. TOUCH-TONE dialing signals, which are 
generated at a subscriber station, are stored and converted to dial pulses 
at aSXS office before the telephone number is processed. On the other 
hand, dial pulses, which are generated at a rotary dial station, directly 
drive a SXS switch during the dialing interval itself; therefore, the 
telephone number is almost completely processed at the end of dialing. 
These operational procedures account for much of the difference be- 
tween the SXS estimates; however, unknown differences in calling dis- 
tances, routing patterns, and toll switching may also contribute to the 
3.9 second difference in average setup times. There are no significant 
differences between rotary dial and TOUCH-TONE results for the other 
types of switching. 

The average setup time for attempts originating from customer sta- 
tions served by SXS machines is significantly higher than the averages 
for those originating through No. 5 crossbar (5XB), No. 1 crossbar (1 XB), 
and ESS machines. In the previous paragraph it was stated that the dial 
pulses generated from a rotary dial station directly drive a SXS switch, 
and hence, at the end of dialing the dialed number is almost completely 
processed by the originating SXS machine. This would tend to shorten 
the interval from end of dialing to ring or busy. However, the same SXS 
procedure occurring at the terminating switching office tends to lengthen 
the setup interval as is illustrated by the SXS and CDO estimates in the 
lower portion of Table VII. An estimated 50 percent of all DLSO-observed 
call attempts which originate through SXS machines also terminate 
through SXS machines. This terminating percentage is 20 to 30 per- 
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centage points higher than for the other types of switching machines. 
In addition to this influence, which tends to increase the length of the 
call setup interval, there is a second influence which comes into play: the 
calling distance distribution for DDD attempts originating from cus- 
tomers served by SXS machines is significantly biased towards longer 
calling distances relative to the distributions for attempts originating 
from customers served by the other types of switching machines. In the 
previous subsection it was shown that call setup time increases as the 
calling distance gets longer. Thus, the long setup times for rotary dialed 
calls handled by SXS machines are not caused by the originating SXS 
machines. They come about because of call destination and calling dis- 
tance characteristics associated with the attempts. However, SXS ma- 
chines at the terminating ends of calls do cause significantly longer 
setups. 

In addition to the SXS findings above, the results in Table VII indicate 
that calls placed by customers served by panel machines have rather slow 
setup times and calls placed or received by customers served by ESS 
machines have rather fast setup times. Terminating ESS machines 
generally provide immediate ringing when a desired line is not busy, 
while other machines generally do not. This difference in procedures 
tends to shorten setup times for calls terminating through ESS machines. 
Other causes which underly these characteristics have not been uncov- 
ered. The precaution of not totally attributing these findings to the panel 
or ESS machines is exercised because of the absence of routing and toll 
switching information. 


V. CONCLUSIONS 


This survey has confirmed that both network performance and cus- 
tomer behavior play a role in determining the overall completion prob- 
ability and call setup time experienced by telephone customers. In the 
area of completion performance, the percentage of equipment blockages 
and failures has been steadily declining in recent years due to extensive 
efforts throughout the Bell System to provide adequate facilities to 
achieve low blocking and to detect and correct equipment failures. As 
a result, customer-dependent failures now account for at least 85 percent 
of all failures to complete. Thus the greatest potential for improved 
completion performance lies in reducing such failures. This may be 
brought about, for instance, through the offering of new services made 
possible by the advent of Electronic Switching Systems which will serve 
to reduce the occurrence of called customer did not answer and called 
line busy failures. | 

One manifestation of the dominance of customer behavior on com- 
pletions is the important influence of class of subscriber service. Many 
trends by day of week, time of day, and calling distance can be directly 
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attributed to changes in the relative proportion of business and residence 
traffic. 

The total call attempt time is made up of several components, most 
of which are customer-controlled, such as dialing time and time from 
start of ringing to answer by the called customer. The one network- 
dependent component, time from end of dialing to network response, 
represents about 29 percent of the total call attempt time. Calling dis- 
tance and type of local switching have both been shown to affect this 
interval. The advent of common channel interoffice signaling is expected 
to reduce this component from its current value of 10.9 seconds to 2 or 
3 seconds, which will reduce the overall call setup time by 20-25 per- 
cent. 

Customer holding time after the onset of ringing and after a positive 
indication of failure to complete has a significant effect on total call at- 
tempt time and hence on network load. Programs directed toward en- 
couraging faster customer abandonment or toward releasing network 
facilities more promptly after receipt of busy of no circuit/reorder signals 
would reduce this load on network facilities. 
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A recent Bell System customer loop survey was Jointly conducted 
by AT&T and Bell Laboratories. Bell Labs developed the sampling plan 
and performed the data analysis while AT&T assumed responsibility 
for the sample selection, data collection, and all intermediate data 
processing prior to analysis. This paper presents the principal physical 
composition and calculated transmission statistical characterization 
of customer loop plant as defined by that survey. Comparisons of data 
obtained from this survey and a similar survey made in 1964 are also 
presented. 


Il. INTRODUCTION 


A recent Bell System customer loop survey was jointly conducted by 
AT&T and Bell Laboratories. Bell Labs developed the sampling plan 
and performed the data analysis while AT&T assumed responsibility 
for the sample selection, data collection, and all intermediate data pro- 
cessing prior to analysis. 

This paper presents the principal results of this Bell System customer 
loop survey. This survey provides a statistical characterization of both 
the physical composition and the calculated transmission characteristics 
of customer loop plant. Comparisons of data obtained from the survey 
and a similar survey made in 1964 are also presented. 

The characterization of the Bell System customer loop plant as defined 
by the 1964 loop survey was published in 1969 by P. A. Gresh.! 

The current Bell System survey consisted of a simple random sample 
of 1100 main stations selected from the population of all Bell System 
main stations. Lines serving official telephones, foreign exchanges, dial 
teletypewriter exchanges (TWX), and special services were, however, 
omitted from the survey, as it was felt their design would not be repre- 
sentative of customer loop plant. 
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Fig. 1—Geographic distribution of sampled loops, 1973 customer loop survey. 


Analysis of the 1100 sampled main stations consisted of encoding and 
keypunching the raw data to enable computer processing. A number of 
computer programs were then accessed which in addition to tabulating 
the physical composition of subscriber loops also calculated and evalu- 
ated the transmission performance of the customer loop plant. 


ll. DESIGN OF THE SURVEY 


Definition of the sampling population, sampling plan, and sample size 
were based on the statistical model used in the 1964 loop survey. The 
population for the current survey was defined as all 57,293,521 Bell 
System main stations as of January 1, 1973. A simple random sample 
was drawn with the size determined so as to provide precision equal to 
that obtained in the 1964 survey. The precision was measured by the 
width of the 90% confidence interval for the estimated mean working 
length. The desired confidence interval (at 90 percent confidence level) 
of +450 feet on the average cable distance to the sampled main stations 
dictated a sample size of 1100 randomly selected main stations. The 
actual sample produced a +490 foot confidence interval. The survey had 
a wide geographical dispersion with every state except Wyoming con- 
tributing to the sample and heavier contributions from the metropolitan 
areas (Fig. 1). 


lll, LOOP SURVEY RESULTS—PHYSICAL COMPOSITION 


Data obtained in the 1973 loop survey included a detailed physical 
makeup diagram of each of the sampled loops. Distributions of such 
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Table | — 1973 Loop Survey: Summary of main station 
characteristics 


Significant 
90% Confidence level for 
limits on mean difference 
Main station Mean (+) of means, 
quantity 1964 1973 1964 1973 % 
Working length 10,613 ft 11,413 ft 476 ft 490 ft 99 
Total bridged tap 2,478 ft 1,821 ft 172 ft 113 ft 99+ 
Working bridged 
tap 228 ft 121 ft 74 ft 53 ft 95 
Airline distance 7,758 ft 8,410 ft 386 ft 395 ft 95 
Working length 
airline distance 1.50 1.51 .03 05 <80 


Note: Drop wire excluded except when individual lengths exceed 400 feet. 


items as type of cable construction, composition by gauge, pair size of 
cables and grade of service (individual or multiparty service) distribution 
were generated in addition to such physical quantities as working length 
to main station and total bridged tap. These quantities are defined in 
the illustration on the next page. Similar data was gathered in 1964 and 
a comparison of the two surveys had been made to detect changes in 
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Fig. 2—Working length to main station. 
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plant composition over the intervening 9 years between the two sur- 
veys. 

Table I gives the means and confidence limits for the principal physical 
properties of both the 1964 and 1973 surveys and the level of significance 
for the observed differences. Cumulative distributions of these quantities 
are presented in Fig. 2 through 5. The distribution of working bridged 
tap is not given, since 93 percent* of the sampled main stations were 
served by loops having zero working bridged tap and consequently the 
distribution is not particularly useful. 

As Table I indicates, the estimated average working length from the 
central office to the main station is 11.4 kilofeet with 90 percent proba- 
bility that the true mean value lies within +490 feet of this estimate. 
Notice that this is 800 feet longer than the average in 1964 and that this 
observed increase can statistically be said to be a true increase. 


Loop Sketch Showing Main Station Quantities 


—_— ———— 6600 FT——_____» 


a oe, =. ew 

<——___—» 7600 FT 7800 FT 8200 FT 8350 FT 
: Z Zs = 

-— 7% SAMPLED MAIN ANOTHER 

7850 FT STATION PARTY 
Main Station Quantities 

Working length to sample (WL) 7800 ft 
Total bridged tap (TBT) 800 ft 
Working bridged tap (WBT) 400 ft 
Airline distance to sample (AL) 6600 ft 
Working distance/airline distance (ratio) 1.238 


Analyses of principal loading characteristics of loop plant are given 
in Table II. As indicated, the percent of loaded loops has increased from 
16.4 to 22.9 percent. The majority of this 6.5 percentage point increase 
is attributable to an increase in loops requiring loading (main stations 
served by loops over 18 kilofeet in length). A second contribution to the 
increase is attributed to the fact that greater care has been taken to as- 
sure that loops in excess of 18 kilofeet are loaded (long nonloaded loops 
have decreased from 1.5 percent to 1.0 percent). The final contribution 


* 92 percent of the sampled main stations were individual lines and 14 percent of the 8% 
party-line stations were working alone. 
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Fig. 4—Airline distance to main station. 
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Fig. 5—Ratio of working length to airline distance. 


to the increase in loaded loops is the fact that the percentage of main 
stations served by short loaded loops has increased (from 4.6 percent to 
6.7 percent). This increase in short loaded loops can not be accounted 


Table Il — 1973 Loop survey: Summary of principal loading 
information 


1973 


90% Change 
Confi- Since 1964 





dence Level of 
; 1964 Interval Signifi- 
Main stations served by %MS %MS (+%) cance 
Loaded loops 16.4 22.9 2. 99+ 
With H88 15.5 22.6 2.1 99+ 
With loading other than H8 0.9 0.3 0.2 93 
Less than 18 kf central office to main 
station 
and loaded* 4.6 6.7 1.2 97 
Having load spacing deviations 
exceeding +500 feet 5.3 8.6 1.4 99+ 
exceeding +120 feett NA 18.5 1.9 — 
Nonloaded loops 83.6 77.1 2.1 99+ 
Long nonloaded loops (>18 kilofeet) 1.5 1.0 0.5 <80 


+ The normal confidence intervals for percentages close to 0 are not meaningful and so 
are indicated by (—). 


* Based on distance from CO. to sampled main station and does not necessarily indicate 
improper design. 
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for by an increase in PBX/centrex services (PBX/centrex service has in- 
creased from 1.8 percent in 1964 to 2.5 percent in 1973). 

Survey data analysis also provides information on loop plant com- 
position by gauge, pair size, type of construction (aerial, buried, or un- 
derground), and grade of service as a function of distance to the sampled 
main station. 

The sampled loops were examined at 1 kilofoot intervals, beginning 
at the central office, to obtain the various distributions presented in Figs. 
6-9. The gauge distribution at 10 kilofeet, for example (Fig. 7a), was 
derived by determining the gauge of each loop sample at the 10 kilofoot 
point of all 494 loops which extend beyond 10 kilofeet. A number of 
conclusions can be drawn from these analyses. First, inspection of pair 
sizes as a function of distance from the central office (Fig. 6a) reveals 
decreasing pair size with distance. But also notice that there has been 
a shift toward the use of larger pair size cables for longer distances as 
compared to 1964 (Fig. 6b). For example, at 30 kilofeet only 47 percent 
of the loops are in pair sizes of 100 pairs or less in 1973 while in 1964, at 
the same distance, nearly 75 percent of the loops were served by cables 
of 100 pairs or less. Similarly, examination of the gauge distribution (Fig. 
7a) shows a transition to coarse gauge with increasing distance (at 20 
kilofeet only 30 percent of the sampled loops consist of finer than 22 
gauge cable). It can also be seen (Fig. 7b) that there has been a shift since 
1964 toward use of the finer gauge cables at further distances from the 
central office. For example at 10 kilofeet, 28 percent of the loops sampled 
in 1973 were contained in 26 gauge cables while in 1964 only 15 percent 
of the loops were contained in 26 gauge cables at this distance. 

As shown in Fig. 8a, longer loops are predominately either aerial or 
buried (at 35 kilofeet, 52 percent of the sampled loops were constructed 
with buried facilities and 48 percent with aerial). Figs. 8b—d show that 
in the intervening 9 years between the two loop surveys there has been 
an increase in buried cable coupled with decreasing use of aerial facilities 
with greater distance from the central office. For example, in the 1973 
survey at the distance of 20 kilofeet, 36 percent and 49 percent of all plant 
was buried and aerial respectively, while in 1964 the values were 12 
percent for buried and 72 percent for aerial (aerial facilities consist of 
both wire and cable). Noteworthy also is the grade of service distribution 
(Figs. 9a—c), which shows with greater distance from the central office 
the increasing use of individual party service coupled with more uniform 
distribution of two-, four-, and eight-party service than was prevalent 
in 1964. 


IV. LOOP SURVEY RESULTS—TRANSMISSION PERFORMANCE 


The 1973 loop survey data also provided sufficient information to 
derive the equivalent “T” network for each loop which was then used 
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Fig. 6—(a) Pair size distribution. (b) Pair size distributions, 1973 vs. 1964 surveys. 
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Fig. 7—(a) Gauge distribution. (b) Gauge distributions, 1973 vs 1964 surveys. 
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Fig. 8—Type of construction. (a) Distribution. (b) Buried cable. 


to derive the calculated transmission loss at nine voice-band frequencies. 
In the 1964 loop survey, the loss of each sampled loop was determined 
by two methods: theoretical calculations as described above, and actual 
measurements. Comparison of the calculated and measured transmission 
losses showed that calculated losses based on the physical composition 
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Fig. 8—(c) Aerial (cable and wire). (d) Underground cable. 
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Fig. 9—Type of service. (a) Distribution. (b) Individual lines and two-party. (c) Four- 
party and eight-party. One-, two-, four-, and eight-party include residence service only. 
Business includes PBX, centrex, and coin. 


of the loops were sufficiently close to the measured transmission losses. 
Consequently, actual transmission measurements were not made on the 
sampled loops in the 1973 survey. The distribution of calculated insertion 
losses at 1, 2, and 3 kHz, ascatter diagram of the 1 kHz insertion losses, 
return loss at 3 kHz, echo return loss (equal weighting over the 500 to 
2500 Hz band), and dc resistance are presented in Figs. 10-13. 

The cumulative distributions of calculated 1, 2, and 3 kHz insertion 
losses for customer loop plant is presented in Fig. 10. Note that ap- 
proximately 98 percent of all Bell System main stations are served by 
loops having 1 kHz insertion loss less than or equal to 8 dB with a cal- 
culated mean loss of 3.7 dB. Similarly, at 3 kHz the calculated mean loss 
is 7.5 dB with 95 percent of all main stations served by loops having less 
than 15 dB insertion loss. A scatter diagram of 1 kHz calculated insertion 
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Fig. 11—Scatter diagram of calculated 1-kHz insertion losses. 
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Fig. 13—Calculated resistance to main station. 


loss as a function of loop length is shown in Fig. 11. Table IV presents 
calculated insertion loss at nine voiceband frequencies. Note that mean 
insertion losses across the voice band have increased, principally because 
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Table IV — Insertion loss (calculated with 900-ohm terminations 
at both ends) 


1964 1973 Difference 
90% Confidence 90% Confidence St. 1964-1973 
Freq., Mean, Interval, Mean, Interval, Dev., Significance 
Hz dB dB (+) dB dB (+) dB Level 
200 2.4 0.06 2.9 0.16 3.2 99+ 
300 2.5 0.07 2.8 0.10 2.1 99+ 
500 2.7 0.08 3.0 0.11 2.3 99+ 
1000 3.5 0.10 3.7 0.14 2.8 94 
1500 4.3 0.13 4.6 0.16 3.2 98 
2000 5.3 0.16 5.4 0.18 3.6 <80 
2500 6.2 0.18 6.4 0.21 4,1 <80 
3000 7.3 0.21 7.5 0.23 4.7 <80 
3200 7.9 0.23 8.3 0.26 5.2 93 


of an increase in the use of finer gauge cable and the resulting increase 
in average loop resistance (Fig. 13). 

Calculated return loss performance is presented in Table V for nine 
frequencies and the cumulative distributions for the 3 kHz return loss 
and echo return loss is shown in Fig. 12. Note that it has remained es- 
sentially unchanged since 1964. 
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Fig. 14—Nonloaded loop input impedances at 1 kHz. RL circles based on 500-type 
telephone set impedance (22-gauge H88 cable termination at the central office). 
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Table V — Return loss (calculated using a 900 ohm + 2uF 
termination) 


1964 ee en Difference 
90% Confidence 90% Confidence St. 1964-1973 
Freq., Mean, Interval, Mean, Interval, Dev., Significance 
Hz dB dB (+) dB (dB (+) dB Level 
200 8.0 0.11 8.3 0.10 2.0 99+ 
300 10.2 0.12 10.3 0.11 2.2 <80 
500 13.4 0.17 13.1 0.16 3.2 97 
1000 15.4 0.30 15.3 0.31 6.2 <80 
1500 13.1 0.27 13.3 0.28 5.6 <80 
2000 10.9 0.25 11.2 0.27 5.3 82 
2500 9.1 0.20 9.1 0.20 4.0 <80 
3000 7.7 0.16 Teal 0.16 3.2 <60 
3200 7.1 0.15 7.0 0.15 2.9 <60 
Echo 11.2 0.15 11.3 0.15 3.1 <80 


The final transmission characteristic to be presented is the loop input 
impedance at the station set. The input impedance of each loop was 
computed by using the equivalent “T”’ network for each loop. 

Central office terminations equivalent to four-wire trunks and two- 
wire trunks were represented by 900 ohms plus 2.16 uF and midsection 
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Fig. 15—Loaded loop input impedances at 1 kHz. RL circles based on 500-type telephone 
set impedance (22-gauge H88 cable termination at the central office). 
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Fig. 16—Nonloaded loop impedances at 3 kHz. RL circles based on 500-type telephone 
set (22-gauge H88 cable termination at the central office). 


input impedance of 22-gauge H88 loaded cable, respectively. Intraoffice 
calls were simulated based on a Monte Carlo technique of selecting 500 
pairs of loops. Using a random number generator a loop was chosen from 
the 1100 loops in the 1973 loop survey data base. The central office ter- 
mination was the input impedance as measured at the central office of 
another randomly chosen loop. 

Scatter diagrams of loop input impedance at the station set are pre- 
sented for 1 and 3 kHz voiceband frequencies. There is wide variation 
in the input impedance of loaded and nonloaded loops; consequently 
the two populations (849 nonloaded and 251 loaded loops) were sepa- 
rated for this study. Return loss circles at 0.1, 1, 2, 3,5, and 7 dB are also 
superimposed on the scatter diagrams. The return loss circles were based 
on the 500-type telephone set as the reference impedance since 84.3 
percent of the sampled loops in the 1973 loop survey were terminated 
with this type of telephone set. Further, since loaded loops have nearly 
15 mA lower loop current than nonloaded loops, the mean loop currents 
of the individual populations were used. 

The scatter diagrams for loops with simulated two-wire trunk termi- 
nations (22-gauge H88 loaded cable) are presented in Figs. 14 through 
17. The dispersion of points can be partially attributed to the use of 
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Fig. 17—Loaded loop input impedances at 3 kHz. RL circles based on 500-type telephone 
set impedance (22-gauge H88 cable termination at the central office). 
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Fig. 18—Median input impedances at all nine frequencies for loaded and nonloaded 
loops (22-gauge H88 cable termination at the central office). 
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Fig. 19—Nonloaded loop input impedances at 1 kHz. RL circles based on 500-type 
telephone set impedances (900 ohm + 2.16 uF cable termination at the central office). 
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Fig. 20—Loaded loop input impedances at 1 kHz. RL circles based on 500-type telephone 
set impedances (900 ohm + 2.16 uF cable termination at the central office). 
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Fig. 21—Nonloaded loop input impedances at 3 kHz. RL circles based on 500-type 
telephone set impedances (900 ohm + 2.16 uF cable termination at the central office). 


REACTANCE, OHMS 


400 g 


200 


- 


-200 


-400 


-600 


~800 


~1000 


-1200 


~1400 
50 250 450 650 850 1050 1250 1450 1650 1850 


RESISTANCE, OHMS 


Fig. 22—Loaded loop input impedances at 3 kHz. RL circles based on 500-type telephone 
set impedances (900 ohm + 2.16 uF cable termination at the central office). 
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Fig. 23—Median input impedances at all frequencies for loaded and nonloaded loops 
(900 ohm + 2.16 uF cable termination at the central office). 
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Fig. 24—Intraoffice simulation input impedances at 1 kHz. RL circles based on 500-type 
telephone set impedance. 
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Fig. 25—Intraoffice simulation input impedances at 3 kHz. RL circles based on 500-type 
telephone set impedance. 


coarser gauge cable than is required to meet the resistance design ob- 
jective of 1300 ohms. Notice that loaded loops (Figs. 15 and 17) are more 
scattered due to their vulnerability to load coil spacing irregularities than 
are the nonloaded loops (Figs. 14 and 16). Median impedances at all nine 
voiceband frequencies for simulated two-wire trunk terminations are 
shown in Fig. 18. The scatter diagrams for loops with simulated four-wire 
trunk terminations (900 ohms plus 2.16 uF) are shown in Figs. 19 through 
22. Similarly, the scatter plots indicate more irregularities present in 
loaded loops than nonloaded ones. Curves of the median input imped- 
ances at all nine voiceband frequencies of loaded and nonloaded loops 
with a four-wire trunk termination at the central office are presented 
in Figure 23. The intraoffice call simulation results can be seen in the 
scatter diagrams presented in Figs. 24 and 25. These reveal the effect 
of connecting together two random loops, one terminated by a station 
set and the other’s impedance calculated at the station set. The median 
input impedances for intraoffice calls are presented in Fig. 26 along with 
median impedances of simulated two-wire and four-wire trunk termi- 
nations of nonloaded loops. 
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Fig. 26—Median input impedances at all nine frequencies for simulated intraoffice calls, 


simulated two-wire trunk terminations (22-gauge H88 loaded cable), and simulated 
four-wire trunk terminations (900 ohm + 2.16 uF cable) of nonloaded loops. 


V. CONCLUSIONS 


The analyses of the 1973 loop survey can be summarized in the fol- 
lowing four principal results. 

(i) The average customer loop length is currently 11.4 kilofeet with 
only 4 percent of the main stations located beyond 30 kilofeet from their 
serving office. The length distributions show a statistically significant 
(99 percent level of significance) increase in average loop length since 
1964, with the average loop length increasing by 800 feet. 

(ii) The average calculated 1 kHz insertion loss of the Bell System 
loop plant is currently 3.7 dB* with 95 percent of all main stations being 
served by loops having a calculated 1 kHz loss of less than 7.5 dB. At 3 
kHz the average calculated loss is 7.5 dB with 95 percent of all main 
stations served by loops having less than 15 dB insertion loss. 


* Anomalies in plant design and record errors will tend to make actual losses slightly 
higher than calculated losses. For example, the 1964 loop survey, which included actual 
loss measurements as well as calculated losses, showed that average measured loss at 1 
kHz was approximately 0.3 dB higher than the calculated value. 
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(111) The percentage of main stations served by loaded loops has in- 
creased since 1964 from 16.4 to 22.9 percent. Sixty percent of the increase 
is a direct result of the increase in the percent of customer loops requiring 
loading (loops longer than 18 kilofeet). The remainder is attributed to 
an increase in the use of loading for loops shorter than 18 kilofeet and 
a reduction in the percentage of nonloaded loops longer than 18 kilo- 
feet. 

(iv) Average dc resistance to the sampled main station has increased 
from 574 to 646 ohms since 1964. This.increase is attributed to the trend 
towards the use of finer gauge cable in the loop plant combined with an 
increase in the average customer loop length. 
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Level Reassignment: A Technique for Bit-Rate 
Reduction 


By B. PRASADA, F. W. MOUNTS, and A. N. NETRAVALI 
(Manuscript received June 17, 1977) 


Three techniques are described for reducing the entropy of a digitally 
coded monochrome picture signal by dynamically changing the 
input-output relationships of a single quantizer. This is done by re- 
assigning the input levels of the quantizer to different representative 
output levels in such a way as to reduce the entropy of the quantized 
output, while keeping the visibility of the resultant quantization noise 
below a certain specified threshold. These techniques are simulated 
both in the pel and the transform domain for a simple differential pulse 
code modulation system and Hadamard transform coding system, re- 
spectively. It is found that they reduce the entropy by about 20 to 25 
percent in the pel domain without noticeably sacrificing the picture 
quality. In the Hadamard transform domain, the entropy of the first 
transform coefficient can be decreased by about 23 percent with little 
change in picture quality. 


I. INTRODUCTION 


It is well known that the statistics of picture signals are nonstationary 
and that the required fidelity of reproduction demanded by the human 
eye varies from picture element (pel) to picture element. Consequently, 
for efficient digital representation of pictures, it is desirable to adapt 
coding strategies to those local properties of the picture signal which 
determine the visual sensitivity to quantization noise. In this study, we 
make use of the spatial masking properties of the human observer to 
adapt the quantization strategies for encoding the picture signal. We 
define spatial masking as the reduction in the ability of a person to vi- 
sually discriminate amplitude errors which occur at or in the neighbor- 
hood of significant spatial changes in the luminance. To this end, we 
borrow, from our earlier work, measures of luminance spatial activity 
in both the transform! and the pel? domain, and the relationships 
(called the visibility functions) between the amplitude accuracy and 
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these measures of spatial detail. We use these visibility functions to 
change dynamically the input-output mapping of a single quantizer to 
reduce the bit rates. This is done by reassigning the input of the quantizer 
to a different representative level than normal in such a way as to reduce 
the entropy of the quantized output, while keeping the visibility of 
quantization noise below a certain threshold. We demonstrate three 
different algorithms for level reassignment and evaluate their potential 
by measuring the entropy of the quantized output for a given picture 
quality. It is worth pointing out that there are other methods of adapting 
the quantizer characteristics. Some of these are discussed in Refs. 3, 5, 
and 6. 


1.1 Summary of the approach 


To illustrate the coding strategy, let ein and @out be the input and 
output to a quantizer, whose decision levels are denoted by X;,i =1,..., 


N + 1 and representative levels by Y;,i = 1,..., N. The quantizer 
input-output mapping in the absence of adaptation is: 
If X; < ein S Xi41, theneour= Yi, 1=1,...N (1) 


In the adaptive quantization, instead of representing the input ejn as 
Y;, we change the representation to Y; (J + i) provided the visibility of 
the quantization noise (ej, — Y;) is below a given threshold. The changed 
representation Y; is chosen, in general, to decrease the entropy* of the 
quantizer output. This is done in several different ways. In one case, we 
increase the frequency of occurrence of the inner levels? of the DPCM 
quantizer. Here we make use of the fact that the frequency of occurrence 
of the inner levels of a DPCM quantizer is generally high, and, therefore, 
making it higher decreases the entropy. In another case, for alternate 
levels Y;, Y; is taken to be Y;-1 or Y;+1, depending on the noise visibility. 
This has the effect of decreasing the frequency of occurrence of alternate 
quantizer levels, and increasing the frequency of occurrence of the other 
levels, and thus reducing the entropy. We refer to such a change in the 
use of the quantizer levels as level reassignment. Spatial detail sur- 
rounding the pel and the corresponding visibility function are used in 
determining whether quantization noise (ej, — Y;) is visible. We have 
studied three different algorithms which differ only in detail from the 
approach outlined above. They are described in Section III. 


* We define the entropy of the quantizer output as 


ni 








N nj 
ea logs 
i=l N 
J 


2, nj L nj 
=1 jJ=1 
where n; is the frequency of occurrence of level Yj. 

t Inner levels are those which are close to level zero. 
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1.2 Summary of results 


We simulated the level reassignment algorithms on a digital computer. 
The efficiency of an algorithm was measured by evaluating the difference 
in entropy between a nonadaptively coded picture and an adaptively 
coded (using level reassignment) picture having the same quality. In the 
pel domain, reassigning all the quantizer levels to the lowest level per- 
mitted by the quantization noise visibility leads to degradation at edges 
in the form of slope overload and, therefore, allows only about 9 percent 
reduction in entropy before picture quality is degraded. To avoid deg- 
radation of these edges, only a few inner levels (three inner levels on both 
sides of zero of a 15-level quantizer) were reassigned. This allows a higher 
threshold for the quantization noise before it is visible and leads to about 
14 percent decrease in entropy. Alternate reassignment, in which certain 
alternate levels are reassigned to one of their two surrounding levels 
depending on the visibility of quantization noise, also allows a higher 
threshold for quantization noise and results in a 25 percent decrease of 
entropy. Delayed reassignment, in which reassignment is done only if 
the reassigning of the present sample to a lower level does not “adversely” 
affect the coding of the next sample, gives a 22 percent decrease in en- 
tropy. Thus alternate reassignment and delayed reassignment appear 
promising. 

In the transform domain, we considered differential PCM coding of 
the first Hadamard coefficient. Using alternate reassignment, it was 
possible to reduce the entropy of the coded first transform coefficient 
by about 23 percent. We also give results on different methods of eval- 
uating the visibility of the quantization noise and the effect of changing 
the picture content. 


ll. MASKING AND VISIBILITY FUNCTIONS? 


We have constructed a simple measure of luminance spatial activity. 
We call this measure the masking function. In the pel domain, the 
masking function at a pel is taken as the weighted sum of the luminance 
slopes at the pel under consideration and at the neighboring pels. Thus, 
at point (i, 7) the two-dimensional masking function M;,j, using a 3 X 3 
neighborhood of slopes is given by 


itil j+1 ie H V 
ge aa eae a (2) 
n=1—-1 t=)— 


where ||(n, ¢) — (i, /)|| is the Euclidean distance between points (n, t) and 
(i, j) normalized by the distance between horizontally adjacent pels; m/,, 
m x , are the horizontal and vertical slopes of the image intensity at point 
(n, t); wis a constant taken to be 0.35 based on the tests given in Ref. 3. 
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Fig. 1—Definition of Hadamard coefficients. A, B, C, D are the pel positions, and H, 
Ho, Hz, H4 are the Hadamard coefficients. 


In this framework, many masking functions are possible and some of 
these are discussed in Ref. 3. 

In the Hadamard transform domain!” using a 2 X 2 transform of a 
block of pels in the same field defined as in Fig. 1, the measure of spatial 
activity in a block is taken as 


H = max(|Hol, |H4)]) (3) 


It is clear from Fig. 1 that H is the maximum magnitude of the average 
line or element difference of pels within the block. 

We hypothesized that the masking of the quantization noise is related 
to the spatial detail as measured by the masking functions (M or H). The 
precise relationship was obtained through subjective tests.* In these 
tests, the test picture was obtained by adding varying amounts of noise 
(to simulate the quantization noise) to all pels (or blocks in the case of 
transform domain) where the measure of spatial detail (M or H) has a 
given value. Such a picture was then compared in an A/B test with an 
unimpaired picture to which the subjects added controlled amounts of 
white noise everywhere in the picture until they found both pictures to 
be subjectively equivalent. 

Visibility functions were derived from such subjective equivalence. 
They are shown in Fig. 2. Figure 2a shows the visibility of noise added 
to a pel as a function of masking function M, at that pel, whereas Fig. 
2b shows the visibility of noise added to Hadamard transform coefficient 
H, asa function of max(|H9], |H4|). Both these visibility functions de- 
crease with respect to their argument, implying that at higher values of 
spatial detail the visibility of quantization noise is lower. Details of the 
subjective tests and visibility functions are given in Refs. 1-3. From this 
work we know that the visibility of noise at a point is proportional to the 
power of noise and the proportionality constant is given by the value of 


* These tests are similar to the ones described by Candy and Bosworth.4 
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Fig. 2—(a) Visibility function for Checkered Girl and Teri using a two-dimensional one 
neighbor masking function. (b) Visibility function for noise added to H, as a function of 
max (|Ho|, |H4]). 


the visibility function evaluated using the spatial detail at that pel. ‘Thus 
visibility functions allow us to judge the visibility of the quantization 
noise. 


Hil. DETAILS OF THE TECHNIQUES 


The entropy-reducing techniques described briefly in Section I have 
been simulated on a computer. In this section we give details of the 
technique as well as the results of the simulation. We have taken two 
different head and shoulders views (Fig. 3a, called Checkered Girl, Fig. 
3b, called Teri) and processed them using our algorithms. We are dealing 
with a still picture, but since the quantization noise in a real television 
system varies from frame to frame even with still pictures, we coded three 
frames of the same picture to which a small amount of random noise (to 
simulate camera noise) was added. Due to the randomness of this noise 
from frame to frame, the noise appeared to move when the three frames 
were displayed in a repetitive sequence. 
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(a) (b) 


Fig. 3—Original pictures of (a) Checkered Girl and (b) Teri. These are 256 X 256 arrays 
with signal having a bandwidth of 1 MHz; each sample is linearly quantized with 8-bit 
PCM. 


In the pel domain, we started with a DPCM coding system which had 
a 15-level quantizer optimized for minimum mean square quantization 
error. Using this quantizer, the coded pictures had noticeable defects, 
but were considered to be of acceptable quality by the authors. The en- 
tropy of the quantized output for Checkered Girl is 2.93 bits/pel and for 
Teri 2.64 bits/pel. In our experiments, the parameters of the coder were 
varied to decrease the entropy while keeping the quality of the coded 
picture similar to that of the 15-level DPCM coded picture. Thus the ef- 
ficiency of the adaptive techniques described in the next section was 
judged by the difference between the entropy of the output of the 
adaptive coder and the 2.93 bits/pel (for Checkered Girl) and 2.64 bits/ 
pel (for Teri) required for the nonadaptive coder. In the transform do- 
main, we applied alternate reassignment to the DPCM coding of the first 
Hadamard coefficient using the first Hadamard coefficient of the hori- 
zontally adjacent block as the predictor. We used a 37-level quantizer 
which resulted in a picture of acceptable quality as judged by the authors 
and entropy of 3.75 bits/block for Checkered Girl and 3.76 bits/block 
for Teri for coding of H. 


3.1 Reassignment to the lowest level 


The 15-level DPCM quantizer characteristic that was used in our 
simulations is shown in Fig. 4. The level numbers are also marked in the 


0 3 8 15 24 33 42 58 
ee os > ao i > a a oo 
QO 1 25 6 11 12 19 20 28 29 37 38 49 50 255 
LEVEL 
NO. 0 1 2 3 4 5 6 7 


Fig. 4—A 15-level symmetric quantizer used for DPCM processing in the pel domain. 


The symbol 


z 
a 


x y 
means that inputs between x and y including x and y are represented by z. 
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Table | — Change in quantizer level distribution from nonadaptive 
to adaptive coder which reassigns all levels to the lowest possible 
level for Checkered Girl. Only non-negative levels are shown. 


ith level of adaptive coder 
2 3 4 5 





I 6 7 
Total 

pel count 
perlevel 26566 10624 3029 1874 1206 678 363 590 
S 0 16825 16825 0 0 0 0 0 0 0 
%& 8 1 14465 5465 — 9000 0 0 0 0 0 0 
=o 2 6161 3424 1286 1451 0 0 0 0 0 
$2 3 3307 532 3802.—Ss«i1474— 999 0 0 0 0 
<2 4 1804 180 22 84 817 701 0 0 0 
su 5 965 78 10 13 43 480 341 0 0 
™2 6 663 27 3 5 8612 21 325 270 0 
= 7 740 35 1 a 3 4 12 93 590 





figure. Let the decision levels for the ith level be X;, Xj+ 1 and the rep- 
resentative levels be Y;. The algorithm for positive 1 proceeds as fol- 
lows: 

If ejn is the input to the quantizer and X; < ej, S Xj41, i.e., without 
adaptation ej, would be quantized as Y;, then, 


(1) Evaluate the quantization error in representing ej, by the next 
lower level, i.e., error = @j, — Y;-1. 
(ti) Compute the visibility of this error by 


evis = lein — i-1|"” -f(M) (4) 


where M is the value of the masking function at the pel, f(-) is the visi- 
bility function, and y is a constant. 

(iit) If the visible error is less than a preassigned threshold then the 
input ej, can be represented as level Y;-1. The above three steps are then 
repeated until the lowest level, at which either the visible error is above 
the threshold or the zero level is reached. 


Similar steps can be followed when ! is negative. Thus the lowest level 
(in magnitude) to which the input can be reassigned such that the visi- 
bility of the quantization is below a given threshold is taken as its rep- 
resentation. In the simulations y was set to 1, and the threshold was 
varied until the lowest entropy was obtained for our picture quality. 
The results of this simulation, for Checkered Girl, are shown in Table 
I. In this table the (, 1)th entry is the number of pels which would be 
quantized by +)th level in a nonadaptive coder but are quantized by +ith 
level in the adaptive coder. The first column gives the total number of 
pels quantized by a level in the nonadaptive coder and the first row gives 
the total number of pels quantized by a level in the adaptive coder. Thus, 
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Table Il — Change in quantizer level distribution from nonadaptive 
to adaptive coder which reassigns the inner three levels to the 
lowest possible level for Checkered Girl. Only nonnegative levels 
are shown. 


ith level of adaptive coder 
0 I 2 3 4 5 6 7 


Total 

pel count 
perlevel 29526 10190 2453 460 1663 827 611 685 
5 0 16950 16950 0 0 0 0 0 0 0 
ws 1 14533 5697 8836 0 0 0 0 0 0 
oo 7229 5483 605 1141 0 0 0 0 0 
2 > 3 3917 1396 749 1312 460 0 0 0 0 
2a 4 1663 0 0 0 0 1663 0 0 0 
So 5 827 0 0 0 0 0 827 0 0 
Re ug 611 0 0 0 0 0 0 611 0 
oS 7 685 0 0 0 0 0 0 0 685 


levels +2 had 6161 pels in the nonadaptive coder out of which 3424 went 
to level 0, 1286 went to level +1, and 1451 remained in level +1, in the 
adaptive coder. Clearly the distribution of pels using different quantizer 
levels in the adaptive coder is more peaked than the corresponding 
distribution in the nonadaptive coder. This causes a decrease in the 
entropy of the coder output. The decrease in the entropy without any 
noticeable change of picture quality is about 0.27 bits/pel which is ap- 
proximately 9.2 percent. 

One of the reasons for such a small decrease of entropy is the feedback 
process inherent in a DPCM coder. If, in quantization of a pel, the 
adaptive quantizer changes its occupancy from level +j to tk (k <j), 
then the differential signal to be quantized for the next pel generally 
increases in magnitude, and there is a good chance that it will occupy 
a higher quantization level than in the case of the nonadaptive coder. 
This, to some extent, hampers our original objective of making the dis- 
tribution of pels in quantizer levels highly peaked. This feedback affects 
the very high levels more severely. Quantizing a “higher level sample”’ 
by a lower level, results in spreading of the picture edges, which is sub- 
jectively very annoying. 

To obviate this, and knowing that the higher quantization levels have 
little effect on the entropy, we used the above algorithm to reassign only 
the inner +1, +2, and +3 levels. The results of this simulation are shown 
in Table II. The Oth level contains more pels than before. This reduces 
the entropy to 2.51 bits/pel which is a reduction of about 14 percent 
without any noticeable sacrifice of picture quality. 

Using the inner +1, +2, +3 level reassignment algorithm, we studied 
the effects of the variation of y in eq. (4). Small y results in good repro- 
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duction of flat areas, but results in severe edge effects similar to the edge 
busyness. A value of y larger than 2, on the other hand, results in better 
edge reproduction but shows granular noise in flat areas. We used three 
different values of y, namely 0.5, 1.0, and 2.0, and varied the threshold 
T such that for each case we got the lowest entropy without changing 
the picture quality. The entropy for y = 0.5 was 2.62 bits/pel, for y = 1.0, 
it was 2.51 bits/pel and for y = 2.0, it was 2.44 bits/pel. Thus, y = 2.0 was 
clearly superior to either 0.5 or 1.0. We used y = 2.0 for all our subsequent 
simulations. 


3.1.1 Alternate reassignment 


Alternate reassignment changes alternate levels (positive or nega- 
tive) to one of the two adjacent levels depending on the quantization 
error visibility. This has the effect of reducing the level occupancy of 
alternate levels, thereby reducing the entropy. Also the highest level 
(positive or negative) was not reassigned for better edge reproduction. 
In our simulation of this scheme only alternate levels i.e., +1, +3, -&5 were 
reassigned. Level 1, for example, was reassigned to either level 0 or level 
2 depending upon both of the following conditions: 


(1) Reassign to level 0 if 


[Cin = Yol|” = lein _ Y.|" (5) 
otherwise reassign to level 2; and 
(it) 
lein — Y,|"f(M) < T (6) 


k = 0 or 2, whichever satisfies condition 1. 


Level 1 was not reassigned if both the above conditions were not met. 
Here again, both the positive and negative levels were reassigned simi- 
larly. This resulted in further improvement in the entropy of coded 
output. We tried alternate reassignment on levels +1 and +3, and +5 
with results as shown in Table III. Note that as expected, the number 
of pels in levels +1, +3, and +5 have decreased significantly. The entropy 
of the quantized output for this simulation was 2.19 bits/pel, which is 
a decrease of about 25 percent. 


3.1.2 Delayed reassignment 


Another way to counter the feedback in the DPCM coder is to wait 
until the next sample for the coding of the present sample and reassign 
the present sample only if the next sample is not “adversely” affected 
by reassigning the present sample. The reassignment of the kth sample 
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Table I!l — Change in quantizer level distribution from 
nonadaptive to adaptive coder which does alternate reassignment 
of levels +1, +3, £5, for Checkered Girl. Only nonnegative levels 

are shown. 


ith level of adaptive coder 
0 1 2 3 4 5 6 7 


Total 

pel count 
perlevel 34802 1845 7868 90 2642 1 852 ~~ 621 
8 0 16274 16274 0 0 0 0 oO 0 0 
es 1 22192 18528 1845 1819 0 0 Oo 0 0 
et 4322 0 0 4322 0 0 0 0 0 
2.2 3 2561 0 0 1727 90 144 0 0 0 
2a 4 1408 0 0 0 0 1408 O 0 0 
SU 5 799 0 0 0 0 490 1 308 0 
“2 6 544 0 0 0 0 0 OO 544 0 
2 621 0 0 0 0 0 oO 0 621 


to the lowest possible level is made only 


(i) If the visibility of error, i.e., |@¢in — Y;|”%-f(M), is less than a 
threshold 7); and 

(ii) If the prediction error for the next sample does not change by 
more than “7”, a specified value. 


In our simulation, we set the threshold 7 at a high value and varied 
T2 between 1 and 15. J. = 5, which implies that reassigning the present 
sample should not change the prediction error of the next sample by 
more than 5, gave the best results. The results of this simulation for 


Table IV — Change in quantizer level distribution from 
nonadaptive to adaptive coder which does delayed reassignment 
of all the levels for Checkered Girl. Only nonnegative levels are 





shown. 
ith level of adaptive coder 
1 2 3 4 5 6 7 
Total 
pel count 


per level 37659 5235 3181 1742 1209 710 504 637 


3S 0 13946 13946 0 0 0 0 0 0 0 
w= S 1 23171 21219 1952 0 0 0 0 0 0 
Reo 7826 2020 3170 2636 0 0 Q 0 0 
oe a 2454 347 97 446 15640 0 0 0 
2a. 4 1405 90 16 90 129 1080 0 0 0 
So 5 804 19 0 8 45 97 635 0 0 
—2 6 591 9 0 1 4 99 70 478 0 

= 7 680 9 0 0 0 3 5 26 637 


~J 
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Fig. 5—A 37-level symmetric quantizer characteristic used for DPCM processing of the 
first Hadamard coefficient. (Only the positive portion is shown.) 


Checkered Girl are shown in Table IV. It is seen that the highest level 
occupancy did not change significantly in the adaptive coder, but oc- 
cupancy in the lower levels changed rather significantly. The resulting 
entropy was 2.28 bits/pel, a decrease of about 22 percent over the no- 
nadaptive coder. It appears that delayed reassignment and alternate 
reassignment perform somewhat similarly in their bit-rate reduction 
capabilities. 


3.2 Reassignment in the Hadamard domain 


Alternate reassignment in the Hadamard domain was simulated for 
the first Hadamard coefficient “H,” of a 2 X 2 block. As mentioned in 
our earlier work,! we did DPCM coding of the first coefficient to make 
use of some of the statistical redundancy not utilized due to small block 
size. It was found that, when all the other coefficients, i.e., H2, H3, and 
H 4, were not quantized, the 37-level quantizer whose characteristics are 
shown in Fig. 5, gave a satisfactory picture without any adaptation. Using 
this quantizer, alternate reassignment was simulated for levels +1, +3, 
+5, +7, +9, +11. The entropy of the quantized output was 2.90 bits/block 
for coding of H;, which amounts to a 23 percent decrease over that ob- 
tainable from the nonadaptive coder. 


3.3 Sensitivity to picture variation 


In order to study the sensitivity of the reassignment algorithms to 
picture variation, we coded Teri using all the above techniques. The 
visibility functions used were those obtained for Teri. The results are 
shown in Table V. It is clear from this table that the percentage decrease 
in entropy using reassignment is a little higher for Teri than for Check- 
ered Girl. For example, using alternate reassignment there is about 30 
percent decrease in entropy. We also coded Teri and Checkered Girl with 
each other’s visibility functions and found relatively insignificant 
changes in entropy. However, the thresholds for the least entropy had 
to be different for the two pictures. Taking the lower value of the 
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Table V — Entropy results with level reassignment for Checkered 


Girl and Teri 

Entropy (bits/pel) 

Checkere 
Algorithm Girl Teri 
Reassigning all levels to lowest level 2.66 2.31 
Reassigning inner +3 levels to lowest level 2.44 2.23 
Alternate reassignment on +1, +3, +5 levels 2.19 1.87 
Delayed reassignment 2.28 1.91 
Hadamard transform coding, alternate* reassignment on 2.90 2.84 


+1, £3, +5, +7, £9, +11 levels 


* Without level reassignment, the entropy of the coded H coefficient was 3.77 bits/block 
for Checkered girl and 3.76 bits/block for Teri. 


threshold and using it for both pictures increased the entropy for one 
of the pictures by about 5 percent. 


IV. SUMMARY 


We have described several coding algorithms to change dynamically 
the input-output relationships of a single quantizer in order to achieve 
bit-rate reduction. The algorithms reassign the input to the quantizer 
to a different representative level in such a way as to reduce the entropy 
of the output, while keeping the visibility of the quantization below a 
specified threshold. The techniques were demonstrated for DPCM coding 
systems both in the pel and the Hadamard transform domain. Although 
we do not discuss it here, it is possible to extend these techniques to PCM 
quantizers. The inherent feedback, which exists in a DPCM system, 
counters the reassignment strategy to some extent. Methods were de- 
vised to decrease the effect of feedback on the bit-rate reduction. It is 
possible to reduce the entropy by about 20 to 30 percent by our adaptive 
techniques, without changing the picture quality. 
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Fiber-Optic Array Splicing with Etched 
Silicon Chips 


By C. M. MILLER 
(Manuscript received August 8, 1976) 


The two-dimensional array splicing concept using aluminum chips 
was reported earlier as a technique that is potentially suitable for 
field-splicing an optical cable containing linear arrays of optical fibers. 
Two arrays are stacked, epoxied, polished, and positioned to forma 
butt-joint splice. This paper reports the use of preferentially etched 
silicon chips to fabricate higher precision arrays which deviate from 
perfect uniformity by only 2.5 um (0.0001 inch) on the average. Average 
losses for splices assembled with these arrays have been in the range 
0.16 dB to 0.32 dB with a yield of 98.8 percent. These average losses are 
close to the anticipated values based on the measured mean offset be- 
tween corresponding fiber axes in the spliced arrays. After several 
combinations of arrays were assembled and loss measurements made, 
the original array configuration was measured to test for array dete- 
rioration, loss measurement repeatability and final alignment re- 
peatability. Only one fiber position showed evidence of contamination, 
and the mean splice loss was repeatable to within 0.02 dB. 


Il. INTRODUCTION 


A two-dimensional array splicing approach for connecting fiber optic 
cable has been reported previously.! This mass splicing technique has 
given encouraging results and was used in the Atlanta Fiberguide System 
Experiment. The method of fabrication of the connector halves involves 
the following operations:!? 


(t) Ribbons containing linear arrays of fibers are prepared by re- 
moving all ribbon and coating material from the fibers at the end of the 
ribbons. 

(it) A stack is formed by interleaving chips and layers of fibers to 
form a two-dimensional array (see Fig. 1). 

(iii) The array is potted to maintain the geometry. 
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Fig. 1—Stacking fiber ribbons. 


(tv) “Good ends” on all fibers in the array are obtained by grinding 
and polishing the array. 

(v) An array splice is made by butting two arrays, prepared as pre- 
viously described; the arrays are held in position using the unoccupied 
grooves on top and bottom of the array for final alignment. 


The purpose of this paper is to give a progress report on the array 
splicing method for connecting groups of optical fibers. Improvements 
in the chips and in the fabrication of the arrays are described. A method 
is given for determining a “figure of merit” for a two-dimensional array. 
Four array-splice loss tests are described and correlation of loss and offset 
is presented. . 


Il. IMPROVEMENTS 


Several improvements have been incorporated in the process for 
fabrication of two-dimensional arrays. These include improved tools and 
fixtures for assembly and end preparation. A stack of 12 prepared rib- 
bons interleaved between 13 chips can be assembled without microscopes 
or micromanipulators in approximately 10 minutes. An improved, 
simplified polishing fixture has replaced the previous version.? This 
fixture uses a negative chip mounted to the inside wall of the holding 
fixture to accurately position the array. 
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0.005 in. 


Fig. 2—Silicon chip cross section. 


The most significant improvement in the array splicing technique has 
been the introduction of preferentially etched silicon chips.? These chips 
were developed and produced by C. M. Schroeder® at Western Electric 
Engineering Research Center using photolithographic techniques ap- 
plied to (100) oriented silicon wafers. Oxide masks are generated on both 
surfaces of the silicon wafer and the wafer is submerged in a basic solu- 
tion such as potassium hydroxide. The etch rate in the unmasked region 
in the (100) direction is much greater than in the (111) direction. Groove 
angles are dependent only on the accuracy of slicing the wafer relative 
to the crystal axis, which is on the order of 1°. Groove opening is de- 
pendent only on the mask accuracy which is approximately +1 micron. 
Alignment between top and bottom grooves is accomplished by using 
a dual exposure channel with near collimated ultraviolet light. Finished 
wafers are laser scribed and broken into individual chips. Figure 2 con- 
sists of two photographs of the cross section of a typical silicon chip. A 
companion paper® describes the chip fabrication process in greater detail. 
Occasionally groove bridges occur (perhaps in 1 percent of the chips) 
which must be detected by visual inspection and discarded. Figure 3 
shows a bridged groove on a silicon chip. 
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Fig. 3—Bridged groove on a silicon chip. 
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Fig. 4—Array splicing. 


Connectors assembled with silicon chips have given greatly improved 
array alignment. Figure 4 shows three views of a connector end at dif- 
ferent magnifications. The completed splice shown in Fig. 4 is one of the 
splices used in the Atlanta Fiberguide System Experiment. 


lll. ARRAY ALIGNMENT MEASUREMENTS 


An accurate method has been developed for determining an alignment 
“figure of merit” which characterizes the degree of uniformity of a 
two-dimensional array. A Leitz-Leica Microscope with a micrometer 
stage is used to measure the coordinates (x,y) of each fiber in the array 
relative to a corner fiber which is used as the reference origin (0,0). The 
accuracy and resolution of these measurement is 1.25 micron. A small 
dark spot which is located at the center of the fiber core* (as seen in Fig. 
4) is used to position a set of cross-hairs accurately and quickly. The 
coordinate data is taken from the microscope’s digital readout and typed 
into a computer data file for processing. Iterative programs are used to 
generate a uniform array of “best fit.” The flow chart of Fig. 5 maps the 
main loops in the computer programs. 

Consider an array of measured data points (x,y) as compared to a 
perfectly uniform array (u,v). The uniform array is defined by a value 
for the uniform spacing of points along both coordinate axes (x,,7,), the 


* The dark spot is thought to be due to loss of germania or the inside surface of the 
preform which produces a dip in the index of refraction. 
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origin relative to the origin of the measured data (x,,y.) and the angle 
between the u-axis of the uniform array and the x-axis of the measured 
data (0). 

The measured data is transformed by translation and rotation into 
the coordinate system of the uniform array by the following transfor- 
mations. 


Xt = Ax + ayoy + by 


Ye = AgiX + Age + ba 


where 
Q11 = A92 = cos 8 
Tv 
aio =cos(—-—9 
= € ) 
vie 
ao, = cos (— + a) 
21 ce 
and 


by = —Xo cos 9 — yg sin 6 

bo = x, sin 8 — yo cos 0. 
The magnitude and angle of the vector, q;, joining the ith point of the 
uniform array with the corresponding ith transformed data point is 
calculated as simply ; 

la) = Gy =a)? + O40)? 

and 
Yt; — Vi 
Xt; — Uj 


t 


Z2q; = arc tan 


The parameters of the uniform array, 9, x9, Yo, X; and y, are now changed 
systematically according to the flow charts in Fig. 5 by using two separate 
computer programs. When this “number-crunching” is finished, the 
output gives the values of 0, x9, yo, x; and y; that define a uniform array 
of “best fit.” That is, the mean value of |q;| is a minimum for this uni- 
form array. This minimum mean value of |q;| is used as a “figure of 
merit” for the connector being characterized. The x and y components 
of q; (the offsets) are calculated and printed out along with the standard 
deviation of the minimum |q;| and the standard deviation of the x and 
y components of q;. 

Several additional calculations are made as a check on the measured 
data to eliminate gross errors. A histogram is calculated which gives the 
number of offsets, |q;|, between o to a, o to 2c, 20 to 3c, and greater than 
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30. Typically, a few offsets will be in the 20 to 30 interval and seldom any 
greater than 3c. Also the angle of the offset, Zq;, places it in one of four 
quadrants. The distribution of offsets among the four quadrants is 
printed out as a check to see if the offset vectors are more or less uni- 
formly distributed among the four quadrants. 

The primary use of these alignment calculations is for use as a “figure 
of merit” for an individual connector half. From the outset, it was realized 
that assembly and measurement of splice loss for a 144-fiber array splice 
would be tedious and time-consuming. This method of alignment 
characterization provides a basis for alignment comparison of connector 
halves, while requiring only one connector half to be assembled and 
without requiring splices to be assembled or losses to be measured. 


IV. ARRAY ALIGNMENT RESULTS 


Table IJ lists the alignment results for selected connector halves. A 
further description of these connector halves is needed. Connectors 34 
and 35 were the last 12 X 12 aluminum chip arrays! assembled. The 0.3 
mil to 0.4 mil value for the minimum mean offset is typical for 12 X 12 
aluminum arrays. (All offset data are in mils, 0.001 inch.) 

Silicon chip arrays have from the outset exhibited greatly reduced 
minimum mean offsets. Connectors 38 through 41 had mean offset values 
close to 0.1 mil. Connectors 40 and 41 were the input and output fanout 
connectors for the Atlanta Fiberguide System Experiment. Connectors 
45 and 46 were the connector halves on the fiberguide experiment cable. 
These connectors mated with fanout connectors 40 and 41 to provide 
access to the cable, and performed satisfactorily. Splice loss data for these 
specific splices are unavailable due to the fact that they appeared in 
series with other splices in the system configuration. 


V. SPLICE LOSS TEST CONNECTORS 


Connectors 47 through 50 were fabricated on the ends of two 0.8 meter 
long ribbon stacks for the purpose of measuring splice loss for these 
improved silicon chip arrays. Fiber diameter variations were measured 
to be approximately +2.5 percent for the fiber used in these connectors. 
Connectors 47 and 49 were epoxied at the same time in the holding vise. 
The stack height for connector 47 was less than that for connector 49 by 
0.4 mils and this “slop” may account for the slightly higher minimum 
mean offset. Connectors 48 and 50 were fabricated separately and gave 
improved offset characteristics. 


Vi. LOSS MEASUREMENT PROCEDURE 


Array splice losses were measured using the setup shown in Fig. 6. The 
beam from a Spectra Physics He-Ne laser was expanded and focused with 
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Connector 
number 


Array 
size 


12X12 
12X12 
12 X 12 
12X12 
12x 12 
12 X 12 
12X12 
12 X 12 
12X12 
12X12 
12X12 
12 X 12 


Material 


Table | — Connector data 


Min. 
mean 
offset, 

mils 


0.33 

0.39 

0.094 
0.104 
0.105 
0.104 
0.129 
0.114 
0.179 
0.135 
0.106 
0.095 


Var. 
of 
offset, 
mils 


0.20 
0.18 
0.05 
0.07 
0.08 
0.05 
0.06 
0.075 
0.072 
0.061 
0.051 
0.052 


Xs, 
mils 


9.03 

9.02 

9.02 

9.015 
9.011 
9.021 
9.021 
9.018 
9.001 
9.019 
9.017 
9.012 


Yss 
mils 


9.51 

9.80 
11.40 
11.33 
11.87 
11.85 
11.75 
11.804 
11.604 
11.658 
11.667 
11.591 


Xo, 
mils 


—0.05 
0.01 
~—0.01 
0.03 
0.09 
0.05 
—0.06 
—0.13 
0.03 
—0.01 
0.02 


Yo» 
mils 


—0.19 
—0.17 
—0.06 
—0.04 
0.03 
0.04 
—0.14 
0.024 
—0.05 
—0.03 
—0.13 
—0.01 


millirads 
—1.4 
—2.5 
—0.62 
—0.53 
1.55 
0.70 
0.44 
1.79 
—1.39 
—-1.12 
1.89 
0.68 
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Fig. 6—Loss measurement setup. 
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a 10X, 0.25 NA, lens onto the core of a fiber in the array. The reflected 
image from the fiber core was observed on the viewing screen (Fig. 6) and 
the laser spot positioned above the center of the core approximately 
midway to the core-clad boundary. The exact spot in this region was not 
critical since input power to the splice was found to be close to maximum 
as long as the spot was about midway between core center and clad- 
ding. 

Input power measurements were obtained by exciting each fiber of 
connector 47 and measuring the power exiting connector 48. The splice 
was then assembled by placing connector 49 (or connector 50) in position 
against connector 48. Again, each fiber of connector 47 was excited by 
positioning the spot as previously described and the power exiting con- 
nector 50 (or connector 49) was measured. Once again the same vicinity 
on the fiber core almost always corresponded to maximum transmission 
through the splices. This technique, by contrast with the technique of 
maximizing input and output power, allowed a larger time constant (1 
or 3 seconds) to be used on the lock-in amplifier, thereby reducing noise. 
fluctuations. The time required for a full 12 X 12 array loss measurement 
was reduced from 4 hours with the previous technique to less than 2 
hours with the new technique. The fact that only two negative losses 
occurred (—0.03 dB and —0.06 dB) compared to 130 positive losses be- 
tween 0 and +0.1 dB is indicative of measurement accuracy and re- 
peatability substantially better than 0.1 dB. 


Vil. LOSS MEASUREMENT RESULTS 
7.1 Case 1—Connector 48 to connector 49 

Figure 9 is a splice loss histogram for case 1. In this case, the ribbons 
were organized so that a particular fiber was spliced to its mating section. 
This splicing of “identical” fibers was done to minimize the effect of 
mismatch due to core diameter variations. 
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Fig. 712 X 12 array splice, connector 48 to connector 49. 


Three fibers were broken within connectors 47 or 48 as indicated by 
a low input power to the splice and one fiber was broken in connector 
49 or 50 as indicated by a large loss. These four fibers were not included 
in Fig. 7 due to definite scattering, indicating breaks. Two fibers also had 
low input power indicating a break in connector 47 or 48 although the 
losses for these splices were 1.05 dB and 1.29 dB respectively. It is sus- 
pected that another break occurred (1.26 dB loss) in either connector 
49 or 50 as this was the only other loss above 0.4 dB. 

These seven high losses are due to breaks and if they are excluded from 
the data, a mean loss of 0.18 dB results with a variance of 0.09 dB. Seven 
breaks in fabricating four connectors involving 576 fibers corresponds 
to 98.8 percent yield, and, although this is high, it is believed that im- 
provements in the yield can be made. 


7.2 Case 2—Connector 48 to rotated connector 49 


Figure 8 is a splice loss histogram for the case of rotating connector 
49 by 180° before joining it with connector 48. Fiber diameter mis- 
matches occur in this case and have contributed to the higher mean loss 
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Fig. 8—12 X 12 array splice, connector 48 to inverted connector 49. 


and variance. Of the five losses above 1 dB in Fig. 8, three can be at- 
tributed to breaks in fibers, leaving two losses over 1 dB unexplained thus 
far. 

Considering the possible effects of +2.5 percent core diameter varia- 
tion—if splice loss increases with the square of the ratio of the core di- 
ameters when the transmitting core is larger,+ then a worst-case effect 
would be 0.42 dB. This worst-case effect may account for the two losses 
over 1 dB mentioned earlier. Since no core diameter mismatch effect 
occurs when the transmitting core is the smaller, the mean effect is 
probably less than half the worst case or perhaps on the order of 0.1 dB. 
Thus the core diameter mismatch could account for most of the increased 
mean loss. 

7.3 Case 3—Connector 48 to connector 50 

Figure 9 is the splice loss histogram for this case. All losses ieee than 
1 dB have been previously identified as breaks. In this case, improved 
alignment accuracy appears to be present, since a significant portion of 
the mean loss could be due to diameter mismatch effects. It is unclear, 
at this time, how mismatch losses and alignment losses combine. From 
Table I it is seen that connector 50 has the least value of minimum mean 
offset and the least 6; therefore, improved splice loss is expected. 
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Fig. 9—-12 X 12 array splice, connector 48 to connector 50. 


It is also seen in Table I that y, for connector 50 is significantly less 
than y, for connector 48. This difference accumulates throughout the 
array and a connector height measurement shows that connector 48 is 
0.52 mils larger than connector 50. When the connectors are placed in 
the final alignment fixture, a greater force is applied to the thicker 
connector; however, it was unclear how much deformation occurred and 
to what extent connector height differences affect array splice losses. 
This case shows that 0.52 mils can be accommodated with little or no 
degradation in the splice. 


7.4 Case 4—Connector 48 to connector 49 


This case is a repeat of case 1. The connectors had been handled, as- 
sembled and disassembled several times since the first set of measure- 
ments. This test was to see if any noticeable deterioration had taken 
place and to check the repeatability of the final alignment procedure. 

Figure 10 shows that the mean loss is reduced slightly from case 1. All 
losses greater than 0.9 dB occur in positions which have been identified 
as broken fibers. It is interesting to note that the three losses between 
1 dB and 1.3 dB in case 1 are between 0.7 dB and 1.0 dB in this case. One 
fiber position measured 0.15 dB in case 1 and increased to 0.89 dB in case 
4 possibly due to contamination or end damage. 
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Fig. 10—12 X 12 array splice, connector 48 to connector 49. 


Vill. CORRELATION OF LOSS AND OFFSET 


A computer program has been written that calculates this mean value 
of offset between two sets of fiber position data. Two fiber positions that 
are well separated within the array are selected on the basis of low splice 
loss at these points. One position is used as the origin, the other is used 
to fix the position of the x-axis. This operation defines a coordinate 
system common to both arrays. Each set of fiber position data is trans- 
formed by the previously described set of transformations, and the mean 
offset is calculated as before. After generation of a computer data file 
containing the measured splice loss data, a correlation coefficient, p, is 
calculated. 


> lal -& — lal -b 


(2 lal? = nlail?) (2 - nl?) 
i=l i=1 
where |; = splice loss data and J; = mean splice loss. 
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Table II — Splice data 
Splices Probable [@;|, mils I;, dB p Loss limit, dB 
5 


1. 48 to 49 0.16-0.18 0.18 <0.2 0. 

2. 48 to inv (49) 0.18-0.20 0.32 <0.15 1.0 
3. 48 to 50 0.32 0.16 0 0.5 
4. 48 to 49 0.16-0.18 0.16 <0.2 0.5 


To eliminate the predominating effect of high-loss fibers, a loss limit 
is set to exclude high losses caused by fiber breaks. Table II lists the re- 
sults of the mean offset and correlation coefficient calculation. It is seen 
in Table II that the correlation coefficients are small, indicating little 
or no correlation between offset and loss on a fiber-by-fiber basis. This 
had not been the case with aluminum chip splices which exhibited cor- 
relation coefficients of typically 0.7. The reduced correlation coefficient 
for silicon chip arrays indicates that alignment accuracies and the re- 
sulting splice losses are approaching the limit of measurement resolution 
as applied on a fiber-by-fiber basis. 

Another approach is to compare the mean offset to the mean loss. 
Figure 11 is a measured loss versus offset curve® for a single fiber similar 
in index profile and diameter to the fibers used in these arrays. Nor- 
malizing a mean offset of 0.17 mils (from Table II for connector 48 spliced 
to 49) with respect to fiber core radius yields a mean offset of approxi- 
mately 0.16 core radius. From Fig. 11, this corresponds to a transmission 
of 96 percent or 0.18 dB loss. The measured mean loss was 0.18 dB (ex- 
cluding fiber breaks), in agreement with the value predicted from the 
measured mean offset. 

For connector 48 spliced to inverted connector 49, the mean offset was 
0.19 mils (0.176 core radius) which corresponds to 95 percent transmis- 
sion or 0.28 dB. The 0.04 dB by which the mean measured loss exceeded 
0.28 dB is probably attributable to core diameter mismatch. 

Connector 48 spliced to 50 does not follow the same trends as the 
previous two splices, due to the 0.52 mil overall connector height dif- 
ference and the supposed unequal deformation during final alignment. 
The 0.32 mil (0.30 core radius) mean offset should cause a 0.75 dB mean 
splice loss, from Fig. 11. The measured mean loss of 0.16 dB suggests that 
significantly unequal array deformation is taking place in this splice as 
previously supposed. 


IX. CONCLUSIONS 


Etched silicon chips have enabled uniform arrays to be fabricated with 
an average error of 0.1 mil without complex tools or microscopes. Ex- 
ceptionally low mean splice loss data have been obtained which corre- 
lated well with mean offset. 
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Fig. 11—Measured transmission vs. transverse offset. 


Fiber diameter variations of 42.5 percent caused an increase in the 
mean loss on the order of 0.1 dB; however, it is not known to what extent 
array uniformity was affected by these fiber diameter variations. Two 
arrays with an overall height difference of 0.52 mil had very low mean 
splice loss, indicating that this height difference was equalized within 
the splice holding fixture. After several combinations of arrays were 
assembled and loss measurements made, the original array configuration 
was measured to test for array deterioration, loss measurement repeat- 
ability and final alignment repeatability. Only one fiber position showed 
evidence of contamination, and repeatability was within 0.02 dB for the 
mean. 
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Accurate Silicon Spacer Chips for an Optical- 
Fiber Cable Connector 


By C. M. SCHROEDER 
(Manuscript submitted July 1, 1977) 


Assembled silicon array connectors have been fabricated with low 
splice loss. The silicon chips can be manufactured with submicrometer 
repeatability and accuracy. Improvements in chip thickness variations 
are anticipated which in turn should reduce the overall groove opening 
variations to a point where individual chip selection is not required. 


I. INTRODUCTION 


Low-loss multifiber splices have been obtained with the use of accurate 
silicon spacer chips. These spacer chips are used in a unique splicing 
technique developed by Bell Laboratories! and used in the Atlanta Fi- 
berguide experiment. This multifiber cable splice is a stacked array 
consisting of two properly prepared cable connectors butted end to 
end. . 

A cable array connector is a laminate sandwich of silicon alignment 
chips with “V” grooves on both top and bottom surfaces interleaved with 
optical fibers epoxied to form a two-dimensional array (Fig. 1). A com- 
pleted connector may consist of up to 144 optical fibers positioned by 
up to 13 alignment chips. Each fiber in the array must be positioned 
accurately to its counterpart in a mating array to obtain low splice loss. 
The transverse misalignment of the mating fibers should be no greater 
than one-tenth the core radius or typically 2.5 micrometers.” This can 
only be obtained with alignment chips having high dimensional accuracy 
with respect to thickness, groove geometry, and position. 

The feasibility of the multifiber splice was shown with alignment chips 
manufactured of aluminum. The aluminum chip, however, could not be 
manufactured repeatedly with the high dimensional accuracies required. 
An improved chip was developed, that met the above goals, using (100) 
oriented silicon* and photolithographic techniques. This improved chip 
has demonstrated a significant reduction in the overall array splice 
loss. 
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Fig. 1—Schematic of a multifiber array connector. 


Il. PHOTOLITHOGRAPHIC AND PROCESSING TECHNIQUES 


The use of silicon for alignment chips offers several advantages. The 
material is readily available, it can be handled and processed relatively 
easily, and it can produce precise V grooves to photomask accuracy. 

When a (100) surface slice of silicon covered by an oxide mask is 
submerged in a basic solution, the etch rate is much greater in the (100) 
direction than in the (111) direction (Fig. 2). This anisotropic etching 
of the silicon results in a precise V groove to an angle of 70.53° with the 
reaction self-stopping at the point where the (111) planes intersect. The 
groove opening is determined by the opening in the oxide mask. The 
(110) plane is used as a reference to align the photomask features parallel 
to the intersection of the (100), (111) planes. 

The processing of the silicon first starts with the selection of a 5 cm 
diameter (100) surface oriented wafer polished on both surfaces to a 
thickness of 0.25 mm. Thickness variations and surface defects should 





Fig. 2—Crystallographic planes used to manufacture alignment chips. 
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Fig. 3—Processing techniques. 


be minimized for best results. The wafer is cleaned, then placed into a 
tube furnace for a 1-micrometer oxide growth. The wafer is then coated 
on both surfaces with a positive photoresist, placed between two preal- 
igned photomasks and exposed to ultraviolet radiation (Fig. 3). The 
exposed wafer is removed, then developed leaving open windows of sil- 
icon oxide which are in turn etched in a buffered hydrofluoric acid so- 
lution. The remaining photoresist is removed leaving a silicon wafer 
covered by an oxide mask which ‘is then placed into a potassium hy- 
droxide solution for preferential etching®. A cross section of a typical 
silicon alignment chip is shown in Fig. 4. 

























































| 
| 
| 
| 
0.2286 mm 9.52 mm 
0.254 mm = 
eel | 
| 
| 
ee 
woojpoooong 12 GROOVES 





Fig. 4—Cross section of typical alignment chip with dimensions. 
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Fig. 5—Data on groove opening variations. 


tl. ALIGNMENT CHIP EVALUATION 


A listing of variations that can cause splice loss in a multifiber con- 
nector in their present order of importance is as follows: 


(¢) Alignment chip thickness variations 

(it) Fiber diameter variations 

(iit) Groove opening variations 

(tv) Transverse misalignment of the top and bottom grooves 
(v) Aperiodic grooves 


Excluding fiber diameter variations, the remaining list must be con- 
trolled by alignment chip accuracy. 

Measurements on the photomask features and the chip profile were 
made with the use of a Hewlett Packard model 5526A laser interfer- 
ometer. The object whose features are to be measured is placed onto a 
traversing stage under the cross hair of a microscope. A reflecting cube 
attached to the stage monitors the movement of one of the two paths of 
the interferometer. This movement is compared to a fixed length of the 
second path with the changes displayed on a monitor. The repeatability 
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Fig. 6—Splice loss as a function of offset. 


of the data was found to be 0.2 micrometer, limited chiefly by the vernier 
acuity of the eye. 

Measurements on the photomask features indicate the required di- 
mensions were well within the 1 micrometer stated accuracy for the 
masks with line-width variations not measurable across the mask face. 
Profile measurements on an etched silicon wafer indicate the groove 
periodicity of 228.6 micrometers was maintained within the measure- 
ment error. The transverse misalignment of the top and bottom grooves 
was held to 1 micrometer or less and the groove depth maintained at 50 
micrometers. Variations of this type do not significantly add to the 
overall splice loss. The major contributors to array splice loss other than 
fiber diameter variations, was found to be chip thickness variations and 
groove opening variations. 

The groove opening directly affects the vertical position of the optical 
fiber. The groove angle is constant but relatively steep with 1 unit of 
groove opening corresponding to 1.4 units of vertical drop in the optical 
fiber. The data on groove opening variations are shown in Fig. 5. These 
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Fig. 7—Laser-scribed silicon wafer. 


data were taken on both surfaces of one silicon wafer with the cross- 
hatched region being data from surface 1 superimposed onto data from 
surface 2. An overall shift in data of 0.2 micrometers can be observed 
when comparing the two surfaces. This shift was caused primarily by 
unequal exposure times. Maximum groove width variations were found 
to be +0.7 micrometers with a standard deviation of 0.23 micrometers. 
The maximum difference in groove width of 1.4 micrometers would cause 
vertical positioning errors of the optical fibers of 1.96 micrometers for 
a worst-case situation. While this deviation does not exceed the 0.1 core 
radius offset, it does point to a need to minimize groove width variations, 
which are probably caused by variations in wafer thickness with Fresnel 
diffraction playing a part in undercutting at the areas where the wafer 
and photomask did not contact. Thickness variations are at present the 
largest potential contributor to splice loss with total variations of 4 mi- 
crometers encountered on some wafers. 


IV. LOSS DUE TO TRANSVERSE OFFSETS 


Splice loss as a function of offset for small offsets of parabolic index 
fibers is shown in Fig. 6.° K is the ratio of the radius of the receiving fiber 
to the radius of the transmitting fiber. The maximum offsets due to 
groove opening and thickness variations are displayed to indicate the 
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probable loss that could be expected in the worst-case situation. To avoid 
variations of this type an individual chip selection is required at 
present. 

Splice loss data for 144-fiber array connectors using silicon alignment 
chips have been made with the results presented in a companion report 
by C. M. Miller. 


V. SUMMARY AND CONCLUSION 


A laser scribed silicon wafer is shown in Fig. 7. Each 5-cm diameter 
wafer yields 36 usable alignment chips 3.8 X 9.52 mm in size. 

Assembled silicon array connectors have been fabricated with low | 
splice loss. The silicon chips can be manufactured with submicrometer . 
repeatability and accuracy. Improvements in chip thickness variations 
are anticipated which in turn should reduce the overall groove opening 
variations to a point where individual chip selection would not be re- 
quired. 
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On the Phase of the Modulation Transfer 
Function of a Multimode Optical-Fiber Guide 


By I. W. SANDBERG, I. P. KAMINOW, 
L. G. COHEN, and W. L. MAMMEL 


(Manuscript received June 1, 1977) 


We consider the range of validity of a Hilbert-transform approach 
in which the measured magnitude of the modulation-transfer-function 
of an optical fiber is used to compute the fiber’s impulse response. It 
is argued that a key “minimum-phase assumption” can fail to be sat- 
isfied in important cases, and a few closely related experimental and 
analytical results are presented. 


Il. INTRODUCTION 


Pulse dispersion in an optical fiber transmission line limits its infor- 
mation-carrying capacity by limiting the temporal spacing of input 
pulses that can be resolved at the output. The impulse response g(t), by 
which we mean the output power of a fiber excited by a unit impulse of 
optical power, provides the necessary information concerning the dis- 
tortion of the pulse by modal and material dispersion. For a strictly 
monochromatic pulse source, only modal dispersion contributes to the 
distortion. However, with regard to the corresponding measurement 
problems, it is difficult to obtain sufficiently short (<0.5 nsec) and 
monochromatic (< 10 A) input pulses to accurately study fibers with very 
low modal dispersion. Another difficulty is the lack of availability of 
suitable sources that are tunable over a wide range of wavelengths, in- 
cluding wavelengths longer than 1 um which are of interest for practical 
fiber systems. 

One can obtain g(t) from the modulation frequency transfer function 
(MTF) G(w), which is the envelope response of the fiber to an incoherent 
optical signal sinusoidally modulated in amplitude at angular frequency 
w. Personick! has shown that to the extent that certain reasonable ap- 
proximations hold, g(t) and G(w) are a Fourier transform pair. In prin- 
ciple, the MTF can be determined experimentally. The method employed 
in Refs. 2 and 3 uses a xenon lamp and monochromator as a tunable 
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source that can be sinusoidally modulated to high frequency (>1 GHz) 
by an electrooptic modulator. It is straightforward to measure the 
magnitude |G(w)| of the transfer function 


G(w) = |G(w) Jet) (1) 


using available components, including an RF spectrum analyzer. How- 
ever, because the fibers must be long (~1 km) in order to obtain good 
measurement precision, the measurement of 0(w) appears to be formi- 
dable. As f = w/2z varies from zero to 1 GHz, 6(w) varies nearly linearly 
with f from zero to 104m radians for a 1 km long fiber. However, the 
contribution of #(w) to pulse dispersion is due to a nonlinear deviation 
A6(w), on the order of 27, from the much larger linear phase shift 09(w); 
i.e., 


Ww) = Ao(w) + AP(w) (2) 
O9(w) = w/v (3) 


where L is the fiber length and v is an effective envelope velocity taken 
to be independent of w. Hence, direct measurement of phase distortion 
in the presence of the large frequency-dependent 09(w) could be subject 
to large error as v varies with temperature or other environmental factors. 
One is therefore led2:3 to consider methods of mathematical computation 
of #(w) from |G(w)| using Hilbert-transform theory. 

The main purpose of this note is to report on results which indicate 
that unfortunately the Hilbert-transform approach described in Refs. 
2 and 3 is, in general, not a helpful one for the particular problem at issue, 
even though early experimental results suggested otherwise.t 

Roughly speaking, it is known that by using a Hilbert-transform 
relation the phase can be obtained from the magnitude of a transfer 
function provided that the transfer function is “minimum phase.” A 
standard condition (which is by no means sufficient) for “minimum- 
phase behavior” is that the Laplace transform of the impulse response 
that corresponds to the transfer function have no zeros in the closed right 
half-plane. This is in accord with the observation that for an ideal 
waveguide with constant positive delay 7) and transfer function Go = 
e—’70, the phase cannot be determined from a knowledge of the function 
|Go(w)| alone, but a transfer function that has the same magnitude as 


that of the waveguide is Gj(w) = 1 for which the phase is zero for all 
w. tt 


1 The statement on page 1518 of Ref. 2 concerning the possible lack of “approximate 
ua phase behavior” was motivated by the results of the joint work described 
ere. 
tt The mathematical reason that Go is not a “minimum-phase” function (even though 
e~70 has no zeros for Re(z) = 0) is that In|e—?7! fails to satisfy a sufficiently strong growth 
condition in the half-plane Re(z) = 0. See Sections 5.1 and 5.3 for related material. 


100 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1978 


For our purposes, the difference between Go(w) and G1(w) is unim- 
portant, because we are willing to ignore a constant time delay that can 
easily be estimated. More generally, it is reasonable to ask if e#G(w) 
is a minimum-phase function, where 79 denotes the linear part of the 
delay. An early impulse response measurement on an actual fiber sug- 
gested? that this might indeed be the case. However, the further ana- 
lytical and experimental study reported on here shows that nonmini- 
mum_-phase behavior is likely to arise, and can arise, in important actual 
cases. Some additional closely related material is also presented. 

Methods for circumventing the difficulties described in this note are 
under study, and it is expected that they will be described in a later 


paper. 


ll. SOME ANALYTICAL PROPERTIES OF THE MULTIMODE TRANSFER 
FUNCTION 


In the general case, the transfer function of a fiber can be written in 
the form 


Ge) = Sf, psn (4) 


in which the integrator a(r) is a real-valued monotonically nondecreasing 
function of r,t and T, and Ty, which are fixed by the refractive indices 
of core and cladding, are the smallest and largest modal delays, respec- 
tively.t? Often a(r) is normalized so that 


To 
f days 
Ta 


For a fiber that can propagate n discrete modes without mode mixing, 
(4) becomes 


G(w) = >> dje “ti (5) 
J=1 


in which each d; is a positive constant that represents the initial exci- 
tation of the jth mode. Typically, n > 100. Most of our discussion is 
concerned with the important particular case in which (5) holds, and, 
in order to avoid a lack of continuity of the presentation, proofs of the 
results discussed are given in a separate section. We assume that the 7; 
are ordered so that 7; < 79 <... < Tp. 

In (5), G(w) is the generalized Fourier transform of a finite train of 


* Thus, roughly speaking, da(z) in (4) can be replaced with b(7)d7 in which the function 
b(r) is nonnegative and may contain impulses corresponding to discrete modes. See Refs. 
1 and 2 for the relevant background material. 

+tWe are of course assuming that material dispersion can be neglected. 


MULTIMODE OPTICAL-FIBER GUIDE 101 


not-necessarily-equally-spaced impulses. Let H(z) denote the corre- 
sponding Laplace transform. That is, let H(z) denote 


nm 
—go72T] 
» dje~?7i 


j 
for all complex z. Of course G(w) = H(iw) for all w. 

A standard Hilbert-transform method*2° for determining the phase 
6(w) of a transfer function G(w), from the function |G(w)|, when it is 
possible to do so, is to use the formulat 


i) = 6° REO A, (6) 


y? — w? 
which, roughly speaking, amounts to a direct application of the Hilbert 
transform 


Imlftio)] = 72 (° Re (7) 


y2 — w2 


in which f(z) is any complex-valued function of z that satisfies certain 
conditions (such as those described in Section 4.1tt) which include the 
condition that f(z) is analytic for Re(z) = 0. 

The way in which (7) is used to obtain (6) is of course to let f(z) bea 
suitably defined single-valued variant of In[A(z)], in which H(z) is the 
Laplace transform of the time function whose Fourier transform is G(w). 
That is what gives rise to the well-known requirement that H(z) be 
zero-free in the closed right half-plane. 


2.1 The zeros of H(z), and related material 


With regard to the location of the zeros of H(z), according to Propo- 
sition 2 of Section IV: H(z) + 0 for Re(z) = Oif 


di> > dj (8) 
j=2 


and if (8) is not satisfied, then given an arbitrarily small positive number 
«, and any set of n numbers t; < to<...<tn, we have H(z) = 0 for some 
z in the closed right half-plane for a choice of the 7; such that |7; — ¢;| 
< ¢ for allj. In particular, and roughly speaking, unless (8) is satisfied, 
given any set of n delays, there is a set of delays arbitrarily close to those 


t The integrals in (6) and (7) are to be interpreted as Cauchy principal values. 

ttSection 4.1 contains an outline of a proof that (7) holds under certain specific condi- 
tions. The derivation given, for instance in Ref. 4, lacks rigor in that, for example, the point 
s in Ref. 4 is initially assumed to be a point internal to a certain contour, while subsequently 
an expression based on that assumption is evaluated for s on the contour. The main reason 
for including the basically tutorial material of Section 4.1 is that it is used to prove another 
result described in this note. 
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delays such that H(z) is not “minimum phase.” Notice that it is not 
claimed that H(z) has a zero in the closed right half-plane whenever (8) 
is violated.t However, Proposition 2 does imply that whenever (8) is not 
satisfied it is incorrect to assert that H(z) ¥ 0 for Re(z) = 0 when the 
7; are known only to within some positive tolerance «, no matter how 
small ¢ is. Therefore, H(z) is zero-free in the closed right half-plane, and 
that property is structurally stable in the sense indicated, if and only if 
(8) is met. This result suggests that it would not be surprising to en- 
counter “nonminimum-phase” behavior with fibers for which the total 
power in a sufficient number of the modes corresponding to the delays 
T9, T3,.+-, Tn 1S considerable. 

An idealized example in which a somewhat analogous conclusion is 
reached is as follows. Consider an n-mode fiber without mode mixing 
for which the modal delays are equally spaced by 6 sec, so that 7; = 7’o 
+ 76 for each positive integer 7. Assume that the fractional power into 
the jth mode is given by d; = ce¥ for all j, in which y and ¢ are constants 
with c > 0. Then 


, e(y-zd)n =e | 
Hz) = cere = G=28) 
and the condition that H(z) + 0 for Re(z) = 0 will be satisfied if and only 
if y< 0.t+ Of course y < 0 means that modes with larger delays have 
smaller excitation. A similar conclusion is reached for the continuum 
mode-mixing case? in which the integrator a(7) of (4) has a continuous 
derivative that is proportional to e7’. 

While the discussion in the preceding paragraphs suggests!! that there 
are important cases in which (6) cannot be used, it certainly does not rule 
out the possibility that there is some other method for determining 6(w) 
from |G(w)| (which, for example, might Bossy exploit the fact that 
the d; are positive). 

In this connection, consider (5). In order to avoid the necessity of in- 
troducing a function equal to e’“'1G(w), assume throughout the re- 
mainder of this section that 7; = 0 (which of course simply provides a 
normalization‘). 

Suppose, for example, that n = 2 and that dz = (1 — d;) with 0 < dy 
<1and d, # 1/2. Then |G(w)| = d;? + (1 — dy)? + 2d,(1 — dj) cos (wr2) 


t In fact, we show that that claim would be false. 

tt Since (8) is violated when 7 is negative and sufficiently close to zero, we see that a given 
specific H(z) can be zero-free in the closed right half-plane when (8) is violated. 

t This example was suggested by H. E. Rowe. 

tt Little information is available concerning how to accurately specify G(w) for an actual 
fiber using purely analytical methods. Very small geometrical perturbations can have 
significant effects on the impulse response of graded-index fibers.® 

+t It is easy to see that without some such normalization, it is is not possible to determine 
6(w) from |G(w)|. 
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Fig. 1—Measured (deconvolved) impulse responses for (a) a GeO2-SiOz fiber with a < 
Qopt; (b) a BgO3-SiOg fiber with a ~ adopt. 


which clearly is unchanged if d, is replaced with (1 — d1). Hence given 
only |G(w)| for all w, and that (5) with 7; = 0 holds, it is not possible to 
find a unique 4(w). 

With regard to results in the opposite direction, when (5) holds and 
(8) is satisfied, it is true that the phase function 6(w) can be obtained from 
|G(w)| by using (6). This is proved in Section 5.3. 


lil. EXPERIMENT 


For each of two fibers A and B the impulse response was measured by 
injecting 0.4 nsec pulses (2c width) from a GaAs laser (A = 0.9 wm) and 
observing the pulse distortions after propagation through the fibers. 
Fiber A was 1 km long and had a graded index GeO2-SiOz2 core with a 
~ 1.9 and appt (0.9 um) ~ 2.0. Fiber B was 1 km long and had a graded 
index BzO3-SiOz core with a © aopt(0.9 um) ~ 1.8. The measured impulse 
responses are shown in Figs. 1a and b. These impulse-response functions 
were Fourier transformed to obtain |G(w)| for each case. Then each 
|G(w)] together with its phase calculated from piecewise-linear formulas4 
based on (6) was used to calculate a corresponding impulse response. The 
plots are shown in Figs. 2a and b for fibers A and B, respectively. It may 
be seen that the calculated and measured responses agree quite well for 
fiber B with a suitable normalization and translation to bring them into 


+ Analogous examples can be given which hold for all n = 2. For instance, let H(z, y) 
denote the expression for H(z) for the idealized exponential-excitation case mentioned 
above, with c chosen to depend on ¥ so that H(0, y) = 1. It is not difficult to verify that 
we have |H(iw, y)| = |H(iw, —y)| for all w for each y. 
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Fig. 2—Impulse responses calculated from the magnitudes of the fourier transforms 
of the measured impulse responses in Figs. 1(a) and (b). 


coincidence. On the other hand, the measured and calculated responses 
for fiber A do not agree.? 

The fact that the functions of Figs. la and 2a appear to be approximate 
mirror images of one another suggests that all of the dominant zeros 
associated with the transfer function of fiber A might lie in the right 
half-plane. In this connection, we note that a general function H(z) of 
the form defined in Section II can have zeros in both half-planes. For 
example, with H(z, y) as defined in a preceding footnote, the product 
H(z, 1)H(z, —1), which can be written in the same form as H(z), has 
zeros in both half planes. 


IV. CONCLUSIONS 


With regard to the overall problem of determining the impulse re- 
sponse g(t), direct measurement of the phase appears to be difficult and 
the general use of the “minimum phase” assumption to calculate the 
phase does not appear to be justified. 

Methods for circumventing the difficulties described are under study, 
and it is expected that they will be described in a later paper. 


V. APPENDIX: PROOFS 
5.1 Outline of a derivation of a well-known formula 


Let z be a complex variable with real and imaginary parts x and y, 
respectively, and let z* denote the complex conjugate of z. Let f bea 
t An early computational error led to a reversal of the sign of the time scale for the phase 


calculations as reported in Refs. 2 and 3. Thus, it was erroneously reported that g(t) for 
measured and calculated responses matched for fiber A. 
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Fig. 3—Contour in the (x, y) plane. 


complex-valued function of z defined throughout an open subset S of 
the (x, y)-plane that contains the half-plane x = 0 such that: f(z*) = 
f(z)*, Re[f(iy)] is an even function of y, and Im|f(iy)] is an odd function 
of y. 


Proposition 1. Suppose that f is analytic on S, and that |f(z)/z| — 0as 
|z| — ~ in the half-plane x = 0. Then 


Im[f(iw)] = aa = gn ees Ree 


for each real w > 0, in which P denotes a Cauchy principal value. 


Outline of a Proof. Assume that the hypotheses of Proposition 1 are 
satisfied. Let w > 0 be given, and choose « > 0 such that ¢ < w and f is 
analytic for |z — iw| < 2¢. With p any positive number such that p > w 
+ ¢, let C denote the contour shown in Fig. 3 which consists of a semi- 
circular arc of radius p, two semicircular arcs of radius e, and a portion 
of the y-axis. 

By Cauchy’s integral theorem, 


fF) 


Cz—le 


fiw Jas Se 


106. THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1978 


Therefore, 


f(z) 
I eal a 10 
myGe)|'= Qn JC (z — iw)(z + tw) a aD 
Consider separately the following contributions to the right side of (10): 
I, the integral along the y-axis from z = ip to z = —ip excluding the two 


gaps due to the semicircular ¢ arcs, the sum I of the integrals along the 
e arcs, and the integral J3 along the remaining arc of radius p. 
We find at once that 


LL + Se 


Using the fact that the integral of (z — iw)—! over the upper « arc is iz, 
it follows that Ig = % Im[f(iw)| + 6(€) in which 6(€) > 0 as « > 0. Also, 
I3 = B(p) in which B(p) > 0 as p > ~. 
Since Im[f(iw)] = J, + Io + I3, we have 
Im[f(iw)] = 22, + 26(e) + 2B(p) 


from which it is clear that the limit 
lim lim 2/, 
pro «0 


exists and is equal to Im|/(iw)]. This completes the outline of the 
proof.? 


5.2 Result concerning the zeros of H(z) 
Let H(z) denotett 


Ae . 
d dje~*7i 
Ja1 


for all complex z, in which n 2 2, d; > 0 for all j, and the 7; are real 
numbers such that 7; < tg <... < Tp. 

Let f4, to... ,t, denote any set of n real numbers with the property 
that ty <te<... <th. 

Consider the condition 


as > (11) 
J=2 


Proposition 2. We have H(z) # 0 for Re(z) = 0 if (11) holds. If (11) is 
not satisfied, then for any positive ¢ there is a choice of the 7; such that 
|r; —t;| < for all j and H(z) = 0 for some z with Re(z) = 0. 


t For further results concerning Hilbert transforms, see, for example, Ref. 6. 
tt For the reader’s convenience the definition of H(z) is repeated here. 
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Proof. If (11) holds and z = x + ty with x = 0, 


|H(z)| = 





n 
e727 E +y djen*i-70 | 
j=2 





> ent (41 25 d;) >0 


j=2 
Suppose now that 


di< d; (12) 


iM 


and let « be given. Let z = x + iy with y = z/e. Choose 7; = t; — ¢, and for 
each J = 2,3,...,m, choose 7; such that e~(j-7™) = — 1 and |7; — ¢;| 
< «. We have 


H(z) = e7271 | a: - 3 djestind | (13) 
jJ=2 


in which 7; — 7; > 0 fory 2 2. Using (12) and (13), it is clear that there 
is an x = Osuch that H(z) = 0, which completes the proof. 


5.3 A corollary of proposition 1 


Concerning the function H defined in Section 5.2 consider the fol- 
lowing hypothesis. 


Hypothesis: 7, = 0 and (11) is satisfied. 

We notice that if the hypothesis holds, then, for each real w, Re[H (iw)] 
> 0 and we have H(iw) = |H(iw)|e'*), in which ¢{iw) denotes the 
' principal value of tan~! {Im|H (iw)]/Re[H(iw)]}. 


Proposition 3. If the aes holds, then 
In|H 
liv) =p MLR ay 
y?2 — w? 
for w > 0 in which P denotes a Cauchy principal value. 
Proof: Assume that the hypothesis is satisifed. We shall show that 


Proposition 1 can be used. 
Let xo be a negative number such that 


n 
d>> dje~*0rj 
j=2 


Wesee that there is a positive constant 6 such that Re|H(z)] > 6 for Re(z) 
> xo. Let 6(z) denote the principal value of tan~! {Im[H(z)]/Re[H (z)]} 
for Re(z) > xo. 

For each complex z such that z ~ 0, let p(z) denote the principal value 
of In z. Thus, p[H(z)] = In| H(z)| + i¢(z) for Re(z) > xo, and p[H(z)] is 
analytic in z for Re(z) > xo. 
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Reduction of Network States Under Symmetries 


By V. E. BENES 
(Manuscript received May 24, 1977) 


It is a folk-theorem of traffic theory that if all sources have the same 
stochastic behavior, then symmetries of a telephone connecting network 
can be used to lump together equivalent states and to reduce the 
number of equations to be solved for the state probabilities. The 
structural and algebraic bases of this idea, and its connections to sto- 
chastic models, are studied here by means of concepts from lattice 
theory, group theory, and combinatorics generally. It is shown that 
when offered traffic is homogeneous and routing is structurally con- 
sistent, the state equations for certain natural Markov processes 
(representing operating telephone networks) can be substantially 
simplified by restricting attention to “macrostates,” defined as the 
structural equivalence classes of states, of which there are typically 
many fewer than of states. Reduced state equations are then obtained 
for general networks under simple Markovian traffic assumptions. 


I. INTRODUCTION 


Most telephone connecting networks are built in stages of identical 
units, arranged symmetrically in frames, and joined by symmetric 
cross-connect fields (Fig. 1). As a result, their structure has so much 
symmetry that it is possible quickly to identify at least some network 
states that are “essentially” equivalent in that their combinatorial 
structure is the same, and that they differ only in point of renaming 
terminals, switches, links, customers, etc. It has long been known in the 
informal lore of traffic theory that if the traffic sources at the terminals 
have the same stochastic behavior, then these symmetries of the network 
could be used as a basis for lumping together equivalent states and re- 
ducing the number of equations to be solved for the state probabili- 
ties. 

Such a line of thought was first pursued in a formal way by S. P. Lloyd! 
in unpublished work dated September 19, 1955, and we follow it further 
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here. But whereas Lloyd right away considered probability transition 
rate matrices (for Markov processes) which he assumed admitted sym- 
metry operations, we shall instead first relate the relevant symmetries 
to the network graph and the semilattice of states, without reference to 
a probabilistic model, and only later consider a natural traffic model. 
Our approach remains combinatorial as long as possible, and allows us 
eventually to include the effects of routing decisions, to connect the re- 
duction ideas with optimal routing, and to find that certain natural 
transition matrices necessarily admit the network symmetries. In par- 
ticular, we show that a traffic model used in previous work? has a tran- 
sition rate matrix that admits the symmetries of the network graph, 
provided only that the routing matrix used also admits these sym- 
metries. 

An additional practical incentive for the present study is the fact that 
for small networks, such as concentrators, it is often possible to press 
the advantages gained by reducing the states to their equivalence classes, 
and to solve completely the problem of optimal routing. One can then 
devise a circuit or a finite-state machine to mechanize the optimal 
routing policy, as has been done by A. F. Bulfer for the RTA concentrator 
structure.?4 Indeed, historically, it was our attempts to prove some of 
Bulfer’s surprising empirical results that led to a realization that a 
thorough study of structural equivalence of states was valuable and 
necessary.°6 


il. SUMMARY 


Various preliminaries are in Sections III and IV: a model for discussing 
connecting networks, with an account of the role of symmetry. Prior 
results of S. P. Lloyd on the use of groups to reduce state equations are 
described in Section V, and there is a heuristic discussion of some nec- 
essary conditions on such groups in Section VI: they should be groups 
of automorphisms of the semilattice of states that preserve the relation 
of having the same calls up. The symmetry group G, of the network 
graph appears in Section VII; it is used in Section VIII to define struc- 
tural equivalence of states and its associated group G,, and in Sections 
IX and X to define a natural homomorphic image G, of G, into G, which 
usually gives a more economical description of equivalence in terms of 
a group than does G,. The semilattice of states induces a partial ordering 
of the reduced states, described in Section XI. 

The next four sections are devoted to detailed calculations of sym- 
metry groups Gp, G,, and G, for various well-known networks: individual 
2 X 2 switches in Section XII, a random slip concentrator in Section XIII, 
a small three-stage network in Section XIV, and general crossbar net- 
works in Section XV, including frames, cascades of stages, and the No. 
5 crossbar type network. The final four sections consider a simple sto- 
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chastic model for traffic in a network with routing decisions, and show 
how the reduced state equations are derived in this setting. 

It can be concluded that symmetry groups afford a precise definition 
of structural equivalence for network states; this equivalence in turn 
allows a substantial reduction in the number of equilibrium equations 
for state probabilities in suitable stochastic models, thus extending the 
range of computable examples. For modest networks the reduction 
method can be used to perform optimal routing in explicit ways.>® 


Il. PRELIMINARIES 


We shall use a model for the structural and combinatorial aspects of 
a connecting network. This model arises by considering the network 
structure to be given by a graph G whose vertices are the terminals of 
the network, and whose edges represent crosspoints between terminals 
by pairs, with some of the terminals designated as inlets or outlets. Calls 
in the network are described by paths on G from an inlet to an outlet. 
Thus a connecting network v is a quadruple v = (G,J,Q,S) where G isa 
graph depicting network structure, J is the set of vertices of G which are 
inlets, 2 is the set of outlets, and S is the set of permitted or physically 
meaningful states. It is possible that J = Q (one-sided network), that I 
() Q = ¢ (two-sided network), or that some intermediate condition ob- 
tain, depending on the “community of interest” aspects of the network 
v. Variables w, x, y, and z at the end of the alphabet denote states, while 
u and v denote a typical inlet and a typical outlet, respectively. 

A possible state x can be thought of as a set of disjoint chains on G, 
each joining J to 2. Not every such set of chains need represent a state 
in S: wastefully circuitous chains may be excluded from S. The set S is 
partially ordered by inclusions <, where x < y means that state x can 
be obtained from state y by removing zero or more calls. It is reasonable 
that if y is a state and x results from y by removal of some chains then 
x should be a state too; i.e., S should be closed under “hangups.” It can 
be seen from this requirement that the set S of permitted states has the 
structure of a semilattice, that is, a partially ordered system whose order 
relation is definable in terms of a binary operation () that is idempotent, 
commutative, and associative, by the formulax < y iffx =x (1) y. Here. 
for x () y we can simply use literal set intersection: x () y is exactly the 
state consisting of those calls and their respective routes which are 
common to x and y. 

An assignment is a specification of what inlets are to be connected to 
what outlets. The set A of assignments can be represented as the set of 
all fixed-point-free correspondences from subsets of J to Q. The as- 
signments form a semilattice in the same way that the states do, and A 
is related to S as follows: call two states x,y in S equivalent as to as- 
signment, written x ~ y, iff all and only those inlets uel are connected 
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in x to outlets veQ which are connected to the same vu in y, though pos- 
sibly by different routes; the realizable assignments can then be iden- 
tified with the equivalence classes of states under ~, and there is a nat- 
ural map y: S — A, the projection that carries each state x into the 
assignment (x) it realizes, i.e., the equivalence class it belongs to 
under ~. 

With x and y states such that x = y, it is convenient to use x-y to mean 
the state resulting from x by removing from x all the calls in y. Similarly, 
with a and b assignments such that a = b, we use a-b to mean the as- 

_signment resulting from a by dropping all the connections intended in 
b. Note that here x-y, a-b have their usual set-theoretic meaning. 

It can now be seen that the map y is a semilattice homomorphism of 

S into A, with the properties: 


x zy y(x) 2 v(y) 
x2zy=y(x —y) = v(x) -— yy) 
v(x Oy) S$ (x) ON vy) 


y(x) = ¢—x = 0 = zero state, with no calls up 


Not every assignment need be realizable by some state of S. Indeed, 
it is common for practical networks to realize only a vanishing fraction 
of the possible assignments, and the networks that do realize every as- 
signment, the so-called rearrangeable networks, have been the objects 
of substantial theoretical study. Thus the image set y(S) of realizable 
assignments is typically much smaller than the set A it is embedded in. 
A unit assignment is, naturally, one that assigns exactly one outlet to 
some inlet, and it corresponds to having just one call in progress. It is 
convenient to identify calls c and unit assignments, and to write y(x) U 
c for the larger assignment consisting of y(x) and the call c together, with 
the understanding of course that c is “new in x” in the sense that neither 
of its terminals is busy in x. 

We denote by A, the set of states that are immediately above x in the 
partial ordering < of S, and by B, the set of those that are immediately 
below. Thus 


A, = {states reachable from x by adding a call} 
B,, = {states reachable from x by a hangup} 


For c new in x, let Acy = Ax 1) y~[y(x) U c]; Acx is the subset of states 
of A, that could result from x by putting up the call c, because y~!y(y) 
is precisely the equivalence class of y under ~. If A,; is empty then we 
say c is blocked in x: there is no yeA, which realizes the larger assignment 
v(x) U c. It can be seen that with F,, the set of new calls of x that are not 
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blocked, the family {A,,,ceF,} forms the partition of A, induced by the 
equivalence ~. 


IV. HOW DOES SYMMETRY HELP? 


Symmetry in the structure of a connecting network, together with the 
theoretically convenient (but in practice false) assumption of homoge- 
neous or interchangeable traffic sources, leads to simplifications in 
calculating state probabilities, loss, carried load, and other traffic 
quantities. In most cases the simplification occurs because network states 
that have the same structure to within renaming of customers, links, and 
switches are in a definite sense equivalent, and because if traffic sources 
are interchangeable, such equivalent states can turn out to have the same 
equilibrium or transition probabilities. Whether they do or not depends 
on the rest of the traffic model, especially on the rule used to make 
routing decisions: roughly speaking, if the rule is consistent in that it 
opts for analogous routes for analogous calls in equivalent states, then 
equivalent states will (or at least can) have the same probabilities. In such 
cases the state probabilities can be calculated from those of the equiv- 
alence classes, of which there are usually many fewer, by considering, 
in place of the original microscopic stochastic process on the set of states, 
a macroscopic one taking values on the equivalence classes. 

Our problem in this paper is to make all these notions, especially that 
of “equivalent” states, as precise as possible in the general network 
setting. In view of the central role of symmetry it is natural to expect that 
the equivalence idea be expressed mathematically by means of group 
theory, and particularly, in terms of the symmetry group of the graph 
G depicting network structure. Applications to optimal routing in net- 
works and concentrators will appear in later work.®:6 These applications 
are considerably complicated by the following “problem of refusals”: 
It turns out that analytical methods for finding optimal routing rules 
are greatly simplified if, as operator of the network, the telephone 
company is allowed the option of refusing to complete an unblocked call 
if it thinks that this denial of service will improve performance according 
to some criterion of interest; with this added option the task of finding 
out when to decline unblocked calls is part of the routing problem, a part 
which it turns out is usually much harder than actually choosing the best 
route if the call is to go in; however, it is often possible to solve the routing 
problem up to refusals, that is, to specify optimal routes for calls that 
might go in without ruling on whether they go in or not. 


V. PRIOR RESULTS OF S. P. LLOYD 


The relevance of group theory was well understood in 1955 by S. P. 
Lloyd, whose unpublished work! is now sketched. He used a standard 
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method of identifying an equivalence relation on a set with a group of 
bijections of the set into itself. This method considers a group G, of 
correspondences of S onto itself, and describes “equivalent” states as 
follows: G, is said to be transitive on a subset X of S iff 


(i) xeX, geG, = BxeX 
(tt) x,yeX = JgeG, such that gx = y 


X is then called a transitive set. With |X| the cardinality of X, it can 
be shown that |X| divides the order of G,, and that each member of X 
appears exactly 


IG] 
|X| 


times in the array gx, geG,, where x is any element of X. For each xeS 
define rx = {gx:geG,}. It can be seen that each zx is a transitive set, and 
that for any x,yeS we have either rx = wy or rx () ry = ¢. Thus the 7x 
form a partition of S, and so G, induces a corresponding equivalence 
relation = according to x = y iff there is a geG, such that gx = y. Con- 
versely, given an equivalence relation = on S, the set of all bijections g:S 
< S with gx = x forms a group G, under composition which induces = 
in the sense above, and we have the following “summation formula:” 
With a an equivalence class of =, x any member of a, and f areal valued 
function on S, 





|| 
fy) = ~ f(gx) 
yea | G,| geGy i 
Now let x; be a continuous-parameter Markov process taking values 
on S, with a stationary transition rate matrix Q = (q,,), assumed to be 
ergodic. There is then a unique probability vector p = {p,,x¢S} such that 
Pp solves the “statistical equilibrium” equations 


2, Pxdzy = 0, yes (1) 


p is the stationary probability distribution for x;. The group G, and the 
relation = become relevant to eq. (1) when the matrix Q of rates is un- 
affected by the permutations (of its rows and columns) corresponding 
to geG,; indeed, if G, and = express what we mean by saying that 
equivalent states differ only in respect of renaming customers, switches, 
etc., this is the precise way that symmetry affects the problem of calcu- 
lating state probabilities. This relevance is recognized in the fol- 
lowing 


Def. 1: Q admits G, iff geG,, x,yeS > 


Ixy = V(gx)(gy) 
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which leads quickly to Lloyd’s basic 1955 results: 


Theorem 1: Let G, bea group of bijections of S onto itself such that Q 
admits G,, and let E = {zx,x eS} be the set of equivalence classes induced 
by G,. Then 


(t) The projection map 7:x — ax defines a “macroscopic” Markov 


process 7x; on £, with transition rates 


dag = X& Qxy for xea and ack 
yeB 


and stationary probabilities {p,,ack} satisfying 


> Podag = 0, Bek 


(ti) For each xeaeE, 
1 
Px =77Pa 
|| 


Thus if Q admits G,, then equivalent states have the same stationary 
probabilities, and these can be computed from a reduced state equation 
of lower order. 


Proof of Theorem 1: Everything follows from the fact that if Bek, 
then 


»~ Qxy 
yep 


has the same value for all xcacE; so we prove this first. If zea then by 
transitivity of a there is a geG, with x = gz, so that 


DY Qxy = © Weezy = DX Vgeey) = DL Izy 
yeB yeB yeB yeB 


The second equality arises from g~! = @, the third from Q’s admitting 
G,. Since x and z were arbitrary elements of a, the result is proved. It 
implies, by results’ of M. Rosenblatt and C. J. Burke, that wx; is a 
Markov process with transition rates 


dap = X Qxy, any xea 
yep 


and 7x; is ergodic if x, is. Hence it has a unique stationary probability 
distribution {p,,a¢ E} satisfying 


a Pad@agp = 0, Bek 
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Thus (i) is proved; to prove (ii) we show first that 


Pa= X Px 


xe 


where {p,,x«S} is the stationary distribution of x;. We find 


aD PxQap = 2 Dx z Qxy = 3 > Px@xy = 0 


ack xea xeS yeB xeS 


since the inner sum in the last term is always zero. Thus 


be Dy, ae| 


XEN 


satisfy the (reduced) equilibrium equations for 7x;, and so by the uni- 
queness of its (probability vector) solution, 


Pa= X Dx 


xXEX 


Now define q by q, = |zx|~!p,x and consider that 


0= > Pedap = x [ae LD ed Gace 


= 2. Gx LD Qxy = Le 2 QxQxy 

xeS yeB yeBxeS 
However if z«@, there is geG, with y = gz, and (since q is constant on 
equivalence classes) qx = q(g-1x) so that 


x QxQxy = x QxQx (gz) = x  (g-1x)4 (g-Ix)z = 2, IxOx2 
XE 


Thus 


» Ox Ixy 
xeS 


is constant over equivalence classes and so it must be zero. Hence q is 
a probability vector solution of the equilibrium equation for x;, so gq = 
D, i.e., Dy = | xx|—!p ax. 

Lloyd’s theorem accurately captured the relevance of the rate matrix 
Q’s admitting the group G,, and he gave examples of the application of 
his result to small networks, such as individual switches and partial ac- 
cess concentrators, but he did not elaborate on the groups G, to be 
considered. However, in applying such a result to traffic in connecting 
networks we want to be sure that the groups G, we use reflect the intu- 
itive notions of invariance of structure under renaming of terminals, 
switches, etc. The theorem thus leaves us with these important questions: 
What groups G, are appropriate or available for describing equivalence 
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relations = useful in applications to traffic in networks? What traffic 
models give rise to rate matrices Q that admit these groups? 

Furthermore, since we expect the applications we make of Lloyd’s 
theorem to networks to depend on both network symmetry and customer 
interchangeability, not to mention routing, it would be well to have a 
formulation in which these items are clearly separated, as they are not 
in Theorem 1. What we need is a more specific definition of “equiva- 
lence” of states, one independent of probabilistic models, and peculiar 
to the network applications we intend, and one that reflects the idea of 
invariance of structure under renaming. These requirements will be met 
by constructing some groups that are appropriate from the symmetry 
group of the network graph G; Lloyd entertained® such an idea but did 
not describe it in Ref. 1. We shall first argue that certain natural neces- 
sary conditions, to be met by groups considered “appropriate,” imply 
that they should be automorphism groups of (S, <) whose elements 
preserve ~; then we show how such groups arise directly from the sym- 
metry group of the network graph. 


VI. TWO INTUITIVELY NECESSARY CONDITIONS 


We need properties and concepts that help make precise the notion 
of “structurally equivalent” states. Some of these will now be arrived 
at quickly and intuitively. Consider therefore two states x and y that 
differ only in point of renaming terminals and links, but otherwise have 
the same structure. In such a situation we expect to be able to make a 
correspondence between the calls in progress in x and those in y, because 
structural equivalence requires that each call c in progress in x have at 
least one analog in progress in y, playing the role of c in the structure of 
y. This being so, we see in the same way that the elements of B,, (ob- 
tainable from x by a hangup) have a natural correspondence to those of 
B,, going beyond the fact that |B,| = |B,| = |x] = ||: namely, to zeB, 
we assign a state in By obtained by hanging up an analog of the call hung 
up in x to get z, i.e., an analog of y(x — z). 

An exactly similar situation holds for the sets A, and A, of states 
which are respectively above x and y; actually, more is true, since for 
every call c free and not blocked in x there will be an analog (possibly 
more than 1), call it c’, free and not blocked in y, and to every way of 
putting up c in x, i.e., for every z¢€A¢x, there will be a corresponding way 
of putting up an analog c’ in y, and thus a natural correspondence of A, 
with A¢y. It is intuitively clear that A,’ cannot have either more or fewer 
states than A,,, else there would have to be something structurally dif- 
ferent about x and y to account for the discrepancy. 

Let now g denote the partial natural correspondence we have so far, 
with the properties §xX = y, &(By) a Box, and 8(A,(2-x)x) = Ay(gz-gx)gx 
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for zeA,. As we indicated, not even this much of g is unique. Nevertheless 
we suggest that g can be extended in a similar way, and to be defined on 
all of S in such a way as to satisfy the commutation g(B,) = Bg, and the 
distribution g(A,(2—-x)x) = Ay(gz—gx)gx for zeAy. It is not hard to see that 
these conditions are the same as 


xSy—>gx S gy 
x~ y = gx ~ BY 


The map g then takes any state into one equivalent in structure, in such 
a way as to preserve order in (S, <) and also equivalence in the other 
sense of having the same people talking. We can expect that g(S) = S, 
and that g is one-to-one. 

Now an isomorphism between two partially ordered systems (POS) 
is precisely a bijection that preserves order both ways; in our case the 
two POS coincide, so g is called an automorphism of (S, <). Since 
equivalence is transitive, it follows that if there are automorphisms g, 
h such that y = gx, z = hy, then there must be one f such that z = fx, 
namely the composition f = hg. Hence one is led naturally to consider 
groups of automorphisms that preserve ~. Thus while any equivalence 
relation on S can be described by a group of bijections on S, we claim that 
to adequately express what is meant by structure invariant under re- 
naming, the groups of interest for a theory of networks should be auto-: 
morphism groups of (S, <). Thus the first part of this study of equiva- 
lence = between states is the search, for a general network », for a group 
G,, of automorphisms of (S, <) that preserve ~ such that the “usual” 
definition of equivalence via a group, viz., 


x=yiff igeG, 5x =gy 


agrees with what we mean by structural equivalence. 


Vil. SYMMETRIES OF THE NETWORK GRAPH 


Structural equivalence of states rests ultimately on properties of the 
network graph G that are independent of and prior to the choice of inlets 
I, outlets Q, and states S that complete the description of a network vp 
= (G,J,0,S). In an informal way, one might say that the equivalence of 
two states under renaming of terminals and links really depends on what 
it means for the network to “look the same” to distinct terminals. As an 
example consider two arbitrary distinct inputs on the left side of the 
standard No. 5 crossbar type network in Fig. 1. It is obvious intuitively 
that if the frames and switches are identical, and the connections within 
and between frames correspond to complete bipartite graphs in the usual 
way, then the network “looks the same” to two such inlets. The same is 
true of any two interframe junctors, or of any two links from the same 
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Fig. 1—Connecting network. 


or from two distinct frames. We seek to clarify this informal notion of 
“looking the same,” and to develop it into a precise definition of struc- 
tural equivalence in terms of a natural symmetry group for the network 
graph. 

Let us think of the terminals of v as the vertices, and of the crosspoints 
as the edges, of the network graph G. It is clear that G is determined by 
giving a relation N on the set 7’ of terminals such that for ¢,s in T, tNs 
iff there is a crosspoint or edge between ¢t and s. N is the symmetric 
“‘nextness” or adjacency relation that completely depicts the network 
structure. 

Now suppose that we rename the terminals in T according to some 
permutation 7. Most permutations would play havoc with the adjacency 
relation; that is, if t and s had a crosspoint between them, then 7t and 
ts easily might not, and conversely. But there might be some permuta- 
tions other than the identity which preserved N in the sense that for 
every t,seT 


tNs iff (rt)N(zs) 


In this case the permuted terminals have crosspoints between them in 
exactly the same pattern as the unpermuted. It is the existence of such 
an “N-preserving” permutation 7 that we take as the precise meaning 
of “looking the same.” 

The network “looks the same” to two terminals t and s iff there is an 
N-preserving permutation 7 of T into itself such that s = zt. It is ap- 
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parent that these N-preserving permutations form a group, which we 
call the symmetry group of the network graph G. We remark that this 
group may be trivial (if there are no N-preserving permutations except 
the identity), and that in any case it in no way depends on what terminals 
have been designated as inlets or outlets, or on what ways of closing the 
crosspoints are to be allowed as physically meaningful states. 

In most telephone connecting networks the set J of inlets and that Q 
of outlets are fixed sets of terminals, and one is not interested in whether 
the network “looks the same” to an inlet as it does to an intermediate 
link or junctor; in most cases it will not, in any case, because of their 
different functions in operation. So it makes sense to restrict the N- 
preserving bijections we think of as renamings of switches and terminals 
to those which either preserve both J and Q, or else map each onto the 
other. This restriction defines a subgroup Gy, called the symmetry group 
of G for J and Q, which will be used to define structural equivalence of 
states. 


Vill. SYMMETRIES OF THE SET OF STATES 


The set S of states of a network pv = (G,J,Q,S) represents all the ways 
of closing the crosspoints which we regard as physically sensible. It is 
closed under hangups, that is, under removal of a chain from a state; it 
need not be closed under adding new chains, nor even under adding new 
chains which by themselves already represent a state with one call in 
progress. It is convenient, however, to require that sets of chains in a 
structural equivalence class either all belong to S, or that none of them 
does. This requirement of course implicitly assumes that we know what 
structural equivalence is before we choose states for S; it will be seen that 
the definition of equivalence below applies to arbitrary sets of chains 
on G, so the requirement can be met as we choose such sets to belong to 
S. Specifying the set S represents definite choice of the ways in which 
the network with graph G is to be used. 

We shall now use the symmetry group G,, of the network to define 
what we mean by two states’ differing only in point of renaming links, 
terminals, etc. Indeed, it can be seen intuitively that the symmetry group 
G,, of the network induces a natural equivalence relation on sets of chains 
on G, and thus on whatever such sets we choose as states: two sets of 
chains on G are equivalent if there is some group element 7¢G,, such that 
a terminal t is busy in the first iff rt is busy in the second. The incorpo- 
ration of a map 7¢G,, in the definition ensures that simultaneously the 
network looks the same to a terminal ¢ and to its analog rt. For states, 
then, we define structural equivalence = by 


Def. 2: x =y iff JreG, 3 t1,...,tm are the terminals busy in x iff 7, 
...,Ttm are those busy in y. 
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It is apparent that since G,, is a group, = is an equivalence relation on 
S, partitioning S into a set E of equivalence classes a,3,..., and inducing 
a natural projection map 7: S — E such that 


rx = {yy =x}eh 


The elements a of EF will be the reduced network states. As in Section 
V, we see that there is an associated group G, that provides an alternative 
description of =; G, is the group of all bijections g: S <> S which map each 
equivalence class a into itself. In fact, G, is the largest strictly imprimitive 
group of bijections of S onto itself whose sets of imprimitivity are pre- 
cisely the ack: 


Def. 3: G, = {g:S << S 3 g(a) = a for ack} 
Remark 1: As in Section V, we have x = y iff JgeG, 3 gx = y. 


Remark 2: Although G, is a group of bijections which does characterize 
=, it is typically not economical. It turns out that a much smaller 
subgroup of G,, defined directly from G,,, suffices to characterize = as 
in Remark 1. These subgroups appear in Section 10, and they are the 
~-preserving automorphism groups desired in Section VI. 


IX. ACTIONS OF 7eG, ON A AND S 


An assignment a is a correspondence or injection from a subset of I 
into 2. Thus an element 7eG,, acts in a natural way on an assignment acA 
to produce a new assignment ra according to the rule that u and v cor- 
respond in a iff ru and 7v correspond in 7a. Similarly an element 7¢G,, 
acts in a natural way on a set X of chains on the graph G to produce a 
new set 7X of chains consisting of the r-images of chains in X, thus: ¢,, 
...,t)1s to be a chain of X iff 7t1,...,7t, is a chain of 7X. In particular 
a teG,, acts on a state xeS to produce a set of chains 7x, and it is rea- 
sonable to assume, as we have done here, that the choice of S is consistent 
with the symmetry group G,, in that S is closed under the action of any 
TeGnixeS, T6EG, = 7Tx6S. This will ensure that either all the sets of chains 
in a structural equivalence class are states, or none of them is. 


~ Remark 3: Since 7 may map 2 onto J its action on acA may reverse the 
“usual” order of the pairs (u,v)el X 2 to (ru,7v)eQ X I. Therefore we do 
not distinguish between an assignment 


a = {(u,v)el X Q:(u,v)ea} 
from its inverse 
a~! = {(v,u)eQ X I:(u,v) ea} 


or else we specify that when 7J = Q then the action is defined by 


ta = {(rv,Tu):(u,v) ea} 


REDUCTION OF NETWORK STATES 123 


X. HOMOMORPHISM OF G, INTO G, 


The network symmetry group G,, depends only on the basic network 
structure: the adjacency relation N of the graph G that depicts the 
network, and the choice of J and Q, since reG,, are restricted to preserve 
inlets and outlets. The state symmetry group G, describing =, however, 
also depends on what sets of chains on the graph G are chosen as phys- | 
ically interesting or important states. Since G,, is used in the definition 
of =, and thus of G,, it is not surprising that G, and G, happen to be al- 
gebraically related: there is natural homomorphic image G, of G, in G,, 
consisting entirely of ~-preserving automorphisms of (S, <), and de- 
scribing the same equivalence relation =. Thus for most purposes it is 
more convenient to use the “equivalent” subgroup G, than the full 
symmetry group G, associated with = by the standard method. The big 
group G, induced by =, incidentally, is not necessarily an automorphism 
group: consider, e.g., the 2 X 2 switch of Fig. 4 supra; clearly 5 =6 and 1 
= 8, but 1 < 5 and 3 < 6. Thus the map g defined by the permutation (13) 
(24) (56) belongs to G, but is not an automorphism. 


Theorem 2: The action of reG,, on states x defines a homomorphism 
n:G, — G, according to the rule that n(r)x = rx; the image group G, = 
n(G,,) is an automorphism group that preserves ~. 


Proof: For reGy, define n(r)eG, by the condition n(r)x = rx that ¢j,... 
yt; is a chain of xeS iff 7t1,... ,7t) is a chain of n(7)x(=7x). The subset 
n(G,,) is closed under composition, so it is a subgroup G, of G,. Ifx < 
y, then n(7)x S n(7)y for each reG,, so that each n(7) is an automorphism 
of (S, <). It is easily verified that n(7172) = n(71)n(72); preservation of 
~ follows from the definition of action. 


Theorem 3: If every terminal teT is busy in some state xeS, then the 
homomorphism 7 of Theorem 2 is an isomorphism, i.e., it is injective 
(one-to-one). 


Proof: If (71) = n(r2), let x be a state with one call up in which a terminal 
t is busy. Then x consists of a single chain ¢;,...,t; such that t is some 
t; and r,t; = rot; for 1 <j <1. Then 7yt = 79t, and since ¢ is arbitrary we 
have 71 = T9, and so 7 is injective. 


Remark 4: The following conditions are all equivalent: 


(i) x=y 

(ii) For some reGpn, x = TY = n(T)y 
(iii) For some heG,, x = hy 
(iv) For some geG,, x = gy 
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Fig. 2—Small 3-stage network. 


XI. PARTIAL ORDERING OF REDUCED STATES 

It is natural to try to partially order the set E of reduced states ac- 
cording to this idea: an equivalence class a is “above” another 6 if some 
element xea can be reached from some member ye@ by adding more calls, 
i.e., if x = y for some xea and ye’. Formally, we make the 


Def. 4: a2 Biff Jxea,yBax zy 

For this to define a partial ordering it must be reflexive, antisymmetric, 
and transitive. The first is obvious; we prove the other two. Let then a 
> 6B and 6 = a; to show a = G it is enough to show a () 6B ¥ ¢, because 
each is an equivalence class. There exist x,zea and y,we@ such that x = 
y and w 2 z. Since x =z and y = w, there exist automorphisms 71,72G , 
such that 71x = z and roy = w. Hence 71x = 71y and 


Toy =w2z2=7X% 2 71y 


so that 7;~!roy = y. But astate z that is above another 71y and has the 
same number of calls in progress equals it: z = 71y.Soa()\ 6 # @. 

To prove transitivity let a = 6 and 6 2 y, so that there exist xea, ye8, 
ze3, and wey such that x = y, y =z, and z = w. There is an automorphism 
7G, such that ry = z, whence 


TX 2TY=Z22wW 
Hence there is something in a which is above something in yj, i.e., a 2 
y. Thus 2 is a partial ordering on E = 7(S). 
Remark 5: x 2y = 4x = wy 


Remark 6: The partial order = on E need not be a semilattice, and a 
fortiori the projection z:S — E need not be a semilattice homomorphism, 
as is y:S — A, if A is ordered by inclusion. Figures 2 and 3 provide a 
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Fig. 3—Reduced state diagram for small 3-stage network of Fig. 2. 


counterexample, as well as an illustration of the reduced states of a 
network. 


Xll. 2 X 2 SWITCH 


The next four sections are devoted to increasingly complex examples 
that will illustrate the notions we have introduced to make precise the 
idea of structural equivalence. Our first and simplest example is the 2 
X 2 switch shown in Fig. 4. We arbitrarily label the inlets 1 and 2 and the 
outlets 3 and 4. So with J = {1,2} and Q = {3,4} the network graph G isa 
square with the terminals of J on diagonally opposite vertices, and 
similarly for 0. The adjacency relation N consists of exactly the pairs 
(1,3), (1,4), (2,3), and (2,4) together with the results of interchanging the 
first and second members of these pairs so as to make N symmetric. 

The maps of the terminal set T = {1,2,3,4} into itself which preserve 
N and either preserve both J and Q, or carry one onto the other, are ex- 
actly the permutations 

identity, (12), (34), (12)(34), (18)(24), (23)(14), (1428), (1324) 
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TWO-—BY—TWO SWITCH EXAMPLE 








4 2 


Fig. 4—Network graph, adjacency matrix, and states for 2 X 2 switch. 


These 4-permutations form the symmetry group G,, of the 2 X 2 switch 
for J = {1,2}, Q = {3,4}. This group has a nice geometric meaning in terms 
of the network graph, the square in Fig. 4: it consists precisely of all the 
rotations and reflections of the square into itself. The multiplication 
table for this group is given in Table I, along with “generators” A = 
(1423) and B = (12) which, under the relations A4 = B? = identity, BA 
= A3B, identify the group as (an isomorph of) the dihedral group D4, of 
order 8. Every terminal is busy in some state, so Theorem 3 applies and 
we need only calculate G,, to define =, instead of passing through G,,. 
The symmetry group G,, describes the equivalence of states through 
the action table in Fig. 5. Here the columns are indexed by states xeS and 
the rows by group elements 7¢<G,,, and the 7,x entry is the state 7x defined 
by the action of 7 on x, i.e., by replacing every chain (t1,f2)ex by (7t1,7¢9). 
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Table | — Multiplication table of G, for 2 X 2 switch 


f B A2B A2 AB A3B A A3 
g I (12) (34)  (12)(34) (13)(24) (14)(23) (1423) (1324) 
I I (12) (34) (12)(84) (13)(24) (14)(23) (1423) (1324) 
(12) (12) I (12)(34) (34) (1324) (1423) (14)(23) (18)(24) 
(34) (34)  (12)(34) I (12) (1423) (1324) (13)(24) (14)(23) 
(12)(34) (12)(34) (34) (12) I (14)(23) (13)(24) (1324) (1423) 
(13)(24) (18)(24) (1423) (1824) (14)(23) I (12)(34) (12) (34) 
(14)(23) (14)(23) (1824) (1423) (13)(24) (12)(34) I (34) (12) 
(1423) (1428) (18)(24) (14)(23) (1324) (34) (12)  (12)(34) I 
(1324) (13824) (14)(23) (13)(24) (1423) (12) (34) I (12)(34) 





I 



















(12) (12) (34) (56) 


(34) (14) (23) (56) 
(12) (34) (13) (24) 
(13) (24) (24) 

(14) (23) (13) 
(1423) (1234) (56) 





(1324) (1432) (56) 






Tx(=(T)x) 






Fig. 5—Action and isomorphism for the 2 X 2 switch. 


Each row of the table defines a bijection of S = {0,.. . ,6} into itself that 
corresponds to the row index 7eG,,, and thus an isomorphism 7:G,, — 
S7 depicted by 
rt |id (12) (34) (12)(384) (18)(24) (23) (14) 
n(z) lid (12)(34)(56) (14)(23)(56) (13)(24) (24) — (18) 
(1423) (1324) 
(1234) (56) (1432) (56) 


It can be verified that G,, = n(G,) is an automorphism group of (S, S), 
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INLETS / (LINES) 


1 2 3 4 5 6 


OUTLETS 
Q 


(TRUNKS) 


10 


Fig. 6—Random slip concentrator. 


isomorphic to D4, and preserving ~. It remains to check that = defined 
with the help of G,, is in fact the intuitive one that is indicated in Fig. 
4: looking at the values of 7 in the two line table above we easily see that 
5 =6,1=2=83 =4, and no other distinct states are equivalent, so the 
equivalence classes are {0}, {1,2,3,4}, and {5,6}, as we had guessed. The 
group G, of all bijections which preserve these equivalence classes is 
isomorphic to S_ X S4, a group of order 48; in contrast, the subgroup G, 
is sufficient to characterize = and is only of order 8. 


XII. RANDOM SLIP CONCENTRATOR 


The crosspoint assignment for a 6-line to 4-trunk random slip con- 
centrator is shown in Fig. 6. The phrase ‘“‘random slip” is old telephone 
terminology that clearly originated as a description of the even or regular 
way in which the incomplete access of lines to trunks is distributed over 
the network, something like the statisticians’ balanced block designs. 
Figure 6 may also be depicted as the labelled tetrahedron of Fig. 7, with 
the interpretation that an edge (line) ‘‘has access” to the two vertices 
(trunks) that it connects. Figure 7 leads naturally to the network graph 
(Fig. 8) and the adjacency matrix (Table II): just add a vertex of degree 
2 in the middle of each edge, with the same label as the edge. In this 
network every t is busy in some state, so Theorem 3 applies, G, ~ Gp, 
and we need only calculate G,. The reduced states are shown in a con- 
venient representation as a partially ordered system in Fig. 9. 

Thus we seek maps of J = {1,2,.. . ,6} and 2 = {7,8,9,10} into themselves 
which preserve adjacency. It can be seen from Figs. 6 and 7 that every 
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10 


Fig. 7—Labeled tetrahedron. 


5 


Fig. 8—Network graph for concentrator. 


Table Il — Adjacency matrix for concentrator 

2 3 4 5 6 7 8 9 

1 1 0 

1 0 1 

0 1 0 0 

0 1 1 

0 1 0 

0 0 1 
1 1 0 0 0 

0 0 1 1 0 0 
1 0 1 0 1 
0 1 0 1 1 


Re OF OO 


EDGES =LINES 
VERTICES =TRUNKS 
—e—— ARROW POINTS 
TO TRUNK IN USE 


Fig. 9—Reduced states of random slip concentrator. 


permutation of {7,8,9,10} forces a unique permutation of {1,.. . ,6} if ad- 
jacency is to be preserved; graphically, we see that all the maps we seek 
are rotations or reflections of the tetrahedron of Fig. 7 into itself, and 
that all permutations of {7,8,9,10} are allowed. More formally, every 
transposition of Q is allowed and requires a unique pair of disjoint 
transpositions in J according to the table 


(78) forces (24)(35) 
(89) forces (12)(56) 
(9 10) forces (23)(45) 
(7 10) forces (15)(26) 
(7 9) forces (14)(36) 
(8 10) forces (13)(46) 


Since the allowed maps form a group, and every permutation of Q is a 
product of transpositions, every permutation of 2 forces a unique one 
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Table Ill — G, for concentrator 


id. id. 
(78) (24) (35) 
(89) (12)(56) 

(9 10) (23) (45) 
(7 10) (15)(26) 

(79) (14) (36) 

8 10) (13)(46) 
(78)(9 10) (25)(34) 
(89)(7 10) (16) (25) 
(79)(8 10) (16) (34) 

(789) (142)(356) 

(798) (124) (356) 

(78 10) (153)(246) 
(7 10 8) (135)(264) 
(79 10) (145)(263) 
(7 10 9) (154) (236) 
(89 10) (123)(465) 
(8 10 9) (132)(456) 
(789 10) (1463) (25) 
(798 10) (2453)(16) 
(79 10 8) (1265) (34) 
(78 10 9) (1562)(34) 
(7 10 89) (3542) (16) 
(7 10 98) (3641)(25) 


of J, as claimed. It follows that G, is isomorphic to S4, and consists of 
the permutations shown in Table III. The ‘“‘forced” permutations of I 
are of course a subgroup of Sg isomorphic to S4. 


XIV. CLOS NETWORK EXAMPLE 


The crosspoint structure and network graph for the simplest 3-stage 
Clos rearrangeable network are shown in Figs. 10-11. A planar form of 
the network graph is on Fig. 12. At the start of this paper we loosely de- 
scribed the invariances of structure to be studied as those associated with 
“renaming terminals, switches, and links.” The little Clos network now 
under discussion gives us a specific example of what this means: it means, 
of course, permuting these entities while preserving the structure, e.g., 
interchanging switches in a stage while dragging along each switch’s rat’s 


~-—-72 X 2 SWITCHES ~~ _ 
1 ~™ 





Fig. 10—Simple 3-stage Clos network. 
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Fig. 11—Crosspoint structure and network graph for 3-stage Clos network made 
of 2 X 2 switches. 


nest of links that connect to other stages. In particular, the organization 
of crosspoints into switches must also be preserved by these permuta- 
tions of terminals, switches, and links. In the next section (XV) we shall 
describe how this constraint provides an approach to calculating the 
group G,, for general crossbar networks. For the present we show how 
the approach applies to the example under discussion. 

Any 16-permutation that preserves the adjacency or crosspoint 
structure depicted in Fig. 11 must permute the inner terminals or links 
(5,6, ... ,12} among themselves only, and it can do so in exactly the 16 
ways shown in Table IV. It is readily seen that (i) these 16 ways corre- 
spond to interchanging switches, and rotating or reflecting the network, 





Fig. 12—Planar form of network graph of Fig. 11. 
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Table IV — Allowed permutations of inner terminals for 3-stage 
Clos network of 2 X 2 switches 


No. Map Type Action, if simple 

1 id. Ss move no switches 
2 (57)(68)(9 11)(10 12) ss interchange m. sw. 
3 (9 10)(11 12) sm interchange rt. sw. 
4 (57)(68)(9 12)(10 11) sm interchange rt. and m. sw. 
5 (56)(78) ms interchange | sw. 
6 (58)(67)(9 11)(10 12) ms interchange ]. and m. sw. 
7 (56)(78)(9 10)(11 12) mm interchange rt. and 1. sw. 
8 (58)(67)(9 12)(10 11) mm rt., m., and I. 
9 (59)(7 11)(6 10)(8 12) aa rotate about vertical axis 

10 (5 11)(79)(6 12){8 10) aa 

11 (5 10 69)(7 12 8 11) ad 

12 (5 126 11)(7 10 89) ad 

13: (596 10)(7 11 8 12) da 

14 (5 11 6 12)(7 9 8 10) da 

15 (5 12)(6 11)(7 10)(89) dd 

16 (5 10)(69)(7 12)(8 11) dd 


and (ii) that they form a subgroup of Sg. Each of these ways can “go 
with” a number (here always 16) of permutations of outer terminals 
among themselves, to form an element of G,. By writing down the 
possible ways this matching can be done (so as to preserve adjacency) 
we get a brute force way of calculating the group G,, for the 3-stage Clos 
network made of 2 X 2 switches. 

The matching in question has a block structure: the permitted per- 
mutations of links can be partitioned in such a way that to each parititon 
element there corresponds a set of permitted permutations of outer 
terminals any one of which can “go with” each link permutation in the 
element. This block structure is indicated in Table IV by the type 
symbol; all allowed link permutations of the same type can match with 
all the same permutations of the outer terminals. Thus to present G, 


Table V — G, for Clos network example 


Type ss sm 

id. (13 1)(14 16) 

(12) (13 16)(14 15) 

(34) (12)(13 15)(14 16) 
(13 14) (12)(13 16)(14 15) 
(15 16) (34)(13 15)(14 16) 
12)(34) (34)(13 16)(14 15) 

(12)(13 14) (12)(34)(13 15)(14 16) 


(12)(15 16) 

(34)(13 14) 

(34)(15 16) 
(13 14)(15 16) 


(12)(34)(13 14) 
(12)(34)(15 16) 
(12)(13 14)(15 16) 
(34)(13 14)(15 16) 
(12)(34)(13 14)(15 16) 


(12)(34)(13 16)(14 15) 

(13 15 14 16) 

(13 16 14 15) 
(12)(13 15 14 16) 
(12)(13 16 14 15) 
(34)(13 15 14 16) 
(34)(13 16 14 15) 

(12)(34)(13 15 14 16) 
(12)(34)(13 16 14 15) 
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Type 


Type 


Type 


Table V — (Continued) 


ms 


(13)(24) 

(14) (23) 
(13)(24)(13 14) 
(14)(23)(13 14) 
(13)(24)(15 16) 
(14)(23)(15 16) 

(13)(24)(13 14)(15 16) 
(14)(23)(13 14)(15 16) 
(1324) 
(1324)(13 14) 
(1423)(13 14) 
(1324)(15 16) 
(1423)(15 16) 
(1324)(13 14)(15 16) 
(1423)(13 14)(15 16) 


aa 


1 13)(2 14)(3 15)(4 16) 
1 13)(2 14)(3 16)(4 15) 
1 14)(2 13)(3 15)(4 16) 
1 14)(2 13)(3 16)(4 15) 
13)(2 14)(3 15 4 16) 
13)(2 14)(3 16 4 15) 
14)(2 13)(3 15 4 16) 
4)(2 13)(3 16 4 15) 
2 14)(3 15)(4 16) 
14)(3 16)(4 15) 
13)(3 15)(4 16) 
13)(3 16)(4 15) 
4)(3 15 4 16) 


1 
1 
13)(3 15 4 16) 
13)(3 16 4 15) 


—se ea 
mS be me 09 09 
Dmvypnnrndnd 


fli 
Re 
mm mm CO OO 


(1 13 4 15)(2 14 3 16) 
(1 13 4 16)(2 143 15) 
(1 13 3 15)(2 14 4 16) 
(1 13 3 16)(2 14 4 15) 
(1 14 4 15)(2 13 3 16) 
(1 14 4 16)(2 13 3 15) 
(1 14 3 15)(2 13 4 16) 
(1 14 3 16)(2 13 4 15) 
(113315214 4 16) 
(113316 2144 15) 
(113 4152143 16) 
(113 416 2 143 15) 
(1143152134 16) 
(1143 16 213 4 15) 
(1144152133 16) 
(1144 16 2133 15) 


mm 


(1324)(13 15)(14 16) 
(1324)(13 16)(14 15) 
(1423)(13 15)(14 16) 
(1423)(13 16)(14 15) 
(13)(24)(13 15)(14 16) 
(13)(24)(13 16)(14 15) 
(14)(23)(13 15)(14 16) 
(14)(23)(13 16)(14 15) 
(1324)(13 15 14 16) 
(1423)(13 15 14 16) 
(1423)(13 16 14 15) 


(13)(24)(13 15 14 16) | 


(13)(24)(13 16 14 15) 
(14)(23)(13 15 14 16) 
(14)(23)(13 16 14 15) 


ad 
(1 153 13)(2 16 4 14) 
(1 16 3 13)(2 15 4 14) 
(1 15 4 13)(2 16 3 14) 
(1 16 4 13)(2 15 3 14) 
(1 15 3 14)(2 16 4 13) 
(1 16 3 14)(2 15 4 18) 
(1 15 4 14)(2 16 3 13) 
(1 16 4 14)(2 15 3 18) 
(115313216414) 
(1163132154 14) 
(115413 2 163 14) 
(1164132153 14) 
(115314216 4 13) 
(1163 14215 4 13) 
(1154142163 13) 
(1164142153 13) 


(1 15)(2 16)(3 13)(4 14) 
(1 15)(2 16)(3 14)(4 13) 
(1 16)(2 15)(3 13)(4 14) 
(1 16)(2 15)(3 14)(4 13) 
1 15)(2 16)(3 13 4 14) 
15)({2 16)(3 14 4 13) 
6)(2 15)(3 13 4 14) 
6)(2 15)(3 14 4 13) 
2 16)(3 13)(4 14) 

16)(3 14)(4 13) 

15)(3 13)(4 14) 

15)(3 14)(4 13) 


( 
(1 
(1 
(1 
(1 
(1 152 

(1162 

(1162 

2 16)(3 13 4 14) 
2 

2 

2 


1 
1 
15 
15 
16 
16 
16)(3 14 4 13) 


15)(3 13 4 14) 


(1 
(1 
(1 
(1 15)(3 14 4 13) 


it is enough to list the “outer” permutations that match each type of 
“inner”; this is done in Table V. The reduced states are in Fig. 13. 
The multiplication table of the group of “inner” or switch permuta- 
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Fig. 13—Reduced states for 3-stage Clos network made of switches. 


tions can be calculated from the cycle representations of Table IV; it is 
displayed in Table VI, from which it can be seen that elements R, = 2 
Ro = 4, and R3 = 15 satisfy the relations 


(R1R3)? = (RoR1)* = (ReoR3)4 = identity 


which identify the group as (an isomorph of) S2 X D4. This factorization 
is related to the fact that the type of a product is determined by the types 
of the factors in a way that is summarized in Table VII, in which the row 
is the type of the first factor, and the column that of the second. Indeed 
the type symbols themselves form an isomorph of D4 under the “mul- 
tiplication” Table VII. What is more, we can identify this group as the 
subgroup (~ D4) of the link or switch permutation group which restricts 
attention to the motion of the outer switches. Clearly the motion of the 
two middle switches is independent of that of the outer ones; this fact 
shows up in the feature that every type consists of just two permutations, 
and in the factor S2 in the switch permutation group. The possible in- 
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Table VI — Multiplication table for switch permutation group 


id. 
1 2 3 4 5 6 7 8 9 10 Il 12 #18 «14 «215~«216 
id. 1 2 3 4 5 6 7 8 9 10 11 12 18 14 «215 «16 
2 1 4 3 6 5 8 7 10 9 12 11 14 #18 #16 «15 
3 4 1 2 7 8 5 6 11 12 9 10 16 #15 14 «18 
4 3 2 1 8 7 6 5 12 11 10 9 15 16 138 14 
5 6 7 8 1 2 3 4 13 14 #16 15 9 10 12 Ii 
6 56 8 7 2 1 4 8 #14 #18 #15 #16 10 9 #=Ii1 12 
7 8 5 6 3 4 1 2 16 15 13 14 I1 412 #10 «9 
8 7 6 5 4 3 2 1 15 16 14 13 12 11 9 10 
9 10 18 14 #I11 #12 16 15 tI 2 5 6 3 4 8 7 
10 9 14 #13 12 11 15 162 1 6 5 4 3 7 8 
11 12 16 15 9 10 18 #14 8 44 7 8 1 2 6 = 5 
12 11 15 16 10 9 14 13 4 3 8 aL. 2 1 5 6 
13 14 9 10 16 15 11 12 5 6 41 2 7 #8 4 8 
14 13 10 9 15 16 12 11 #6 5 2 1 8 7 3 4 
15 16 12 11 14 #18 10 9 8 7 44 8 6 5 I1 2 
16 15 11 #12 «213 «14 9 10 7 8 3 4 5 6 2 1 
Type ss ss sm sm ms ms mmmmaa aa ad ad da da dd dd 


terchanges of outer switches are summarized in Fig. 14 and they corre- 
spond to the type symbols as follows 


switch interchange 1 r l h Uv x a b 
type symbol 1 sm ms aa mm dd <da_ ad 


The switch interchange symbols of course also form a group ~ D4 under 
the multiplication transferred from the type symbols; indeed {1,h,v,x}, 
{1,r,l,v}, and {1,a,6,v} each forms a sub-” Vierergruppe”’: 


XV. CALCULATION OF G, FOR CROSSBAR NETWORKS 


In connecting networks built out of stages of rectangular crossbar 
switches, a permutation of the terminals cannot preserve adjacency in 
the network graph unless it maps the terminals of each switch onto those 
of some switch, either itself or another one, so as to either preserve or 


Table VII — ‘‘Multiplication’’ table for type symbols. Top row: 
generators identifying with D4. Left column: identification with 
outer switch interchanges. 


A3B AB A A3 A? B A2B 1 
aa dd da ad mm ms sm Ss 
h aa Ss mm sm ms dd ad da aa 
x dd mm SS ms sm aa da ad dd 
a da ms sm mm Ss ad dd aa da 
b ad sm ms Ss mm da aa dd ad 
D mm dd aa ad da ss sm ms mm 
l ms da ad aa dd sm Ss mm ms 
r sm ad da dd aa sm mm Ss sm 
1 Ss aa dd da ad mm ms sm Ss 
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Fig. 14—Possible interchanges of outer switches. 


interchange inlets and outlets. It follows that any such permutation 
determines a permutation of the switches, and that the set G, of such 
induced “switch permutations” is again a group, which we call the switch 
permutation group. The task of calculating G,, for a crossbar network 
is substantially simplified by first finding Gs. 

In the language of group theory, the symmetry group G,, for a crossbar 
network must be an imprimitive group, because it consists of permu- 
tations that either preserve switches or map them onto each other. In- 
deed it can be seen that the map ¢ which assigns to each geG,, the switch 
permutation that g induces is a homomorphism of G, onto G;. Once a 
permitted switch permutation peG, is chosen, each element in ¢~}{p} 
is determined by choosing, for each permuted outer switch, a map which 
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Fig. 15—-Frame with “canonical” complete bipartite cross-connect field. 


assigns the outer terminals of that switch to the outer terminals of its 
image under p. These maps are independent, so we have proved 


Theorem 4: If vis a crossbar network, then 


Gn ~ Gs X J] Sjo) 
0 


where the product is over outer switches 0, |o| is the number of inlets 
(outlets) on o and S; is the symmetric group on k objects. 


For many crossbar networks made of stages it is rather straightforward 
to calculate the switch permutation group G,, because G, turns out to 
be isomorphic to a semidirect product of groups that are determined by 
the way the stages are joined together. 


Example: Frame. It is clear that for a frame made of two identical 
stages, interconnected by the “canonical” (complete bipartite) cross- 
connect field as shown in Fig. 15, G, will consist of all maps that take the 
inlet switches onto the outlet switches and vice versa, together with all 
maps that permute the inlet switches among themselves, and also the 
outlet switches among themselves. If each stage has n switches, then G, 
is isomorphic to the largest imprimitive permutation group on 2n objects 


with two equinumerous sets of imprimitivity. If 1,...,n are the inlet and 
n+1,...,2n are the outlet switches respectively, then 

Gs =F, U oF, 
where F,, is the isomorph of (S,,)? in Se, which permutes 1,... ,n among 
themselves and alson + 1,...,2n among themselves and o = (1n + 1)(2 


n + 2)--(n 2n). Each switch is n X n, so 

Gr a (F, U oF ,) x (S,)2” 
It can be verified that the union F,, U oF, is in fact a group and asem- 
idirect product. The order of G, for the frame is 2(n!)?. For n = 2, G, ~ 
D4, the dihedral group of order 8. 
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Fig. 16—Three-stage Clos network. 


Example: 3-Stage Clos Network. From Fig. 16 it is evident that G, is 
the imprimitive group which permutes the middle switches, and either 
permutes the inlet switches and the outlet switches independently, or 
else maps inlet switches onto outlet and vice versa. If now K,, is the iso- 
morph of (S,.)? in Sy,.42n Which permutes 1,...,2 andm+n+1,...,m 
+ 2n independently, and L,, is the isomorph of S,,. in Sm+2n which 
permutes only n + 1,...,m +n, then for the Clos network of Fig. 16, 


G, = (Ky, U 7Kn) X Ln 


where 7 = (1 2n + 1)--(n 3n), and G, ~ G, X (S,,)2". For the case m = 
n = 2, treated in detail in Section XIV, we have G, ~ So X Dg and G, ~ 
(S2)° X D4. It is now easy to see that the S. X D4 structure of G, arises 
from viewing the outer stages as a frame that yields D4 as in the previous 
example, while the S2 arises from permuting the middle switches. 


Example: Cascade of s stages with “complete bipartite” cross-connect. 
The structure of this network is shown in Fig. 17. The parity of s plays 
a role here, as follows: If s is odd, then under a switch permutation the 
switches of the middle stage can only go into each other, and as we saw 
in the previous example, they contribute to G, a group factor isomorphic 


STAGE NO. 1 2 eee {s/2] +1 eve S$ 


n {s/2] +1 
a 


nxn SWITCHES n [s/2] +n 





Fig. 17—-Cascade of s stages with “complete bipartite” cross connect. 
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n2 +4 2n2 +1 3n2 +1 





Fig. 18—No. 5 crossbar type network. 


to S,. Whatever the parity of s the noncentral switches contribute a 
factor isomorphic to an imprimitive group (not the largest) with 2[s/2] 
sets of imprimitivity; these are partitioned into pairs and the group el- 
ements either map each set onto itself, or each set onto its paired set, and 
never some of each. If the switches are numbered as in Fig. 17, this group 
is essentially the semidirect product 


M, U 6M, 


where M,, is the largest subgroup of S,, which fixes all the central 
switches and permutes vertically within every other stage, and 


[s/2] 
6= JJ (kn+1(s-—1—k)n+1)-(kn +nsn — kn) 
k=0 


Finally 
Gs ~ (Mn U 0Mn)(S,)s~ 2821 


Example: No. 5 Crossbar. With the switches numbered as in Fig. 18, 
let P,, be the isomorph of (S,,)? in S4,2 which for k = 0,1,2,3, permutes 
kn? +1,...,(R + 1)n? among themselves in such a way as to map frames 
onto frames, outer stages onto outer of the image frame, and inner stages 
onto inner stages of the image frame, and inner stages onto inner stages 
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of the image frame. Let 
o = (1 38n2 + 1)(2 8n2 + 2)--(n? 4n2) 
X (n2 + 1 Qn2 + 1)(n?2 + 2 Qn? + 2)---(2n2 3n?) 
Then 
G, =P, U oP, 
Gn = (Pn U oPp) X (Sp)? 


The remainder of this paper derives the reduced state equations for 
the stochastic model of Ref. 2. 


XVI. ROUTING OF CALLS 


We shall use a routing matrix R = (rx) as a convenient formal de- 
scription of how routes are chosen for calls. The class of routing matrices 
R can be described thus: for each x«S let I, be the partition of A, in- 
duced by the relation ~ of “having the same calls up,” or satisfying the 
same assignment of inlets to outlets; it can be seen that II, consists of 
exactly the sets A., for c free and not blocked in x; for YeTI1,, r, for yeY 
is to be a probability distribution over Y, that isr,y 2 Oand DiyeyPrxy = 
1; ry is to be 0 in all other cases. 

The interpretation of the routing matrix as a method of choice is to 
be this: any YeIL, represents all the ways in which a particular call c (free 
and not blocked in x) could be completed when the network is in state 
x; for yeY, rzy is the chance (or fraction of times) that if call c arises in 
state x it will be completed by being routed in the network so as to take 
the system to state y. The distribution {r,,,yeY} indicates how the call- 
ing-rate due to c is to be spread over the possible ways of putting up this 
call. Evidently, such a description of routing could be made time-de- 
pendent, and extended to cover refusal of unblocked calls as an option; 
we do not consider these possibilities here. The problem of choosing an 
optimal routing matrix R has been worked on at some length, in Refs. 
9 and 10; its relations to state reduction are described in forthcoming 
sequels” to the present paper. 


XVII. STOCHASTIC MODEL 


We now recall? a stochastic model for the traffic offered to a network. 
A Markov stochastic process x; taking values on S can be based on these 
simple probabilistic and operational assumptions: 


(t) Holding times of calls are mutually independent variates, each 
with the negative exponential distribution of unit mean. 
(ii) If uis an inlet idle in state xcS, andv =~ u is any outlet, there is 
a conditional probability Ah + o(h), \ > 0, as h — 0, that u at- 
tempt a call to v in (t, t + h) if x, = x. 
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(vt) A routing matrix R = (r,,) is used to choose routes, as follows: 
If c = {(u,v)} is a call free and not blocked in x, then the fraction 
of times that the system pass from x to yeA,, if c arises when x; 
= x Is just ryy. 
(tv) Blocked calls are declined, with no change of state. 
It is convenient to collect these assumptions into a transition rate 
matrix Q = (qxy), the generator of x;; this matrix is given by 
1 if ycB, 
Arzy if yeA, 
— |x| — As(x) ify = x, with s(x) = |F,| 
0 otherwise 


Cxy 


and the associated statistical equilibrium (or state) equations take the 
simple form 


[|x| + As(x)]px = XS py +X YX PyryxrxeS 
yeAx yeBy 
where {p,,x¢S} is the asymptotic distribution of x;. 
The remainder of this paper takes up the problem of correctly writing 
some reduced state equations for the stochastic model we have de- 
scribed. 


XVII MULTIPLICITY 


In an (unreduced) state x belonging to a structural equivalence class 
a there may be several calls in progress whose termination would carry 
the system into an equivalence class 8. For example, in Fig. 13, for a state 
x in a = 7 there are two calls whose ending would yield a state from 3, 
and there is one call leading to 2. Since we have assumed a unit hangup 
rate per call in progress, the downward transition rates between adjacent 
unreduced states is always unity; but as soon as we lump states into 
equivalence classes a,8,... these rates in effect add up. A similar situa- 
tion obtains for new calls: there may be several calls free and not blocked 
in a state xea whose completion could (under a suitable choice of route) 
lead to a state in equivalence class £. 

We call this phenomenon multiplicity, and we need an account of it 
in order to define transition rates between reduced states, and to cor- 
rectly use® and understand the reduced state equations. 

The “hangup” matrix H = (h,,) associated with a network v is a matrix 
of zeros and ones that tells which states can be reached from which by 
a hangup: 


Ney = XycBy = XxeAy 


It can be seen that H admits the automorphism subgroup G,, Le., 
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that 


hay = hegxy(gy) geG, 


FH is so-called because if we make the standard traffic assumption that 
all calls in progress terminate independently of any past history at unit 
rate, then H is precisely the part of the transition rate matrix due to 
“hangups.” 
The transition rate, due to hangups, from a state x to an equivalence 

class 8 is 

DL Axy 

yeB 
We claim that this number is the same for all x in an equivalence class 
a, and that it is the hangup rate from a to 8. There are two ways of 
proving this result, which will be put in the following form: 


Lemma 1: The numbers |B, ~ 6| and 


De Axy 

yeB 
are equal; they assume the same value for all x in an equivalence class 
a, and represent the hangup rate hag from a to B. 


Proof: Equality is obvious. The numbers are zero unless a covers @ in 
the partial ordering on E. Let xea, zea, z = gx where g is an automor- 
phism of (S, <). If now yeB, () 8, then gyeB,, = Bz, and gyef because 
@ is an equivalence class. Conversely if yeB, () 6, then g~!yeB, (1 8, by 
the same argument, so B, () 8 = g(B, (1 8). The result follows because 
£ is a bijection. 


Alternative proof: With xea, zea, z = gx as above, consider that where 
w is any element of 8 


| Box ‘a B| => XycBex 
yeB 








18| [8] 
vr weBgy h x Ww 
elie Gel 
|B 
= | SS hg = ey = Bel 8 | 
IG, | £G, vB 


The multiplicity question is more complicated for new calls than it 
is for hangups; this is because only one state can arise by removing a call 
c in progress, whereas possibly one of several could arise by completing 
a call c which is new and is not blocked; as a result, the calculation of 
transition rates between reduced states depends on routing decisions 
as well as on multiplicity. For by itself multiplicity will only provide the 
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calling rates which might take the state into a reduced state 8: whether 
a particular transition occurs depends, however, on what state x the 
system is in and what routing policy is being followed. We shall see, 
though, that if the routing is consistent in the sense that it routes anal- 
ogous calls analogously in equivalent states, then these calling rates will 
depend only on wx and §, i.e., on the reduced states in question. 
Mathematically this idea of consistency is expressed by the condition 
that the routing matrix R = (r,,) admit the group G,,; if c is a call free 
and not blocked in x, and 7¢G,, then in the terminology above, c and rc 
are analogous calls, x and 7x are equivalent states, and 7(A¢x) = A(ze)(rx); 
so consistency amounts to asking that we have ry = r(7x)(ry) for yeAcx, 
ie., that R admit G,. 

The matrix N = (n,,) associated with a network v is a matrix of zeros 
and ones that tells which states can be reached from which by adding 
a new call; it is obviously the transpose of the hangup matrix H: 


Nyy = hyx = XxeBy = XycAx 


Thus N admits the automorphism subgroup G,,, and we have this analog 
of Lemma 1, with a similar proof. 
Lemma 2: The numbers |A, (1) 6| and 

Lo xy (= Nag) 

yeB 
are equal; they assume the same value gg for all x in an equivalence class 
a, and they represent the number nag of calls free and not blocked in x 
which could be put up so as to lead to a state of 8. 

To calculate the actual transition rate from a state x into an equiva- 
lence class 6, we must use the routing matrix R. Notice that R is just N 
with enough of the ones reduced (but still = 0) so that r,, summed over 
yeA-,x is unity for each x and ceF,. 

For the traffic model assumed in Section XVII run according the 
routing matrix R the transition rate from a state x into an equivalence 
class 6 that intersects A, is 

> Ds ryy (= » rsy ) 


ec freeinx yeBCAecx yeBN\Ax 
c not blocked 


Notice that replacement of each r,y by 1 in this expression increases its 
value to precisely 


Nap = |Ax () B| = [{c:Bemr(Acx)}| 


Lemma 3: If R admits G,, then the numbers 
> r = 
yeliag 2 Pa) 
are the same for xea. 
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Proof: Like that of Lemmas 1 and 2. 


We put 
Aq = {@cE: B covers a in the partial order of E induced by <} 
By = {@cE: a covers @ in the partial order of E induced by S} 


These are the reduced analogs of the sets A, and B,, useful for writing 
the state equations. It can be seen that 


Aa = U @(Ax) 
By = U w(B,) 
Let us set 
hap BeBy 
Map BeAg 
Gap = 
—[As(x) + |x|] xea, a=8 
0 otherwise 


We can now informally describe the reduced equations as follows: 
when FR admits G,, they correspond to equilibrium equations for a 
Markov process 7x; on E whose transition rates up and down are Arag 
and hag respectively, with the obvious necessary convention on the 
“diagonal.” In terms of the notation introduced above, these reduced 
equations are 


~Pcdaa = » Pphgat r Pz PsrBa (2) 
BeAq BeBy 


with daa = —[As(x) + |x|] for any xea. 


XIX. REDUCTION OF STATE EQUATIONS 


We can now give precise explicit conditions under which the original 
“microscopic” equations of state can be exactly replaced by the less 
numerous “macroscopic” or reduced equations for the probabilities of 
equivalence classes. It will be shown that if the routing matrix admits 
G,, then the original and the reduced equations imply each other, and 
that the state probabilities {p,,x«S} are constant over equivalence classes 
a. Only the equilibrium case is considered. 


Theorem 5: Let the routing matrix R admit G,.' Then 
(t) The “microscopic” transition rate matrix 
Q = H +)R — diag[As(x) + |x|]xes 


also admits G,. 
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(it) The projection map 7:x — ax defines a “macroscopic” Markov 


process 7x; on the = — equivalence classes with transition rate matrix 
_ {hap t+ Nag a #6 
das = eas 
Qxx a B, XE 


(wii) If Da = dXoxea Px, then the equilibrium equations >"; Pxqx, and 
>a Padag imply each other. 

(iv) If{p,, wexS} solves the reduced equations, then {p,, x«S} defined 
by 

= Pa - 
Dx = Jal’ XEQ 

solves the original equilibrium equations. 
Proof: (t) is clear from the hypothesis that R admits G,,, from Lemmas 
1 and 8, and from the fact that s(x) = s(7x), |x| = |rx| for reG,. To 
prove (ii) it is enough, by (Lloyd’s) Theorem 1, to prove that 


se Qxy 


yeB 


are constant for xea. For a = 8 this follows from invariance of s(-) and 
|-| under 7eG,,. If 8 covers a in the induced partial order on 7S, then 


2. day = De tay = Map 
yeB 


Here we have used Lemma 3. Similarly if a covers 8, then for xea 
De Ixy = » hyy = [Br M Bl = 
yeB 
by Lemma 1; thus (iz) is sak 
To prove the reduced equations from the original ones, we sum the 
latter over xea, to get 


X pslrs(x) + [x[]= XL py tA XL YX dyrye 
xEea xea yeAy xew yeBy 


Since s(-) and |-| are constant over a, the left-hand side is just —podaa- 
For the first term on the right, argue thus: since hy, = 0 unless yeA, 


and 
Sa 


BeAg yeB 
sums over (at least) all yeA, if xea, we have 


» ae UF = 2, 2, Py 2 Iyx 


xea yeAy BeAg yeB xXEa 
=D DY pyl|By MO a| 
BeAa yeB 
mS Pphea 
BeAq 


REDUCTION OF NETWORK STATES 147 


using Lemma 1. Similarly, since ry, = 0 unless h,y = 1, we find 


 S Dylyx = aed Dd Py! yxltizy 


xea yeB, xea BeBu yeB 


= » Do Dy ds Pry 


BeBa yeB xEQ 


= D Psrga 
BeBa 

This proves the reduced equations from the original; the original ones 
can be proved from the reduced by reversing the identities used above. 
Neither argument depends on the fact (iv), to be proved next, that px 
is constant over xea. 

To prove (iv) we shall set py = p;x/|7x| in the unreduced equations 
and obtain an identity. Substitution in eq. (2) and multiplication by |a| 
gives, for xea, 


palds(x) + |x]] = lal © P™+arlal 5 P™ ry, 
yeA, | wy] yeB Bx |7 Ty 


= |a| x. PByaI 2 ys 


a3 


+ Ale| Ear 2 Ty 


ial 


la 
= 2. PB oy Dy Ney) 
BeAgq |8| yeB |G,| geGy er 


[a| 
+A > Pe 2 D Px(gy) 
BeBa eee ou 


|a| 
=) op »~ 
BeAq Hattie | geG, 


1 — fal 
RON Oo Dey a ae 
BBa i aa | eG, 


= oP DD hys 


BeAq Tl yeB xea 


ae re Psy 2 De Pix 
BeBa res B xen 


= DY pshgat XX dL DergBa 
BeAa BeBa 


Since the left side is —P a9 qa, the substitution has resulted in the reduced 








y(gx) 


Paes (gx )y 
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equations. Since the solution p of the original equations is unique, (iv) 
is established. 
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Elastic State of Stress in a Stalpeth Cable 
Jacket Subjected to Pure Bending 


By G. S. BROCKWAY and G. M. YANIZESKI 


(Manuscript received June 21, 1977) 


The stresses in the plastic jacket of a slightly bent telephone cable 
are analyzed within the linear theory of elasticity. The jacket is con- 
sidered to be bonded to the underlying corrugated steel by a flooding 
compound. The constraining effect of the steel results in a three-di- 
mensional state of stress that differs substantially from the predictions 
of elementary beam theory. For the thin Jackets typically used on 
telephone cables it is found that the stress state 1s essentially biaxial, 
the axial and circumferential normal stresses being at least an order 
of magnitude larger than the others. On the tensile side, the stresses 
are closely approximated (at any given point) by those in the well- 
known biaxial strip experiment, in which the principal stresses are in 
proportion by the Poisson’s ratio of the plastic. The compressive side 
ts likewise in biaxial compression, and there the flooding compound 
is subjected to tensile stresses even before the onset of any jacket 
buckling. The results confirm the validity of previous approaches to 
the effects introduced by imperfections and indicate further that the 
probability of spontaneous cracking ts increased by the adherence of 
the jacket to the soldered steel layer. 


I. INTRODUCTION 


The selection and development of plastic jacketing compounds for 
multipair cables depend to a large extent on cable behavior during 
bending. Cable jackets are expected to be relatively flexible for ease of 
handling and installation while at the same time surviving large strains 
(up to 15 percent) without cracking, splitting, or severe wrinkling. 
Temperature extremes encountered in the field render these criteria even 
more stringent. 

Recently completed analyses have led to easily performed laboratory 
tests for the screening of candidate compounds with regard to some of 


151 


these requirements. Now, for example, the relative influence of various 
sheath-grade plastics on the bending stiffness of cables can be evaluated 
by conducting ordinary tensile tests.! In addition, the relative sensitivity 
of compounds to low temperature and high strain-rate cracking can be 
determined through impact tests on notched specimens. Still, there 
remains the observation of slow crack growth at high temperatures 
during bending and the occurrence of wrinkles in cable jackets during 
duct installation at low temperatures. None of these phenomena nor how 
they are affected by cable parameters such as jacket moduli and thick- 
ness, flooding compound tackiness, and the depth of the corrugations 
in the underlying steel, is presently understood. 

This paper is devoted to a study of the state of stress in the jacket of 
a Stalpeth cable subjected to classical pure bending. The problem is 
treated within the framework of the linear theory of elasticity, which 
supposes small strains and rotations and an elastic material. Although 
the bending strains in telephone cable are frequently large, our analysis 
is intended to provide insight into the circumstances at incipient 
cracking, buckling, or yielding of the cable jacket. Standard techniques 
from linear viscoelasticity theory2? can be applied to the elastic results 
given here to account for the time dependence inherent to plastics. 

The primary emphasis in the present investigation is directed toward 
estimating the effect on the stress field in the jacket resulting from the 
constraint imposed on its inner surface by the underlying cable structure 
(see Fig. 1 for a detail of the Stalpeth construction). For slight bends, 
the influence of the soldered steel shield dominates that of the wire core, 
and we shall, therefore, consider the pure bending of a plastic jacket 
bonded to a corrugated metal shell. The analysis is further simplified 
by the realization that the corrugation wave length used in Stalpeth cable 
is one to two times larger than the jacket thickness, while the valleys 
imprinted on the jacket’s inner surface are reasonably shallow because 
of the presence of the flooding compound. We are thus afforded the 
privilege of averaging field variables over a corrugation wave length to 
arrive at a boundary value problem involving a uniform cylindrical ge- 
ometry. The amplification of jacket stresses created by the corrugation 
imprints can be deduced from the results obtained here together with 
published concentration factors in the usual way. . 

We begin in the next section by formulating the relevant three-di- 
mensional field equations and show that they can be reduced through 
a change in dependent variables to plane strain equations. The approach 
taken is reminiscent of the scheme used in elementary elasticity for the 
St. Venant-bending of cylinders with irregularly shaped cross sections, 
but differs in that here the strain field, rather than the stress field, is, 
supposed to conform to the elementary Bernoulli-Euler theory. . 

Next, the boundary conditions at the jacket-steel interface are con- 
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POLYETHYLENE JACKET s- 


N 
SOLDERED CORRUGATED STEELS \ 
WITH THERMOPLASTIC FLOODING “\ 
\ 






CORRUGATED ALUMINUM SHIELD ~ 
N 


Fig. 1—Stalpeth cable. 


sidered in detail. The steel is assumed to have negligible stiffness in the 
corrugation direction and to otherwise obey the usual hypotheses on the 
deformation of thin shells. Integration over a corrugation wave length 
then permits the stresses and displacements of the jacket to be related 
to those of the steel shield. 

The appropriate (plane strain) boundary value problem for the stress 
state in the sheath having been set up, we find that it admits an ele- 
mentary solution in closed form through the introduction of an Airy 
stress function. The full three-dimensional stress and displacement fields 
_are then calculated. 

When these results are applied to Stalpeth cables as presently man- 
ufactured, it is found that, on the tensile side of the cable, the jacket is 
essentially in a state of biaxial tension. That is, the ratios of the axial and 
circumferential normal stresses to the shear and transverse (thick- 
ness-direction) stresses are.of the order of the diameter-to-thickness ratio. 
Furthermore, if F denotes the Young’s modulus of the plastic jacket, D 
the cable diameter and p the bend radius, then the longitudinal and hoop 
stresses are shown to be in approximate proportion to the bending 
stresses E-D/2p from elementary beam theory. Of particular interest is 
the conclusion that these constants of proportionality vary appreciably 
only with the Poisson’s ratio of the jacket, at least for the ranges of cable 
size, corrugation geometry and jacket moduli encountered in Stalpeth 
applications. With Poisson’s ratio chosen in the typical range for plastics 
(0.3-0.5), the longitudinal stress varies from about 110 to 130 percent 
of the beam theory stress, while the hoop stress factor ranges from 30 
to 45 percent. The stress state on the tensile side of the cable is thus 
closely approximated by the biaxial strip experiment discussed at the 
end of the paper. The importance of biaxiality in bending has been 
recognized previously.4 

All of the results outlined above for the tensile side of the cable apply 
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as well on the compressive side, except that the stress state is, of course, 
one of biaxial compression. It is also observed that the flooding com- 
pound is subjected to a tension on the compressive side that exceeds that 
which currently used compounds are likely to support, even for very 
small bend radii. 

The implications of these results on the cracking to which we alluded 
earlier are discussed qualitatively in the final section of the paper. For 
the time being, we remark only that the biaxiality tends to increase the 
likelihood of spontaneous cracking. Finally, it should be mentioned that 
the jacket stresses generated in the bending of cables other than Stalpeth 
can differ drastically from those obtained here. Indeed, the biaxial state 
of stress produced by the constraining effect of the soldered steel layer 
would not be present in those designs that allow relative motion between 
the jacket and underlying metallic layers. 


ll. REDUCTION OF THE PROBLEM TO ONE OF PLANE STRAIN 


Consider a cylindrical shell of inner radius r; and outer radius ro. 
Referring to Fig. 2, choose cylindrical coordinates (r,6,z) in the obvious 
way, and let the shell be subjected to pure bending in the y-z plane. The 
center line of the bent shell is then a circle of prescribed radius p. 

Recall the field equations of linear elasticity in cylindrical coordi- 
nates:° 

Strain-displacement relations 

















Ou; 1 OUg Ou, 
= ,¢¢=—-(—+u,),¢ =— 7; 
= or 8 r le ) Oz (1) 
1 /du, OU du, . OU; ou, 1 du; 
2 3 = + — . 2 = + : 2 =o _ 
pian ey us) ae oe oe Ba ee ee Oe 


Equations of equilibrium 








Oo,  _1Oorg | o- — 06 
ies -—————= 0 2 
or sr 06 y 22) 
re) 10 ey 
—2rh 4 = SOG | SO 9 (2b) 
or r 06 
O0rz I 0092 Ooz = (2c) 
or r o0 Oz 
Hooke’s law 
Qu(1 + ve, = 0, — vlog + oz) (3a) 
Qu(1 + veg = a9 — v(o, + oz) (3b) 
2u(1 + vyez = oz — vo, + a9) (3c) 
2uere = Oro, 2MErz = Frz, 2MEoz = Coz (3d) 
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—P < 





Fig. 2—Choice of cylindrical and Cartesian coordinates. 


Here, the symbols u, ¢, and o, suitably subscripted, stand for components 
of displacement, strain, and stress, respectively. The constants » and 
v represent the shear modulus and Poisson’s ratio for the cable jacket, 
while the coefficient 


E = 2u(1 + v) (4) 


of ¢,, €9, and ¢, in the first three equations of (3) is Young’s modulus. 

Next, make the assumption that the annular cross sections remain 
plane and normal to the shell’s axis during bending.* It then follows from 
a familiar geometric argument’ that 

r sin @ 
Erz = €9z = 0, eg = —_ (5) 
p 

With no consequent loss of generality, we assume further that the cross 
section at z = 0 remains in the x-y plane and undergoes no rotation about 
the z-axis nor overall rigid translation in the x-y plane, so that 


u,(r,0,0) = 0 (6a) 


U19(1-9,0,0) = wy (ro ; 0) = ug(ro,7,0) = 0 (6b) 


We show next that the formulae (5) and (6a) require the in-plane 
stresses a,, a9, and o;9 to be associated with a state of plane strain. That 
is, we establish the existence of plane displacements ti, and tg, inde- 
pendent of z, that generate strains é,, é9, and é-9 obeying the plane strain 

* The viewpoint here is similar to that of the traditional St. Venant “semi-inverse” 


method.® The solution thus obtained is justified by exhibiting stress distributions on the 
ends of the cylinder that support the assumed deformation. 
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form of Hooke’s law:° 


Qué, = (1 — v)o, — vag (7a) 
2Queég = (1 — v)o9 — vo, (7b) 
2ue-e = Org (7c) 
We first confirm (7) and then prove the existence of the requisite dis- 
placements. 
Since ¢, is known from (5), we have from (8c) that 


oz = Qu(1 + v)e, + v(a, + 9) (8) 
Thus, o, may be eliminated from eqs. (3a), (3b) for e,, eg to obtain 
€r = €p — VE, €9 = €9 — VE (9) 
where é,,é9 are given by (7a), (7b). The last of (7) is, of course, satisfied 
by taking 
€ro = Ero (10) 
To see that the strain field é,, é9, and é,-9 is indeed generated from (7) by 


in-plane displacements tw, and tg, independent of z, observe that the 
displacement field* 


= (pr2+2z2)sind _ (vr2—z2)cos6 _ rzsin#@ 
LS Bp SS, 8S LD 
2p 2p p 
has an associated strain field 
€r = €9 = —Vez, Ez = Ez, Eg = Ez = €9z = O (12) 
Thus, the defining equations 
Uy = Ur + Uy, Ug = Ug + Up, Uz = UZ, + UL, (13) 


for u,, tg, and uz, together with (1), (9), and (12), reveal that the dis- 
placements U,.,%,9 satisfy the strain-displacement relations for €,,€9, €-9. 
Moreover, when (13) is combined with (7c), (5), and (12), there re- 
sults 

Ollr | Oz _ RL 0, Ou: 
Oz or oz r OO Oz 


Integration of (14) subject to the constraint (6a) yields 








=0 (14) 


and the desired result is established. 


* This displacement field (11) is that which the shell would exhibit were it hollow (see 
Sokolnikoff,® Article 32, for example). 
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The three-dimensional problem is now reduced to determining the 
plane-strain elastic state with displacements (t,,t/9), strains (é,,€9,é,9), 
and stresses (0;,09,0;9), subject to boundary conditions on the cylindrical 
surfaces r = r; and r = ro. The remaining stresses o,z,0¢;2,09; are then 
provided by (8) and 


Orz = O92 = O (15) 
which follows from (5) and (8c). 


I. THE BOUNDARY CONDITIONS 


As previously mentioned, the intent of this paper is to account for the 
influence of the soldered steel shield on the plastic jacket. Since the 
longitudinal (z direction in Fig. 3) stiffness of the corrugated steel is small 
for small longitudinal extensions,® the constraint imposed on the inner 
surface of a Stalpeth jacket is essentially confined to the x-y plane, 
provided the jacket field variables are interpreted as averages over a 
corrugation wave-length. 

Moreover, since the corrugation depth H (see Fig. 3) is small compared 
to the radius r; of the cable, the shield will deform in approximate ac- 
cordance with Kirchoff’s hypothesis of classical shell theory.? In par- 
ticular, denoting by u and vu the circumferential and radial displacements 
of the midsurface of the shield (see Fig. 3), one has the relations?® 


ur(é,0) =v (16a) 





2é dv 
(£0) =ut —— 16b 

SET oe («- 3) ee) 
for the radial and circumferential displacements of a particle at a dis- 
tance £ from the midsurface (again, see Fig. 3). Requiring the displace- 
ments of the steel to conform to those of the jacket at £ = H/2 results 
in 


v(O) = u-(7ri,9) 


17 
Hou, ou 


2h; 
eee = ug(r;,0) + 
Ce 7) ata) 2r,; -H 00 


Next, take the shield to have an (in-plane) modulus E, and use (16), 
(17) to compute the tension T and moment M in the shield averaged over 
a corrugation wave length @ (refer again to Fig. 3). There results 








(r;,0) 





7 (0) = 224 oo(r;,0) a Mie) 
(18) 
2 
M(0) = ae | trl 8) + a (r,0) | 
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Z- 
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(a) (b) SECTION A-A 


(c) 





Fig. 3—Corrugation details. 


where A and J are respectively the area and area moment of inertia about 
the midsurface of a corrugation wave length (see Fig. 3a), and ¢g is the 
circumferential strain in the jacket. 

Finally, recall that for the shield to be in equilibrium under the in- 
terfacial loads o,,o,9 exerted by the jacket, T and M must satisfy!° 


(:- 3) = T-—_ atr=r; 
ey r; - H/2 do? 

Hy _ aT __1 dM 
(25) a ae a 


(19) 
atr=rj; 


The sought-after boundary conditions at the inner surface of the jacket 
can now be extracted from (18) and (19) with the aid of (5), (9), (11), and 
(13). Neglecting terms in H/2r; compared to unity, one has 








E,A AH 1 d2M E;A vr; . 
=I; er €9 — Gaza aa) =~ sin 0 (20a) 
E,A 0é9 AH 1\dM_ E,Apr; 
Org + Sy (SS SS 6 20b 
eee ap G -) dd ep. ve 
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where 








EI 07t, 
= — 4 20 
r;2e ( J ae) ( °) 
The remaining boundary conditions 
Or = og = 0 atr=Przo (21) 


insure the absence of load on the exterior of the cable. 


IV. SOLUTION OF THE PROBLEM 


According to the plane theory of elasticity, there exists an Airy stress 
function ¢ generating the stresses o,, o9, and o,9 through the rela- 
tions® 


Or 


10¢. 1 0%¢ 
== + 
ror. r? 062 


(22) 


re) (- ot) 
org = - a (7 
ue or \r 06 


In addition, ¢ must satisfy the biharmonic equation 
V4¢ = 0 (23) 


Michell?! has given a general form for Airy’s function for annular do- 
mains subject to the requirement that the stresses and displacements 
be single-valued.* Examining Michell’s solution (see also Fung,° p. 246) 
and the boundary conditions (2), it seems worth trying ¢ in the form 


@ = (ar? + b/r) sin 6 (24) 


To determine the constants a and b in (24), substitute first into (22) 
and impose the boundary conditions (21) to conclude that 


b=ar,4 
and hence 
o, = —2af(r) sin 0, og = 2ag(r) sin 0, o-6 = 2af(r) cos 0 


f(r) = r(ro4/r4 — 1), g(r) = r(ro4/r4 + 3) 


Before enforcing (20), we compute the displacements and strains gen- 


(25) 


* Although Michell’s argument for the generality of his solution is sketchy, his result is 
nevertheless correct. This can be established directly with the aid of Muskhelishvili’s!2 
complex valued stress functions. 
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erated by the stress state in (25). Toward this end, substitute (25) into 
the plane strain form of Hooke’s law to arrive at 


A 


ran ‘ 
é- = ——6,(r) sin 0, 9 = — Ep(r) sin 9, é-9 = —= f(r) cos 0 
a 


= 


(26) 


5 ro4 i i 

Ear (“2+ 4-1), 6o=r (“2+ 3-40) 
r ic 

If eqs. (26) are now incorporated into the strain displacement relations 

(1) and an elementary integration is performed, it is found that 


4 an ~ 
tn = 5 | = tv — re + A sin#@+B cos 0 
Qu r2 
oe (27) 
Ug = -2 | 4 (—4y)r2+ Al cos6—Bsiné+ Cr 
QuLr2 
where A, B, and C are constants fixed by the condition (6b). In fact, (27), 


(13),-(11), and (6b) give 
B=C=0,A= (% + 2(2-3)) 10? (28) 
pa 


We are now in a position to determine the unknown constant a. Ob- 
serve first that (27) and (20) require that the moment M in the steel 
vanish identically. Thus, the second of (20) follows from the first pro- 
vided 


org 
od 


a condition that is met for all r by the stress field in (25). Equations (25), 
(26), and the first of (20), therefore, reveal that the boundary relations 
(20) are all satisfied provided 





=o, atr=r; 


uvr;/p Quer; 
Sate eee Fae 29 
° 6 o(r;) + Sf(r;)’ E,A | (29) 


The result in (29) may be rewritten in the form 


pv 
a= 7 , 
} (30) 
A=3- a 474s [T 1] 
rj 
with the aid of (25) and (26). The stress field is now completely deter- 
mined since by (8) and (25), 


sin 6 


4 
r= (1 rene —\ 2ur 
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or, using (30) and (4), 
2 
eee oe a 
ACL +v)4 p 


By the same token, the displacements are given by (11), (13), (27), (28), 
and (30) in the form 





2 2 2 2 
“ea 2 eas Geer A+ 207-3) | sin 0 
Eger Dhan 


Qu Lr? 
aro? [ro2 r2 
= ——* | *- 4+ (5-4 -\) — (32) 
is 2u [ ( : yn 


2 
+2240 4200-2) cos 6 
Vo 


The solution to the original three-dimensional problem has thus been 
found. Although there is no available uniqueness theorem encompassing 
the present problem (the boundary conditions (20) are nonstandard), 
it can be shown from Michell’s!! representation for the Airy function 
that the solution obtained here is unique except for certain peculiar 
choices of the elastic and geometric parameters of the jacket and shield. 
The physical significance of these instances of non-uniqueness is related 
to the important sheath-buckling phenomenon that has been observed 
in telephone cables. Since the present analysis is insufficient to ade- 
quately describe the crucial z-dependence of the ripples, we shall explore 
this issue no further in this paper. 


V. DISCUSSION OF THE STRESSES IN THE JACKET 


Observe first [eq. (31)] that the axial stress o, varies linearly with the 
distance y = r sin 6 from the neutral plane (y = 0), just as in the ele- 
mentary Bernoulli-Euler theory. In addition, (25) reveals that o, and 
og are likewise antisymmetric about the neutral plane, co, being com- 
pressive and og tensile when co, is tensile and vice versa. The shear stress 
og, aS expected, is symmetric about the neutral plane. 

The variation in o, (and hence o;») and o» across the thickness of the 
jacket is illustrated in Fig. 4, where the dimensionless functions 


FBAR = —f/4r,, GBAR = g/4r, (33a) 
[Eq. (25)] are plotted against the dimensionless radius 
RBAR =F =r/ro (33b) 


for 0.5 <7 < 1. It is seen that as long as the jacket thickness is less than 
10 percent of the cable radius (0.9 < r < 1), the transverse stress a; is less 
than 10 percent of the hoop stress og, which is itself essentially constant 
across the jacket. 
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RBAR 
Fig. 4—Radial dependence of o, (FBAR) and og (GBAR). 


With a view toward examining the influence of the properties of the 
jacket on the stress state, let ¢,™ and o9” stand for the maximum values 
of the axial and hoop stresses, respectively. Then (25) and (31) together 
with (30) and (4) give 


Ev 
ms Ly 2, Trg a GAR i 
a9” = a9(ri,7/2,z) (+ yap 27?) 
(34) 
Zz 2" 09 ) (1 4+ pv) 0 


The ratio of the maximum bending stress in the cable to the prediction 
of elementary beam theory is thus 
m 2 
Ge bs gig cae 
Ero/p (1+ v)d 





SIGMAX = (35) 
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Further, the degree of biaxiality created in the jacket by its constrained 
inner surface is described by the parameter 
m S 
ALPHA = 22 = vg(ri)/ro 

on™ (L+v)d\ + 4v? 

In Fig. 5, the dimensionless stresses SIGMAX and ALPHA are plotted 
versus Poisson’s ratio v for a 3-inch diameter Stalpeth cable conforming 
to the current Bell System design. In computing the parameter S that 
enters the formula (30) for A, the area A of a half-corrugation wave length 
has been approximated by 


(36) 





Aa 
A=Vm?+1 fis, m= s (37) 


which is the area of the parallelogram shown in Fig. 3c. Also, the modulus 
E's was assumed to be that of steel and the jacket modulus F was sup- 
posed to be 45,000 psi, representative of low-density polyethylene at 
room temperature. !3 

The curves shown are quite insensitive to variations in jacket modulus, 
SIGMAX changing by less than 44 percent and ALPHA by less than 2 
percent when F£ is increased to 90,000 psi. The results indicate that for 
a Poisson’s ratio v = 0.35, the bending stress in the jacket is about 12 
percent higher than the prediction of the elementary theory while the 
hoop stress is roughly 1/3 of the axial stress. 

The formulae (35) and (36) can also be used to investigate the effect 
of cable size on the state of stress in the jacket. Once the cable diameter 
is chosen, all other geometric parameters (e.g., jacket thickness, corru- 
gation height) are fixed. This information was incorporated into (35) and 
(36) resulting in equations relating SIGMAX and ALPHA to cable diameter 
and the elastic constants of the jacket. For the range of Stalpeth cable 
currently manufactured, SIGMAX and ALPHA are independent of cable 
diameter. 

The conclusions already reached as well as others of interest may be 
inferred more easily from the approximations of the previous equations 
arising from assuming the jacket to be thin. Toward this end, let 


t=Pro—T; (38) 
so that, for any n 
An n t t\2 
ee, = (1-<) =1-n—+0|(-) | ast/r, > 0 (89) 
ee To To Yo 


If the approximation in (39) is applied to the formula (30) for A, one 
gets 


N= 4(1 — v) +: 4t/ro + 4S8t/r. + O[(t/ro)?] 
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Fig. 5—Dependence of bending and hoop stresses on Poisson’s ratio (cable diameter 
= 3 inches, EF = 45,000 psi). 


But from (30), (37), and (39) 


SB tha (1-1) ool 
ro Esvm2+1tsro EstsVm24+1 Ps Ty 


whence 





had (1-vt So ) + O(t/ro) 


1lt+p (40) 
— w+ vt Et 
9 m2 + LlEsts Vm24+1E£Ests 
Equations (40) and (35) give 
2 
v 
=) 4 FE OG (41) 
SIGMAX i- +S, (t/ro) 


Similarly, 
g(ri) =A4 O[(t/ro)?] 
r 


oO 
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SIGMAX 





ie 6—Dependence of bending and hoop stresses on the modulus parameter S, (v = 
0.35). 


so that, using (40) as well, ALPHA may be approximated from (36) by 


v 
+S, 


The dependence of the nondimensionalized stresses on Poisson’s ratio 
is thus seen from (41) and (42) to have the form indicated in Fig. 5. 
Moreover, (41) and (42) exhibit no explicit dependence on cable radius, 
though there is a slight implicit dependence entering through So, since 
t and m are governed by cable diameter. 

The approximate relations (41) and (42) have been used to generate 
the curves in Fig. 6. The parameter S, for present cable designs is very 
small (1072 to 1073) which accounts for the insensitivity of the curves 
in Fig. 5 to variations in jacket modulus in contrast to those in Fig. 6. In 
fact, the ratio Et/E,t, is so small for all conceivable cable applications 
that S, may as well be neglected. If this is done, then (41), (42) can be 
used in conjunction with (25), (31), (35), and (36) to provide 


E 
1-p? 

ro sin 0 
eS 


p 





ALPHA = + O(t/r,) (42) 





Oz ~ €, 09 ~ Voz 


(43) 


All other stresses zero 


as a good approximation to the jacket stresses. 
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VI. IMPLICATIONS ON CABLE PERFORMANCE 


One important result that emerges from the present investigation is 
that (43) (for fixed 0,0 < 6 < x) is the stress state that arises in the so- 
called “biaxial strip” experiment (Fig. 7). When a strain ¢ > 0 is imposed 
in the y direction, the strip is in plane stress (o, ~ 0) throughout and 
plane strain («, ~ 0) in a region near the center. The circumferential and 
longitudinal directions in the cable thus correspond, respectively, to the 
x and y directions in Fig. 7. 

The approximate plane strain condition in the bent Stalpeth cable 
is created by the adherence of the jacket to the steel. Should this con- 
straint remain intact during continued bending of the cable, the analogy 
with the strip experiment also continues. In that event, the elongation 
at break in the biaxial strip experiment would be an important material 
parameter. Another biaxial experiment* (equal principal stresses) in- 
dicates that some low-density polyethylenes exhibit an ultimate elon- 
gation biaxially that is substantially reduced from the uniaxial value, 
say to 20 percent or perhaps even less at lower temperatures or higher 
rates. Since jacket strains typically reach 15 percent during duct in- 
stallation, failure due simply to biaxiality is indeed a concern. The strip 
experiment should, therefore, be instituted as a materials screening 
test. 

Even if the material exhibits a high elongation in the biaxial strip 
experiment, however, the cable jacket may fail because of the localization 
of deformation in the neighborhood of imperfections, such as corrugation 
imprints or surface scratches. In these cases, failure occurs at bend radii 
for which the present analysis is applicable to points in the jacket away 
from the imperfection. Thus, if the geometry and orientation of the 
imperfection are known, concentration factors may be used to determine 
under what conditions the transition to highly localized deformation 
occurs. The critical parameters at which this transition occurs are cur- 
rently found from impact tests on relatively narrow notched bars.! The 


4 
\ \ = -t 
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Fig. 7—Biaxial strip experiment (clamps pulled in y direction). 
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biaxiality evidenced here should have little effect on the conclusions 
drawn from such a procedure, since the notches create a complex state 
of stress which is largely independent of width.4 

In addition, the growth of sharp internal or surface flaws under the 
influence of the far-field stresses calculated here can be analyzed with 
the aid of a recently developed viscoelastic theory!® coupled with stress 
intensity factors available from the literature.!® Theoretical failure times 
would then be compared against the duration of loading to give a prob- 
ability of cracking failure for a given flaw population. Here, we remark 
only that the probability of premature fracture is obviously enhanced 
by the constraining effect of the steel, since the severity of loading for 
those flaws not oriented in the @ direction is increased by the biaxial- 
ity. 

In addition to the effect on fracture performance, biaxiality can also 
have a substantial effect on the yielding behavior of plastics.417 But in 
view of the fact that cable jackets, including Stalpeth jacket, are im- 
printed with corrugation valleys, localized yielding behavior is more 
relevant. 

The failure mechanism on which the present analysis sheds the most 
light is that of sheath buckling. Firstly, on the compressive side of the 
cable, (25) indicates that the jacket exerts tensile loads on the flooding 
compound. Estimates from the preceeding formulae reveal the magni- 
tude of this stress to be as high as 200 psi during installation. It would 
therefore seem that only an exceptional flooding compound would suc- 
cessfully restrain the jacket at low temperatures and suggests a need for 
the mechanical characterization of these compounds. The situation is 
further aggravated by the biaxiality and the fact that the stresses vary 
in proportion to the modulus. This accounts for the occurrence of 
buckling at low installation temperatures. The temperature at which 
buckling will occur for a given duration of loading can be estimated from 
this analysis and elementary viscoelasticity theory. This buckling 
analysis is presently being pursued. 


Vil. CONCLUSIONS 

We have shown that: 

(t) The state of stress in the jacket of a bent telephone cable is es- 
sentially biaxial and constant across the thickness. 

(ii) At any fixed point in the jacket of a Stalpeth cable, the stress state 


is essentially that in the biaxial strip: 
Hips a a fo sin 9 
ae a a F 


with all other stresses and strains negligible. Thus, the maximum 
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bending stress is 110-130 percent of the elementary beam theory pre- 
diction while the maximum hoop stress is 30-45 percent of the bending 
stress. 

(tii) The flooding compound on the compressive side of the cable is 
subjected to tensile and shear stresses as high as 200 psi. 

(tv) Practical variations in cable design parameters such as corru- 
gation depth and frequency, steel thickness and jacket thickness have 
no significant effect on the global stresses due to bending. 


Having the results of this as well as previous investigations in mind, 
we have concluded that: 


(t) The currently used techniques! for evaluating the notch sensi- 
tivity of plastics are applicable to the bending conditions encountered 
in cable installation. 

(it) The failure strength of flooding compounds can play an impor- 
tant role in preventing jacket buckling. For this reason, any anticipated 
change in flooding material should be thoroughly examined with respect 
to tensile strength and adhesive strength. 

(iit) The probability of spontaneous cracking is increased by the 
adherence of the jacket to the soldered steel layer. 
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A Generalization of Takagi’s Theorem on 
Optimal Channel Graphs 


By F. R. K. CHUNG and F. K. HWANG 


(Manuscript received June 6, 1977) 


A channel graph, also called a linear graph,°* is a multistage graph 
with the properties that (1) each of the first and the last stages consists 
of a single vertex (denoted by I and O respectively); (ti) for any vertex 
v ~ Tor O, vis adjacent to at least one vertex from the preceding stage 
and at least one vertex from the following stage. In a switching network, 
the union of all paths connecting a fixed input terminal to a fixed 
output terminal can usually be studied as a channel graph by taking 
each switch as a vertex. In comparing the blocking probabilities of two 
channel graphs with the same number of stages, we say one is superior 
to another if its blocking probability is less than or equal to that of the 
other under any link occupancies. Takagi proved a basic theorem in — 
showing one type of channel graph is superior to another. In this note 
we present a more powerful result which includes Takagi’s theorem as 
a special case. 


|. INTRODUCTION 
A graph is called a multistage graph if its vertex set can be partitioned 
into subsets V,,..., V;, for some number s, and its edge set into subsets 


Ky,..., £,-1 such that E; connects V; with V;+;. A channel graph, also 
called a linear graph,° is a multistage graph with the properties that (z) 
each of the first and the last stages consists of a single vertex (denoted 
by J and O respectively); (zi) for any vertex ~ J or O, v is adjacent to at 
least one vertex from the preceding stage and at least one vertex from 
the following stage. In a switching network, the union of all paths con- 
necting a fixed input switch to a fixed output switch can usually be 
studied as a channel graph by taking each switch as a vertex and each 
link as an edge. 

In a switching network, a link can be in either one of two states, busy 
or idle, depending on whether it is part of a connection carrying a call. 
A path from a fixed input switch to a fixed output switch is blocked if 


171 


(a) (b) 
Fig. 1—Channel graphs in Takagi’s theorem. 


it contains a busy link. A pair of switches is blocked if every path between 
them is blocked. The same notion of “blocking” applies to the study of 
channel graphs and therefore we can talk about the blocking probability 
of a channel graph. 

A series-parallel channel graph is a channel graph which is either a 
series combination or a parallel combination of two smaller series-parallel 
channel graphs with an edge being the smallest such graph. Channel 
graphs which are not series-parallel are often called spider-web channel 
graphs. Recent studies have shown, either by analysis or by simulation 
(see Refs. 3 and 9, for example), that spider-web channel graphs can 
sometimes significantly reduce blocking probabilities over series-parallel 
channel graphs for given switching network hardware. In particular, 
Takagi®9 gives a useful theorem which compares the blocking probability 
of the spider-web channel graph in Fig. 1a and the series-parallel channel 
graph in Fig. 1b. 

In Fig. 1a, the connection between the two middle stages can be viewed 
as a complete bipartite graph on m and n vertices. In Fig. 1b, the con- 
nection between the two middle stages can be viewed as a matching of 
m pairs, each pair joined by n multiple edges. 

The above theorem is the basis of Takagi’s work®.? on optimal channel 
graphs which has been widely quoted in the literature (see Refs. 1, 2, 4, 
6, 7, 10, 11, and 12, for example). In this note we present a more powerful 
result which deals with a much larger class of channel graphs and in- 
cludes Takagi’s theorem as a special case. 


ll. TAKAG’S THEOREM 


Takagi’s comparison of two 4-stage channel graphs as shown in Fig. 
1 actually has broader applications than it appears. The extension is 
made possible by interpreting each edge in the 4-stage channel graphs 
as a reduction of a multistage graph. The only requirement is that the 
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(a) (b) 
Fig. 2—Two 6-stage channel graphs. 


multistage graphs represented by the edges in the same set E; are iso- 
morphic. For instance, the two 6-stage channel graphs shown in Fig. 2 
can be reduced to the two 4-stage graphs shown in Fig. 3. 

Note that a vertex in the 4-stage graphs can represent a group of 
vertices (from the same stage) in the 6-stage channel graphs. Further- 
more, two disjoint edges in the 4-stage channel graphs can come from 
two nondisjoint subgraphs of the 6-stage channel graphs. Finally, an edge 
in the 4-stage channel graphs can have more than one state where a state 
is basically a distinct subset of nonblocking paths in the corresponding 
multistage graph. 

By associating a probability distribution to the joint states of the edges, 
the blocking probability of a channel graph can be computed. Let D be 
a collection of probability distributions on the joint state of the edges. 
Then an s-stage channel graph G is said to be superior to another s-stage 
channel graph G’ with respect to D if given any member of D, the 
blocking probability of G never exceeds that of G’. Let P(X;), i = 1, 


Fig. 3—Two 4-stage channel graphs. 


E3 


E, 
COS 
oe 


E E% E3 


ee 


Fig. 4—The mapping from Fig. 2 to Fig. 3. 
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Fig. 5—Channel graphs in main theorem. 


...,«, denote the probability that an edge of FE is in state X;, let P(Z;), 
j=1,...,y, denote the probability that an edge of Eis in state Z;, and 
let Y(t,j7) denote the blocking probability of a path from I to O which 
contains an edge of FE; in state X; and an edge of E3 in state Z;. Finally, 
let S be a joint state of the edges in E3. Then Takagi proves: 


Takagi’s theorem. The channel graph of Fig. 1a is superior to the channel 
graph of Fig. 1b for arbitrarily given S, P(X;), and Y(i,j) under the fol- 
lowing assumptions: 


_ (i) The states of the edges are independent. 
(ti) n =m. 


In the next section, we give a generalization of Takagi’s theorem. 


lil. THE MAIN THEOREM 


Consider the two channel graphs in Fig. 5. 

In Fig. 5a, each vertex in stage 2 has n’ edges connected to n’ distinct 
vertices in stage 3, and each vertex in stage 3 has n edges connected to 
n distinct vertices in stage 2. Furthermore mn = m’n’,n’ < mandn S 
m’. Figure 5b is the same as Fig. 1b. 

Using the same notation as in Section II, we now prove: 


Main theorem. The channel graph in Fig. 5a is superior to the channel 
graph in Fig. 5b for arbitrarily given S, P(X;) and Y(i,j) under the fol- 
lowing assumptions: 


(t) The states of the edges are independent. 
(tt) m’ =m. 


Before we prove the theorem we state a lemma proved by Takagi® 
which is a generalized version of Hélder’s inequality. 
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Lemma. If a;; 2 0, and 


- a1 
1 6; 


for b; > 1, and d; 2 0, then is following inequality holds: 


TMs 


> Qv TI Ul (a3) < TI (= nal)” (1) 


1=1 jJ=1 j=l 


Proof of the theorem. Let S be an arbitrary state on the set of edges 
between stage 3 and stage 4 in Fig. 5. Suppose under S, z; edges between 
stage 3 and stage 4 are in the state Z;,j = 1,...,y. In Fig. 5a, let wz; be 
the number of edges joining the kth vertex in stage 2 to a vertex v in stage 
3 such that the edge between v and O is in the state Z;. Then it can be 
easily checked: 


D wp = Nz; (2) 
k=1 
y 
Do Wa = 1’ (3) 
jJ=l1 
and 
y 
2, 4 =m (4) 
J=1 


Let Y(i,j) denote the blocking probability of a path from J to O which 
contains an edge in the state X; between stage 1 and stage 2 and an edge 
in the state Z; between stage 3 and stage 4. Furthermore, let P(X;) de- 
note the probability of X; and let x be the number of possible states X;. 
Then the blocking probabilities of the channel graphs in Fig. 5a and 5b 
for the given state S, denoted by B, and Bz, are 


Ba= TI | > POX) tl YG" 


k=1 


and 
o= Tr] & Pexpyaa | 
Using (2), (3), and (4), we have 
wpj/n 
y m’ x Rj/n m! wpRj/n 
= ft { fh Ec iota |= ft {ft | > pexnves] "| 
=] k=1 (j=1 
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(a) (b) 
Fig. 6—Example 1. 


(a) (b) 
Fig. 7—Example 2. 


(a) (b) 


Fig. 8—A counterexample. 


(a) (b) 


Fig. 9—A reduction of Fig. 8. 
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If m’ = m, then n’ = n. By substituting A; = P(X;), aij = Y(t)” and 
b; = n/w;; in (1), we obtain 
By> TI { x: P(X) tl | YGie} = Be 
k=1 
If m’ > m, thenn > n’. Define 





Y@yt+1)=1 and byyi= 
—n 


Then 6,41 > 1 and 


ytl : , 
> ee oe eae 
j=1 0; jar n n 
Now 
m’ fytl x is bj 
Bs= TE LT & Pri? | | 
k=1 (j=1 Li=1 


m’ x yt1 
> Tt |X PO) TT vain 
k=] j=l 


= fl [ P(X) ql [Yt] = Be 
where the inequality is obtained by making a rene substitution as 
before. The proof is complete. 

We note that by setting m’ = n (hence n’ = m) in the given theorem, 
we obtain Takagi’s theorem immediately. 


IV. DISCUSSION 


Two examples to which Takagi’s theorem does not apply while our 
theorem does are shown in Figs. 6 and 7. In both figures, the channel 
graph (a) is superior to that of (b). The comparison in Fig. 6 is especially 
useful since the degrees of corresponding vertices are exactly the 
same. 

Next we ask can we generalize our theorem in the direction that the 
states of the edges between stage 1 and stage 2 can also be dependent. 
We conjecture the theorem will still be true but no proof is known yet. 
Certainly we cannot hope to imitate the proof given here. A ready 
counterexample to this approach is illustrated in Fig. 8. 

If the two states we assume are such that the first edge between the 
first two stages and the first edge between the last two stages are blocked, 
then the two channel graphs can be reduced to the channel graphs in 
Figs. 9a and b, respectively. Clearly, the channel graph in Fig. 9b is su- 
perior. But from Takagi’s theorem or our theorem, we know that the 
channel graph in Fig. 9a is superior to that in Fig. 9b. 
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A general computer-aided design method for use with electromag- 
netic devices such as ringers, relays, and solenoids is described. The 
method is demonstrated by applying tt to the design of polarized bell 
ringers. A lumped-element model with electrical, magnetic, and me- 
chanical portions is used in the analysis. First, interaction equations 
are derived using a Lagrangian formulation applied to a simple model. 
Second, the model ts refined by subdividing the iron members and in- 
cluding more leakage paths. An electrical circuit analysis program 
assembles the equations for the electromagnetic portion of this more 
complete model and produces a subroutine that solves these equations. . 
A computer program has been written to predict the effects of changing 
motor parameters. The versatility and usefulness of the design tech- 
nique has been demonstrated by applying it to the Bell System 
TRIMLINE® telephone set ringer to achieve major design improve- 
ments. 


I. INTRODUCTION 


In 1952 the Bell System introduced the standard 500-type telephone 
set which uses a “universal” two-gong ringer (C-type) to meet a wide 
variety of service conditions.! A universal ringer must: (1) have two 
voltage sensitivity modes that ensure adequate operation under worst- 
case conditions and provide protection against cross ringing on party 
lines, (ti) be electrically polarized to protect against bell tapping due to 
dialing transients, (iii) have high impedance (~8 kohms) at ringing 
voltages so that multiple ringers can be bridged across the line, (iv) have 
high impedance (~120 kohms) at voice frequencies to prevent speech 
signal attenuation, and (v) have two coil windings for use with multiparty 
ringing circuits.” 

Recently a project was undertaken to design a miniature single-gong 
ringer that would have the low cost, reliability, sensitivity, and sound 
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Fig. 1—Ringer with rocking-armature motor (gong and clamp plate removed). 


output of the C-type ringer but would be significantly smaller. The small 
single-gong P-type ringer commonly used in TRIMLINE® telephones 
and other new telephone sets did not optimally meet all of these objec- 
tives and previous attempts to meet them were unsuccessful. To achieve 
these objectives, effort was channeled toward a rocking armature type 
motor®4 that drives a single clapper and fits under the gong. A mathe- 
matical model was essential for the design process in order to perform 
parametric analyses and design optimization. The mathematical model 
had to realistically account for flux saturation of materials and flux 
leakage paths in the design, both of which have significant influence on 
the performance of a compact ringer motor. 


ll. RINGER OPERATION 


The structure of a ringer with a rocking armature motor is shown in 
Fig. 1 and an exploded view of the motor is shown in Fig. 2. In the absence 
of coil current a bias spring (not shown) returns the armature to the 
position shown in Fig. 1. The permanent magnet flux passes through 
both armature gaps and nearly saturates the shunt member. The ringer 
is electrically polarized since only coil current of one polarity will cause 
operation. Coil current in the operate direction increases the flux in gap 
A and decreases the flux in gap B causing clockwise armature rotation 
which drives the clapper away from the gong. As the current returns to 
zero the armature returns to its normal position and drives the clapper 
into the gong. For negative coil current the shunt member of the mag- 
netic circuit has a relatively low reluctance so that most of the coil flux 
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passes through the shunt instead of the working gaps. The low reluctance 
shunt results in the ringer having a high electrical inductance which is 
necessary for resonance of the ringer circuit at the frequency of the 
ringing power supply, 20 Hz. 

The ringer motor consists of a bobbin-wound coil, a pole piece as- 
sembly, a permanent magnet, a pivot pin, and an armature of low carbon 
steel. Figure 2 shows one version of the rocking armature motor which 
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Fig. 2—Rocking-armature motor. 


is described in this paper. In this version the pole piece assembly consists 
of two lamination stacks sandwiched between two silicon steel lamination 
straps and held together with three aluminum rivets. 

Basic motor performance is usually studied using measured static 
torque curves that show the relationship between blocked armature 
torque, armature position, and dc coil current. Figure 3 shows two ide- 
alized torque curves, one with the armature blocked in the nonoperate 
position (gap B closed) and the other with the armature blocked in the 
operate position (gap A closed). A single cycle of very slowly varying coil 
current is shown below the torque curves to illustrate the quasistatic 
operation of the motor. A major objective of the mathematical model 
described here was the prediction of static torque curves from basic 
motor parameters such as dimensions, material magnetic characteristics, 
and armature displacement. 
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Fig. 3.—Ringer motor static torque curves and quasistatic operation. 


lil. MATHEMATICAL MODEL 


Successful methods of modeling electromechanical devices such as 
relays and solenoids have appeared elsewhere.®»* The electromechanical 
ringer is a complex device and a complete analysis would involve the 
solution of Maxwell’s three-dimensional field equations in space occu- 
pied in part by nonlinear material with memory. Fortunately, because 
ringers are low-frequency and low-velocity devices, the following sim- 
plications can be made for design analysis: 


_ (t) The field problem can be approximated by a network where each 
volume of space occupied by a uniform material is represented as one 

or more lumped elements. Each element represents a physical effect in 
the uniform member, such as reluctance or loss. 

(it) All magnetic members except the permanent magnet are as- 
sumed not to have memory (no hysteresis). 

(iit) Stray flux paths through air can be adequately fepiesented by 
lumped leakage reluctances. 


A lumped-element model for the magnetic portion of the motor of Fig. 2 
appears in Fig. 4. This model uses the common flux-current analogy 
between magnetic circuits and electric circuits. Iron members are 
modeled by nonlinear reluctances and are characterized by normal B-H 
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Fig. 4—Lumped element model (equivalent circuit) for magnetic portion of rocking- 
armature motor. 


curves. The magnet is modeled using its demagnetization (de-mag) curve 
and recoil permeability. A function was fitted to reference data for the 
B-H characteristics. The fitting process and results are described in 
Appendix A. 

Leakage flux paths are modeled by lumped reluctances. The number 
and magnitude of these reluctances were determined from available 
formulae and prototype measurements as described in Appendix B. 


IV. EQUATION FORMULATION 


A consistent way to develop the equations for a system that has me- 
chanical, magnetic, and electrical energy is to use the Lagrangian for- 
mulation. The procedure is illustrated in Appendix C by applying it to 
the simplified model of Fig. 5. The equilibrium equations for the sim- 
plified model are a set of five differential equations, one electrical, one 
mechanical, and three magnetic, the simultaneous solution of which 
describes the dynamic performance of the motor. 

The equilibrium equations are reduced to the static case since only 
quasistatic motor performance is considered here. In the resulting set 
of nonlinear equations the current source becomes a magnetomotive 
source NI in one of the magnetic equations and the mechanical equation 
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Fig. 5—Simplified lumped element model for rocking-armature motor. 


gives the static torque on the armature. See Appendix C, eq. (21). Ar- 
mature torque t, is given by 


ob: OP b3 oP, (ba + bp)? OPng 
2P? 00 =2P? 08 oP 7 O60 


where the ¢’s and P’s are the gap fluxes and permeances respectively. 
Magnetic permeance is the reciprocal of reluctance. 

It is apparent from this equation that for accurate calculations of 
torque, appropriate functions for gap permeance must be determined. 
This is done in Appendix D. The gap fluxes are found by solving the three 
magnetic equations for specified values of coil current and armature 
displacement. The gap fluxes are then used in the torque equation above. 
To achieve accurate calculations of circuit fluxes it was found necessary 
to include up to 30 elements in the magnetic portion of the model. 

Having used the Lagrangian formulation to determine the terms of 
mechanical-magnetic interaction for the simplified model, circuit 
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analysis techniques were then used to generate the equations for the 
electromagnetic portion of the more complex model shown in Fig. 4. This 
step permitted increasing model complexity to account for significant 
measured effects without tedious derivation of new equilibrium equa- 
tions. 


V. SOLUTION OF EQUATIONS 


The next task is to find the solution to a set of simultaneous nonlinear 
equations. Since linear equations are easy to solve, at least in principle, 
the set of nonlinear equations may be solved most easily by finding a 
sequence of solutions to related linear equations that converges to the 
solution of the nonlinear equations. In general, closed-form (exact) so- 
lutions do not exist for nonlinear systems, so some such iterative method 
of solution must be used. 

The procedure adopted is well known as Newton’s method. The 
nonlinear system is expanded in a Taylor series about some trial solution. 
Retaining only the linear terms in this expansion, the resulting set of 
linearized equations is solved to yield a new approximation to the solu- 
tion. Successive solutions converge to the solution of the original non- 
linear system. 

The magnetic portion of the more complex model was analyzed using 
the Circuit Analysis Program for Efficient Computer Optimized Design 
(CAPE COD).’ From a topological input description of the magnetic cir- 
cuit, CAPE COD produces a subroutine that solves the network equations. 
This routine is combined with the interaction terms derived from the 
Lagrangian formulation to provide the complete device model for the 
static case. ; 


5.1 Sequence of analyses 


The first calculation procedure determines a maximum operating 
point for the magnet. For this calculation, the magnet is described by 
its de-mag B-H curve, and the independent coil current source is set to 
zero. A tentative operating point is found by iteration, as described 
above..Once this maximum operating point is determined, the operating 
point of the magnet is lowered until the ringer will just begin to operate 
at a specified current. The results of this analysis are then printed. 

A graph for the de-mag B-H curve for the magnet and the “load line” 
of the rest of the circuit as seen by the magnet is produced. For this 
calculation of the load line, a value for H is selected, and the slope of the 
B-H curve for the magnet is set to zero. The solution of this system gives 
the corresponding value of B on the load line for the selected H. The 
independent source associated with the coil is again set to zero for these 
calculations. 
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Fig. 6—Typical program output showing the predicted effect of decreasing the cross- 
sectional area of a magnetic shunt member. 


To calculate torques, the magnet characteristic used is that of the 
recoil line through the final operating point. The independent source 
associated with the coil is set to the appropriate ampere-turns, and the 
calculations proceed as before. The torque on the armature is calculated 
from the flux through the gaps at the armature. The displacement of the 
armature and the ampere-turns of the source are varied to produce a 
family of curves. Figure 6 is a sample of the plotted program output 
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Fig. 7—F lux versus coil current, predicted and measured. 


showing both the permanent magnet and torque characteristics of a 
motor and the predicted effect of changing a single parameter, the 
cross-sectional area of the shunt. 


VI. VERIFICATION OF MATHEMATICAL MODEL 


Measurements were made on an early prototype of the rocking ar- 
mature motor in order to verify the mathematical model. The solid 
curves of Fig. 7 show the predicted flux versus coil current in key mag- 
netic members. Measured fluxes, some of which are indicated, showed 
good agreement with predicted values. Agreement was particularly good 
for fluxes ¢g and ¢», which determine armature torque. 
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Fig. 8—Predicted and measured torque curves. 


The comparison of predicted to measured armature torque is shown 
in Fig. 8. Agreement between these two curves is very satisfactory. A 
comparison of the predicted to measured effect of changing two im- 
portant variables, armature displacement and magnet strength, is shown 
in Fig. 9. The measured change is considered to show adequate agree- 
ment, especially in the important midrange currents where ringer bias 
adjustments are made. 


Vil. OPTIMIZATION 

The mathematical model of the ringer motor aided first in identifying 
key parameters and second in finding a set of parameter values which 
optimize motor performance. The General Purpose Optimization 
Package (GPOP)’ was combined with the ringer analysis program to 
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Fig. 9—Predicted and measured effects of parameter changes. 


perform that optimization. First a criterion function based on the desired 
motor characteristics was developed. Second, the program determined 
a set of values for the parameters that improves the performance of the 
ringer relative to a set of initial values provided to the program. 


7.1 Desired motor output 


Basic motor performance is evaluated through use of static torque 
curves. The two idealized torque curves of Fig. 3 are shown in Fig. 10, 
with important points indicated. The two “stick” torques are defined 
to occur at zero current with the armature on either pole. They are in- 
dicated on the figure by the labels T, and Ty. Torque T, must be either 
negative (counterclockwise) or, if positive (clockwise), must be small 
enough to be overcome by a practical bias spring. This ensures that the 
armature always returns to its nonoperate position in the absence of coil 
current. 

Torque 7, must be large enough to drive the clapper against the gong 
with sufficient velocity. At the same time Ty must not be so large that 
excessive armature impact noise and wear are produced by the armature 
striking pole B. 

Of primary importance is the slope of the nonoperated torque curve 
immediately above turn-on current /;. This is called the motor torque 
factor and should be as high as possible to facilitate ringer sensitivity 
adjustments and to maximize ringer reliability. The slope between J2 
and /3 is usually lower than that between J, and J2 but must be high 
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enough to prevent reliability problems in high ringing-voltage situations 
with the ringer bias spring in its high-tension position. 

Another important criterion is the inductance of the ringer. The 
change in flux through the coil between peak currents was taken as an 
approximation to the inductance. 
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Fig. 10—Idealized torque curves. 


7.2 Criterion function 


From these diverse requirements a function is formed that expresses 
in a single number a measure of performance for a ringer built to match 
any given set of parameters. The form of the function selected involves 
the use of penalty functions. The idea is fairly simple. To a function 
describing how good the performance is, functions are added that are 
zero if a given constraint is satisfied and nonzero if the constraint is 
violated. Appendix E shows the function used. 


Vill. RESULTS 


8.1 Rocking armature motor 


The modeling technique was used to evaluate variations of the basic 
rocking armature motor so as to best meet the diverse requirements. 
Based on the theoretical results obtained, laboratory models of ringer 
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designs were built and evaluated in terms of the Bell System ringer re- 
quirements. Figure 11 shows a motor design which meets requirements 
and appears simple to manufacture. Iron straps are used to sandwich 
the lamination stacks and also to reduce flux densities in the assembly. 
A short Alnico 8 magnet is used which allows the use of a straight ar- 
mature. Figure 12 shows the predicted and measured results of an op- 
timization run. As shown, a significant reduction in pole A stick torque 





Fig. 11—Ringer motor with straight armature. 


is achieved, but with a sacrifice in torque factor. Further evaluation and 
refinement of this rocking armature design was not pursued since effort 
was turned toward improving the P-type ringer as described below. 


8.2 Bell System P-type ringer 


The analysis techniques developed for the rocking armature motor 
were applied to the Bell System small P-type ringer. The objective was 
to investigate potential improvements to the torque curve with simple 
design changes. Figure 13 shows the motor structure which is a single- 
ended rocking-armature type with air gaps on opposite sides of the ar- 
mature. Figure 14 shows the equivalent circuit used to analyze it. It was 
predicted and then verified with measurements that a significant im- 
provement im torque factor is achieved with two basic design changes. 
First, the magnetic bias force is reduced by moving the pivot pin closer 
to the magnet center. This results in higher magnet flux after the mag- 
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Fig. 12—Predicted and measured effect of optimized parameters. 


netization adjustment. Secondly, the flux density at a critical point of 
saturation is reduced by increasing the cross-sectional area of the pole 
piece. The measured results are shown in Figure 15. 


IV. CONCLUSIONS 

A general method of computer-aided design for electromagnetic de- 
vices has been presented. The method was applied to the design of two 
polarized bell ringer motors. A mathematical model and computer 
simulation of the rocking armature motor made it possible to accurately 
and efficiently study motor performance as a function of the many pa- 
rameters involved. The simulation methods were also applied to the 
P-type ringer commonly used in the TRIMLINE telephone and resulted 
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Fig. 13—-P-type ringer motor. 
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Fig. 14—Lumped element model of P-ringer motor. 


in a significantly higher torque factor. This improvement in addition 
to others not covered here have led to a modified P ringer virtually 
meeting the original goals of the rocking armature ringer develop- 
ment. 
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Fig. 15—P-type ringer motor measured torque curves. 


APPENDIX A 
lron member characterization 


The characteristics of an iron member are modeled using the element’s 
length, width, and the normal B-H curve for the material. Several 
methods have been used for representing this curve: a table of (B,H) 
pairs may be given, simple hyperbolas may be fit to a range of the curve, 
a series of exponentials may be fit to the curve, a function of a polynomial 
divided by another polynomial (rational function) may be fit to the data. 
For this analysis H is given as the ratio of two polynomials in B. 

A family of functions with numerator and denominator polynomials 
of varying orders was fitted to the data. The function selected from this 
family provided a good combination of low order, close approximation 
to the data being fitted, and smooth behavior of the derivative. Figure 
16 shows the selected curve fit and the data points used in the fitting 
process. 


APPENDIX B 
Leakage elements 

Leakage elements are reluctances which represent stray flux paths 
through air from two points in the magnetic circuit. Leakage elements 
placed across nonlinear magnetic elements reduce the sharpness with 
which the elements go into saturation and are essential to obtaining 
agreement between predicted and measured fluxes. 
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ngB° +484 +383 + n2B2 +48 +19 


H=B* 
m4B4 +m3B? + m2B2 +m,B +1 


Ng = 10.97018 my = 11.34651 
ny = 31.94147 mp = —22.93222 
Ny =-97.93912 m3 = 14.37377 
n3 = 88.21910 mg = —2.94200 
1g = —35.45198 

ng = 5.72501 


8B, TESLA 


B, TESLA 


O FITTED DATA 





0.001 0.01 0.1 1.0 
H,A/mm 
Fig. 16—BH curve fit for silicon steel. 


The number of leakage elements and their placement in the model 
can be observed by sprinkling iron fillings on a prototype as shown in 
Fig. 17. Precise numerical values of these leakage elements are difficult 
to obtain. However, the method of “estimating the permeances of 
probable flux paths” from Ref. 8 has proven useful in obtaining first- 
order estimates. Also, low-frequency inductance measurements of partial 
motor assemblies can be used to estimate leakage around the coil. Flux 
measurements on the physical model are then used to refine the initial 
estimates. Once determined, the leakage reluctances remain constant 
unless major changes in geometry are made. 


APPENDIX C 
Lagrangian formulation 


The first step for the Lagrangian formulation is to select an appro- 
priate set of generalized coordinates which define the state of the system. 
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Typically these are flux linkage, node voltage, capacitor charge, inductor 
currents, and physical displacements. For this ringer analysis, coordi- 
nates representing linkage, node voltage, and angular displacement of 
the armature where selected. 

The second step of the procedure requires formulation of the Lag- 
rangian function and a Rayleigh dissipation function in terms of the 
generalized coordinates. Since most power dissipation occurs in winding 
resistance and since the technique becomes burdensome with individual 
hysteresis and eddy current loss elements, a single electrical loss element 
was assumed to account for the total system dissipation. 
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Fig. 17—Prototype ringer with iron filings. 
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The third step applies Lagrange’s formula to the energy functions to 
produce a set of simultaneous differential equations. Primarily this in- 
volves taking partial derivatives and collecting terms. 

The advantage of this technique is that it systematically accounts for 
the interaction between the magnetic and mechanical portions of the 
system. The disadvantages are as follows: first, the choice of generalized 
coordinates is not necessarily obvious; second, hysteresis and other losses 
are difficult to fit into the formulation. 

The Lagrangian J is defined as the difference between the total sys- 
tem coenergy function 7’ and the total system energy function VV as 
defined below. Choosing a nodal formulation for the electrical portion 
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of the system and using the functional dependencies prescribed in Ref. 
9 we have: 


L(0,0,,\,t) = T’(6,0,d,t) — V(O,A,t) (1) 


where @ is the angular displacement of the armature, d is the flux linkage, 
t is time and dots above quantities denote differentiation with respect 
to time. 

The system coenergy function is the sum of the mechanical (kinetic) 
and electrical coenergy functions: 


T’(0,0,A,t) = T’(8,0,t) + W.(K,8) (2) 


The system energy function is the sum of mechanical (potential) and 
magnetic energy functions: 


V(8,r,t) = V(G,t) + Wi (A,9) (3) 
For the system of Fig. 5 the mechanical coenergy function is 
: | ee 
T’(0,t) = 5 oe (4) 


where J is the moment of inertia of the armature about the pivot axis. 
The electrical coenergy function is: 


Wi(\y) = 5 ON (5) 


where C is the value of the ringing capacitor. The total system coenergy 
function now becomes: 


: DY ies «all sate 
T'(O,A1,t) = role + A Ch? (6) 
Since no change in mechanical potential energy is assumed, it remains 
to find the magnetic energy function, W,,. The energy W in a lumped 


magnetic element of length @, area A, flux density B and magnetic field 
H can be approximated by the following integral: 


B 
W=2A f H(B’) dB’ (7) 
Since B = ¢/A, 
to) op’ b gp’ 

w= f tH (7) de’= f x (2) de’ 8 
: 7 p A @ (8) 

where €H is called the magnetomotive force (MMF) and is labeled #. 
For linear magnetic elements F is a linear function of flux, i.e., ¥ = 


R¢, where R is the reluctance of the element and defined as ¢/y,A. 
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Therefore the energy in a linear element is: 
¢ 1 
W = f Ro! do! = 5 Re? (9) 


In a system model that includes more magnetic elements than those 
which directly link the coil, it is convenient to use flux (@) for coordinates — 
in the magnetic portion of the circuit rather than flux linkage (A). Both 
have the same dimensions, differing only by the dimensionless quantity, 
number of turns. 

The total magnetic energy function, W,,, can now be expressed by 
summing the energy of magnetic elements: 


1] 1 1 
Wr (Oin8) = = Ra (0% + 5 Ro(0)48 + Re (Or 
1 dm ; ‘ d ; ; 
+5 Rada + f Ab \ida + f * Falbiy) doy 


ge ; P de , , 
+ f Fo($,) dd, + f, %.($.) dg, (10) 


Note that the first three terms are gap energies and are functions of ar- 
mature displacement. The last four terms are the energies in the non- 
linear elements. 

Continuity of flux in the magnetic circuit provides the following 
constraint equations: 


dm = ba + bb 
ba = ba — $e (11) 
he = be + hd 


After substituting the constraint equations into eq. (10) and then using 
eqs. (3), (6), and (10) in (1), the Lagrangian becomes: 


LOMAGit) = SIP + 5 CN = RalO)d2 
1 ‘ect a ; 
= 5 Rol) obe = 5 emg (9) (ba r db) ce 5 Pea (ba ae fc) 
dat bb ? 4 ga— be , , 
7 jf Fm (Gin) doin — § Fada) db 


ct , , de , , 
-f aaa Ne eee Sf Fel) 46. (12) 


The Rayleigh dissipation function D is 
Az _ (Ao = Ar)? _ (Noe = Ai)? 


We, 
os J, Ror QR QR QR (13) 
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where \, is the voltage across resistor R and the voltage across the coil 
dz is equal to Ng. 
Lagrange’s formula can now be applied: 


d sol oL . oD 
— (——-) --——-+——=4; 4 
dt cai) Og; 24; 
For the electrical coordinate \1, we have 
ee 1 ‘ 
Cry as R (Nod. = Ay =0 (15) 


For the mechanical coordinate @ we have 


2 2 2 
sh 4 2 Ra, $5 2R6 , (bat $0)? mg 
2 06 2 086 2 ol] 


For the magnetic coordinates ¢g, ¢,, and ¢, we have 


(Ra + Ring + Rra)dba + Ring dd = Reade 
5 F m(ha + bo) + Falda a de) =0 


Rmgba + (Ry + Rmg)bo + Fm (ba = ho) 


= tg(t) (16) 


+ Folge +oo)=90 (17) 
Fo(ge) + Fo(he + db) i Fala = de) _ Riada + Riadbe 


+o (Née — ha) = Ni(t) 


where Ra, Rmg, and R, are understood to be functions of 6. After the 
permanent magnet has been charged and then stabilized, its second 
quadrant minor loop operation is modeled by a source of constant. 
magnetomotive force, —F;, in series with an internal magnet reluctance 
Rm. Thus, 


Fm (ba + do) = =f + Rn (ba + ho) (18) 


The five equilibrium equations are then reduced to the static case by 
setting all time derivatives to zero and representing constant coil current 
and armature torque by J and ¢, respectively: 


(Ra + Ring + Ria + Rm) ba + (Rng + Rm) bd 
— Reabe + Falga — be) = Fs 
(Ring + Rm)ba + (Ro + Rng + Rn)bo + Felbe + 5) = Fs (19) 
—Rraga + Riabe + Folbe) + Felbe + bb) — Falda — bc) = NI 
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The result is three nonlinear algebraic magnetic equations (19) which 
must be solved for the three unknown fluxes. Fluxes ¢, and ¢, are then 
used in the mechanical eq. (20) for calculating torque developed on the 
armature. For convenience in handling gap functions, the torque equa- 
tion is rewritten in terms of gap permeances (reciprocal of reluc- 
tance): 

$a OPa $5 OPy (ba t+ $b)? OPmg _ 


— fa Soa fb Sb Na Obl Ome = 21 
2P? 00 2P2a00 2P2, av 21) 


APPENDIX D 
Working gap functions 


Since armature torque is determined from gap permeances and the 
partial derivatives of these permeances with respect to armature dis- 
placement, accurate mathematical expressions (gap functions) for these 
permeances must be derived. Gap functions which account for both 
fringing flux and main gap flux between nonparallel surfaces can be 
derived by use of the principles given by Roters.® By introducing the 
“magnetic pivot” concept,!° these functions are modified to account for 
the fact that the mechanical pivot is not located in the plane of the gap 
magnetic surfaces. 

The derivation of the permeance of gap B (P,) will serve to illustrate 
the method. The region between the armature and pole B is divided into 
11 simple geometric elements as shown in Fig. 18. Each element repre- 
sents a region of flux paths. Total gap permeance is the sum of the in- 
dividual permeances, each defined in terms of gap geometry as shown 
below: 


P, (0) = permeance of main gap 


oF» (FR) rp 0 


6 6B, + Ling 
P,(0) = (22) 
Hol (pr a ear) for |0] «< = 
6Bo+Lmg 0B, + Ling By 


P, = permeance of one-half of a semicircular cylinder with length Fy, 
and radius 0B, + Ling 


Po = 0.52 tos (23) 


P3(0) = permeance of a quarter annulus with axial length F,, and 
bounding radii of 6B, + Lmg and 0B, + Lng +t 


2UoFr 


P3(6) = if (1 + nae (24) 
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Fig. 18—F lux paths at gap B. 


P4 = permeance of a semicircular cylinder with length F, and diameter 
OBo + Ling 


P4 = 0.26 uoFs (25) 


Ps (9) = permeance of a half annulus with axial length F;, and bounding 
radii of 1h(0Bo °F Ling) and ',(0Bo + Lg) +L 





Fo 2t 
P;(6) = He= 2] ia) 9 
Se OBo + Ling (26) 


Pg = permeance of a semicircular cylinder with length By — B, and av- 
erage diameter 6(By + B,)/2 + Ling 


Pg = 0.26 wo (Be — By) (27) 


P7 (6) = permeance of a half annulus with length Bz — B, and average 
bounding radii of . 


5 (2) +a 


2 2 2 
and 
5 (2S) bmg yy 
2 2 2 
»(Bo —-B Ot 
P,(6) = #2 2 Wy ae (28) 
VT 
aaa *) + Lng 


COMPUTER-AIDED MAGNETIC CIRCUIT DESIGN 201 


P3(6) = permeance of a spherical quadrant with diameter 0Bz + Ling 
P3(0) = 0.077 Uo (8Bo + Lng) (29) 


Py = permeance of a quadrant of a spherical shell with bounding radii 
of V,(0Bo + Lae) and 1,(0Bo + Line) +t 


Hot 


Po» = ae (30) 
P40(9) = permeance of a spherical quadrant with diameter 0B, + Ling 
P0(8) = 0.077 Ko (6B, + Lmg) (31) 


P ; = permeance of a quadrant of a spherical shell with bounding radii 
of (0B, I Ling) and 44(6B, + Ling) +t 


Mot 
4 
Now we can write P, (0), the total permeance of gap B: 


+ 2P8(0) + 2P9 + 2P30(0) + 2P31, (83) 


Since Py = 2P4 and Pg = Py, 


P, (0) = P,(0) + P3(0) + 8P4 + Ps(0) + 2P.6 + 2P1(0) 
+ 2P3(0) + 2P10(0) + 4Py, (34) 





Py = (32) 


The total permeance of gap B (34 and 22 to 32) has now been expressed 
in terms of ringer motor parameters. The permeance is then differen- 
tiated with respect to 0 for use in the torque eq. (21). Gap A and the 
magnet-to-armature gap are handled in the same way. 


APPENDIX E 
Criterion function 

The first term in the function to be minimized is the reciprocal of the 
torque factor squared. The rest of the terms are constraints imposed by 


means of penalty functions. Refer to Fig. 10 and the definitions in Section 
VII. Let 


R, = 1/(torque factor)? 


Rg =T, 
R3 = Ty 
Ra= 7; 
Rs =T3- T2 


202 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1978 


Rg = estimate of inductance 


Ry 


thickness of shim in gap A 


(torque factor from [2 to I3) 


Rs = roll = 
BrOsnege wen (torque factor from J, to I») 


Then let 
D;= 0, Ring S Rp = Riaz 
D; = Ri — Rin, Ri < Rin; 
Dj = Rimax; — Ri, Ri > Rinax; 
ne criterion function P is formed by taking the weighted sum of the 
P= = A;D;? 


The designer controls the optimization process by judicially choosing 
the weighting factors A; and the constraints Rin; and R max,. 
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The Preparation of Optical Waveguide Preforms 
by Plasma Deposition 


By R. E. JAEGER, J. B. MACCHESNEY, 
and T. J. MILLER 
(Manuscript submitted February 11, 1977) 


Optical fiber preforms are prepared by a technique similar to the 
modified CVD process, except that an RF plasma is used. Optical losses 
as low as 6.5 dB/km at 1.06 um were achieved. Higher reaction ef- 
ficiencies resulted in deposition rates more than three times greater 
than the modified CVD process. 


Il. INTRODUCTION 


The modified chemical vapor deposition (MCVD) process! has achieved 
acceptance as a means of preparing low-loss optical waveguides. A recent 
study of the preparation of silica-clad germania borosilicate optical fibers 
by this process? indicates that it can provide low-loss optical fibers having 
properties which surpass minimum requirements for wide-band com- 
munications applications. Typical fibers produced by MCVD are char- 
acterized by low loss (<4 dB/km at 1.06 um), a reproducible graded index 
profile (a ~ 2.0) and controlled core/cladding ratio, circularity and 
concentricity. 

Although current deposition rates are competitive with other pro- 
cesses, higher rates are desirable. Furthermore, the efficiency, defined 
as the ratio of glass deposited/reactants in, is generally below 50 percent 
for most GeO2 containing compositions. These considerations stimulated 
our investigation into using the ionized environment of a plasma to in- 
crease the rate and efficiency of the MCVD process. Initially both mi- 
crowave and RF plasmas were considered. In the former, ionic species 
are produced which deposit as a vitreous deposit directly on the tube 
walls as the result of a heterogeneous reaction similar to conventional 
CvpD.? This advantage, however, is purchased by the use of reduced 
pressures which dictates a low flux of reactants. As a result, even at 100 
per cent efficiency, the deposition rates in the microwave method need 
be no higher than MCVD and frequently are lower. 
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Fig. 1—Apparatus for preparing preforms. 


In contrast, the use of an RF argon-oxygen plasma operating at a 
pressure of 1 atmosphere provides the potential for higher deposition 
rates as well as increased efficiency. In this type of plasma, temperatures 
in the vicinity of 10,000-20,000°K exist, producing an intense reaction 
region in which highly ionized species occur. Under these conditions, 
reactions between ionized species occur resulting in a deposit on the tube 
wall which is similar to that produced by MCVD. Significant improve- 
ments in the efficiency and deposition rate achieved by the RF plasma 
MCVD process are derived from the ability to produce gas phase reactions 
in a highly ionic environment. 


ll. EXPERIMENTAL 


The apparatus consists of a glass working lathe positioned adjacent 
to an RF generator as shown in Fig. 1. An RF coil, mounted to a motorized 
carriage on the lathe was connected to the generator via water-cooled 
flexible leads. In a typical deposition, a 25 X 28 mm TO8 silica tube was 
positioned in the lathe and a flow of argon was started. All gas flows were 
left to right and were controlled with rotometers. The plasma was ini- 
tiated at a frequency of 4.5 MHz by striking an arc with a graphite rod 
inserted into the tube in the vicinity of the coil. The coil was traversed 
at a constant rate over ~25 cm length and the tube was rotated as shown. 
Gaseous SiCl, and GeCl, were introduced to the reaction zone using an 
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Fig. 2—Loss spectrum. 


oxygen carrier which was bubbled through room temperature chloride 
liquids. BCls gas was directly introduced from a cylinder of liquid BCls. 
SiC], at 92 cc/min and BCls at 10 cc/min were initially introduced into 
the gas stream to produce a borosilicate “barrier layer.” Deposition was 
carried out for 5 minutes before GeCl, at 65 cc/min was added to produce 
a germanium borosilicate core. Core deposition continued for 15 minutes, 
at which time the RF power and all gas flows were turned off and the tube 
collapsed using an oxyhydrogen burner. 


lil. RESULTS AND DISCUSSION 


It is well known that the stability of an RF plasma may be related to 
the diameter of the tube in which it is being sustained. We found that 
under our operating parameters, stable operation could be obtained at 
a frequency of 4.5 MHz by using 25 mm ID tubes. Fundamental consid- 
erations indicate that, in general, it is desirable that the tube radius be 
greater than the skin depth‘ of the plasma, which is proportional to the 
1) power of the ratio of the plasma resistivity to the operating frequency. 
Hence, higher frequencies would be one variation that would simplify 
the use of smaller tubes. 

Figure 2 gives the loss spectrum for a 1100 meter length of a plasma- 
deposited fiber with a 35 um core and a 110 wm OD having a minimum 
loss of 6.5 dB/km at 1.06 um. The peaks are due to overtones and com- 
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Fig. 3—Unetched microstructure. 


binations of the fundamental OH absorption band and the loss at 0.95 
um indicates a water content of ~50 ppm.° Since typical MCVD fibers 
exhibit a minimum loss below 4 dB/km and contain only a few ppm 
water, an attempt was made to determine the extent to which the current 
plasma process parameters contributed to the excess water content and 
loss, by preparing a preform in the conventional manner using an ox- 
yhydrogen torch as the heat source. This resulted in a fiber having the 
same OH content but with a minimum loss < 3 dB/km at 1.06 um, 
suggesting that the excess loss at this wavelength is related to specific 
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process parameters while the OH concentration may be due to the gen- 
eral lack of humidity control in the surrounding environment. 

The unetched microstructure of a core deposit for a 200 wm fiber 
viewed in an optical microscope, shown in Fig. 3, indicates some of the 
limitations encountered under the present operating conditions. The 
noncircularity of the core reflects the difficulty in collapsing a 25 mm 
tube while the rings result from compositional variations in the deposit. 
Examination of an etched core deposit in a scanning electron microscope 
revealed that the dark rings are GeO deficient regions. Thus, the excess 
loss in this fiber may result from scattering caused by inhomogeneity 
of the deposit and/or compositional variations in the core. In either case, 
large excursions of the tube wall temperature during the deposition are 
the probable cause. However, other loss phenomena have not been ruled 
out. 

In order to make very low loss preforms, it may be desirable to heat 
the deposit on the wall to a temperature where fining and homogeniza- 
tion of the glass take place with minimal vaporization of the more volatile 
constituents or heating the tube to temperatures where it will deform. 
Temperature control in the plasma deposition process was attained 
primarily by traversing the plasma ball along the tube at a rate sufficient 
to accomplish these objectives. Additional cooling of the outside of the 
tube to prevent deformation was provided by an air stream directed at 
the region surrounding the plasma. However, other important param- 
eters that affect the wall temperature include the composition and ve- 
locity of the plasma gases and the RF power levels. For example, in- 
creasing the oxygen/argon ratio in the plasma requires increased RF 
power, both of which produce a high plasma temperature that is trans- 
ferred to the walls. The evaluation and control oh these parameters is 
currently underway. 

The principal advantage of plasma deposition is the potential for in- 
creased deposition rates and improved overall reaction efficiency. For 
approximately equivalent reactant flows a rate of deposition of 0.19 
g/min was obtained compared to a rate for MCVD of 0.05 g/min. This 
increase in the deposition rate supports the belief that the overall reac- 
tion efficiency is markedly increased for the RF plasma process. In ad- 
dition, the NA of the plasma-deposited fibers were 0.28 compared to 0.24 
for the conventional process, indicating increased efficiency for GeQ2 
incorporation. Furthermore, reasonable extension of the results achieved 
to date suggest that deposition rates perhaps an order of magnitude 
greater than normally achieved by the MCVD process may be obtain- 
able. 

IV. CONCLUSIONS 

The RF plasma CVD process is capable of producing low loss (< 10 

dB/km) optical fibers at high deposition rates. A better understanding 
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of the effect of various operating parameters on the nature of the depo- 
sition process is necessary to exploit the full potential of this process. 
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ERRATA 


In the Bell System Technical Journal, 56, No. 8 (October 1977), Figs. 
2 and 3 on pp. 1454 and 1455 should appear as follows. 
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Fig. 2—(a) Diagram of a grating resonator in the external transmission configuration. 
(b) The transmission spectrum at ~145.5 MHz for a resonator on YZ-LiNbOs using Ti- 
diffused gratings with A = 12 um and L = 6.48 mm. (c) The calculated transmission 
spectrum for the device in (b) using x = 3.74 cm~! and a = 0.036 em7!. 
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Fig.3—(a) Diagram of cascade of N identical resonators. (b) The transmission spectrum 


at 145.5 MHz of a cascade of three identical resonators with A = 12 umand L = 3.84 mm. 
(c) The calculated transmission spectrum for the device in (b) using « = 3.55 cm=! and 


a = 0.027 em—}, 


Equation (15b) on p. 1509 should be: 
2 _ M2 
Oomn = Yow?CinSn _ —- Ey se Yow? CF n (15b) 
M m=] 
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